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As an oxidizer for liquid fuels, 
Nitrogen Tetroxide has even more 
to recommend it than its high per- 
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ease with which it can be stored 
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advantage. 
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How to Test a Solid Rocket Motor 


In the early stages of research 
for new or improved propellants 
incongruous combinations of mate- 
rials are often brought together in the 
chemistry laboratory—in this case, a 
common drinking straw, a strand of 
copper wire, a small batch of propel- 
lant and a sewing needle. The work 
shown here is a step in an experiment 
to test burning characteristics of a 
new high-energy propellant formula. 
The propellant, having been extruded 
into the straw, is threaded with a 
strand of copper wire. Later the 
“motor” is hooked into a specially 
designed apparatus that ignites the 
propellant and records and evaluates 
its performance. 


From the igniting of grams of fuel 
in the chemistry labs to the test-firing 
of huge motors on the 2,000,000 Ib. 
thrust test stands at THIOKOL’s Utah 
facilities, research paces and supports 
every phase of the design, develop- 
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ment and production of a rocket 
powerplant. 

Through such research is coming 
the propellant formula that will give 
increased range to our ICBM mis- 
siles, and the tremendous power 
needed for spatial exploration. 

In its Rocket and Chemical Divi- 
sions, THIOKOL boasts one of the 
finest scientific research teams in 
rocketry—an expanding group, con- 
tinually in search of new members. 

For qualified chemists, there exists 
stimulating work in polymer re- 
search, in propellant analysis and 
formulation, in high vacuum tech- 
niques, and through the entire spec- 
trum of high energy fuel. 

Doors are open, too, to research 
engineers and scientists with interest 
or experience in: 

Solid State Physics, High Tempera- 
ture Gas Dynamics, Electronics of 
Rocket Exhaust, Electromagnetics, 


Thermodynamics, Magnetohydrody- 
namics, Astronautical Systems 
Analysis, Ion and Plasma Propul- 
sion, Nucleonics, and other phases 
of advanced propulsion work. 
You’ll work in areas that are good 
places to live—Denville, New Jersey; 
Huntsville, Alabama; Marshall, 
Texas; Moss Point, Mississippi; Brig- 
ham City, Utah; Trenton, New Jersey; 
Bristol, Pennsylvania. When you 
think of opportunity for yourself and 
your family, think of THIOKOL. Write 
to: Personnel Director at any plant 
address listed above. 


Thiokol 


is Research to the Core 


THIOKOL CHEMICAL CORPORATION 
Bristol, Pennsylvania 


®Registered trademark of the Thiokol Chemical Corp. for 
its rocket propellants, tiquid polymers, plasticizers and 
other chemical products. 
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Scientists have long been preoccupied with the techno- 
logical problems of Man and the Machine. The increas- 
ingly complex nature of advanced systems has created 
an urgent need to enhance man’s contribution to effective 
systems performance. The complicated nature of this 
relationship requires the skills of psychologists, social 
scientists, mathematicians, and engineers. 


At Ramo-Wooldridge, human engineering, personnel 
selection, individual and system training, display design, 
and communications are successfully integrated into 
systems design and development by the technique of 


a division of 


THE HUMAN FACTOR in 


large-scale simulation. cerning these opportunities write to Mr. D. L. Pyke. 


RAMO-WOOLDRIDGE 


P. O. BOX 90534, AIRPORT STATION « LOS ANGELES 45, CALIFORNIA 


Thompson 


L 


Simulated inputs enable scientists to observe a sdaiiaas as 
it operates in a controlled environment and make possible 
the collection of data on performance, training, human 
engineering, maintenance, and logistics and support. 
Scientists and engineers use this data to assure the design, 
production, and delivery of a unified system capable of 
high performance and reliability. 


Expanding programs at Ramo-Wooldridge in the broad 
areas of electronic systems technology, computers, and 
data processing have created outstanding opportunities 
for scientists and engineers. For further information con- 
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) THE past few years the impetus of the ballistic missile 
programs has led to tremendous progress in the ablation 
field. Successful re-entry and recovery of an ablating nose 
cone for the IRBM 1500-mile range was demonstrated by the 
Army Ballistic Missile Agency in August 1957. Then in 
April 1959 an ablating nose cone, after traveling the more dif- 
ficult 5000-mile ICBM range, was successfully recovered by 
the Air Force Ballistic Missile Division. The degree of suc- 
cess of these and other re-entry tests is, of course, classified in- 
formation, though we must conclude that ablation has been 
proved and that it represents a major advance over the heat 
sink as a means of absorbing aerodynamic heating. The use 
of ablating heat shields for satellite re-entry, rocket nozzles 
and other flight applications can be anticipated in the near 
future. 

It is unfortunate that flight results are not available for this 
article, since the real proof of any material comes only from 
tests under true conditions. A fundamental understanding of 
ablation can nevertheless be obtained by means of laboratory 
experiments combined with theory. The attempt here is to 
set forth the current understanding of ablation by drawing 
upon the available experimental and theoretical results in the 
literature. 

A literature survey shows that many materials have been 
subjected to a large variety of heating tests. Experimental 
facilities include solar furnaces (1),! oxyacetylene torches (1), 
rocket exhausts (2,10), water plasma jets (3 to 5), air plasma 
jets (6 to 9) and other high temperature sources. Data ob- 
tained from solar furnaces and water jets have furnished useful 
information on physical properties and thermal stress be- 
havior of materials; however, because of questionable simula- 
tion or unknown test conditions, it has been difficult to in- 
terpret the meaning of the ablation data. 

As will be seen later, the free stream stagnation enthalpy is 
one of the most important parameters to be duplicated in a 
laboratory test. The rocket exhaust produces gas enthalpies 
appropriate for flight velocities up to about 10,000 fps, but it 
suffers from improper chemical composition for simulating 
flight conditions. This facility has nevertheless yielded much 
valuable ablation data, and the Army Ballistic Missile Agency 
(10) successfully utilized rocket exhausts for the development 
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Reynolds numbers high enough for turbulent flow. Other hot 
air facilities, such as a pebble bed heater (12), produce stream 
energies corresponding to flight velocities up to about 7000 
fps, although no ablation data were available from such 
facilities. 

Available data indicate that two classes of materials have 
received the most attention: Pure plastics and glasses. Many 
composites also contain one or both of these materials as 
major components. Metals, metallic oxides and graphite by 
comparison are, in general, inferior, because they have either 
poor thermal stress resistance, a low heat absorbing capability 
(due to melting and flowing without evaporation) or a high 
thermal conductivity. It is possible that some of these latter 
materials may prove useful if fabricated in composite with a 
material of lower thermal conductivity, for example with a 
plastic. 

In the following sections of this report, brief discussions will 
be presented with the objectives of defining the parameters 
and physical properties which are important for ablating ma- 
terials, and reviewing various ablation theories and com- 
paring them with experimental results for specific materials 
where possible. Two applications will be briefly considered, 
namely satellite re-entry and rocket nozzles. - rs 


Definition of Some Important Ablation — 7 


Parameters 

To gain a feeling for some of the ablation parameters and 
material physical properties of importance, it is instructive to 
treat a simple example. Consider the case of a step function 
heat pulse go applied to a material for a given period of time. 
For the present we do not distinguish between laminar and 
turbulent heating or between locations along a particular 
surface. 

Initially, all of the incident heating is conducted into the 
material until the surface reaches the ablation temperature, 
corresponding to a temperature rise 7,. Ablation starts after 
this heatup period and then follows a second transient period? 
during which a steady-state ablation velocity is reached. In 
the steady-state condition, the temperature distribution is ex- 
ponential and is given by 


of nose cone materials. T = = 
More recently, the arc-heated air jet (11) has proved capa- where = 7 

ble of producing gas enthalpies from below the IRBM energy Vw = ablation velocity : 

level up to the ICBM and satellite energy levels. Although - -y = distance beneath the surface i ¥ 

the air are appears to be the best facility at present for re- _ a thermal diffusivity of the material > i ee 

entry tests, it has not yet been developed to operate at = dy = thermal thickness oi he 


Received July 31, 1959. 
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? This transient period equals the time for a particle of ma- 
terial to traverse a thermal thickness, i.e., £ = 57/V. 
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Utilizing this temperature distribution, the energy stored 
within the material per unit surface area is easily calculated 
tobe 


The energy absorbed by the ablation process is then the 
difference between the total heat input gt and Q,. If the 
material is efficiently utilized, Q, < got, and the ablated 
mass can be calculated by neglecting the energy stored. If we 
define the energy absorbed per unit mass of ablated material 
as Heer, then* 


Q = got ~ pVutHert 


or 


The fraction of the total energy stored by heat conduction is 


then given by 
[4] 


Thus we see an apparent conflict. A large Hes is required to 
minimize the mass loss by ablation, yet at the same time the 
energy, which is inefficiently absorbed by conduction, is in- 
creased in proportion to H.1:. Clearly, it is desirable to have 
a large Her and the lowest possible product of density, con- 
ductivity, ablation temperature and H.s. As an illustration, 
if Hes = 5000 Btu/lb, p = 100 lb/ft*’, T, = 5000 F, q = 
1000 Btu/ft? sec and ¢ = 25 sec, then in order that Q,./Q < 0.1, 
the thermal conductivity must be-less than 10~* Btu/ft sec F. 

The quantities which are difficult to obtain for any given 
material are the heat of ablation, the ablation temperature 
and the thermal conductivity. This paper is concerned pri- 
marily with the first two quantities, but a knowledge of the 
thermal conductivity is equally necessary. An important 
technique for measuring thermal conductivity was developed 
by Sutton (2). He employed thermocouples beneath the 
surface of materials under actual ablating conditions and de- 
duced the thermal diffusivity (and therefore thermal con- 
ductivity) from the measured temperature distributions. 

From the above discussion it is apparent that an ablating 
heat shield does two things: First, the major fraction of aero- 
dynamic heating is absorbed by the material which is ablated 
away, and second, some additional material, usually termed 
insulation, absorbs the stored or conduction energy Q,. De- 
pending upon the physical properties of a particular material, 
the heat transfer rate and the heating time, the insulation 
weight may be either more or less than the ablated weight. 
A material with the largest heat of ablation is not necessarily 
the best. Graphite is a good example to illustrate this point. 
Its heat of ablation is estimated at greater than 10,000 Btu/lb, 
yet because of a high thermal conductivity of 3 x 107? 
Btu/ft sec F, it acts as a heat sink in most re-entry applica- 
tions and will not absorb much of the heating by ablation. 
On the other hand, if graphite could in some way be fabricated 
in composite with a material of much lower thermal conduc- 
tivity, then its heat of ablation could be utilized. 

It should be emphasized that the combined insulation and 
ablation requirements must be considered together for a 
particular application. Strictly speaking, the ablation proc- 
ess is coupled with the internal heat conduction, and the com- 
bined problem is usually more complicated than the constant 
heat pulse case discussed above. However, if only a small 
fraction of the total energy is absorbed by conduction, as 
should be the case, then the ablation and conduction problems 


The definition of the effective heat of ablation H.; is the 
same in most of the literature; however several different symbols 
have been used. It should be noted that H.¢; is defined in terms 
of the heat transfer rate to a nonablating surface (at the ablation 
temperature). Throughout the present report, the heat transfer 
rate go will be regarded as known or prerequisite information. 
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can be decoupled. Generally, it is a satisfactory approxima- 
tion to determine the heat of ablation, the fraction of material 
vaporized and the surface temperature, assuming quasi- 
steady-state ablation for a given heating history. These abla- 
tion characteristics then furnish the boundary conditions 
necessary to calculate the internal heat conduction in a 
straightforward manner. An iteration procedure could also 
be carried out, if more accuracy were required. 

Inasmuch as the effective heat of ablation is prerequisite 
information for the design of an ablating heat shield, most 
of the basic papers in the literature have been concerned with 
the determination of this quantity for various materials. 
The next three sections of this report will be devoted to a re- 
view of the current understanding of the heat of ablation, and 
then I shall consider an example of a combined ablation and 
insulation problem. 

I have arbitrarily chosen to categorize ablation under three 
headings: Melting ablation, subliming ablation, and com- 
bined melting and sublimation. 


| 
Melting Ablation 


The interaction between a particular material and the gis 
boundary layer is, in general, a complicated process. The 
least interaction occurs for those materials which melt and 
flow without evaporation; consequently this class of materia!s 
is considered first. Though melting materials are probably 
the least interesting from a practical point of view, an under- 
standing of melting is fundamental to an overall understand- 
ing of the ablation problem. 

Very little coupling exists between the flowing liquid layer 
and the external gas boundary layer, especially if the liquid 
velocity is small compared with the external gas velocity and 
the melting temperature is much less than the gas stagnation 
temperature. The heat transfer, shear and pressure distribu- 
tion can then be regarded as known and equal to their values 
for a solid boundary at the melting temperature. On this 
basis several authors (13 to 18) have made theoretical studies 
of melting ablation. 


Theories for Melting Ablation 


The most work thus far has been devoted to the stagnation 
point of blunted bodies, although Lees (15), Lew and Fanucci 
(17) and Hidalgo (18) have treated ablation around body 
contours as well. Some of the references mentioned consider 
partial evaporation of the liquid layer; this problem will be 
discussed later. 

Sutton (13) first treated ablation at the stagnation point 
and integrated numerically the liquid layer equations for the 
case of Pyrex glass subjected to hypersonic flight conditions. 
He pointed out the important role played by the viscosity of 
the molten material, namely, the larger the viscosity, the 
larger the surface temperature and heat of ablation. Although 
Sutton’s approach is the most exact of the various theories, 
it is difficult to use for gaining an insight into the problem, as 
is usual with numerical analyses. 

Lees (15) obtained approximate analytical solutions for 
melting ablation under the assumption of constant viscosity 
through the liquid layer. He investigated the stagnation point 
and regions around a body contour to the sonic point, for both 
laminar and turbulent heating. Though not strictly applica- 
ble to glasses, which exhibit a large viscosity variation through 
the liquid layer, he applied the theory to Pyrex. Calculated 
heats of ablation were greater than 1125 Btu/|b at the stagna- 
tion point for re-entry velocities; the contribution from heat 
capacity up to the melting temperature was 750 Btu/Ib, and 
an additional 375 Bt/ulb was absorbed by the liquid. For 
turbulent heating at the sonic point, he calculated about 
twice the additional energy absorbed by the liquid or a total 
of 1500 Btu/Ib. 

Lees was concerned primarily with similarity parameters for 
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melting ablation. He showed that the stagnation enthalpy 
of the free stream gas is one of the principal parameters to be 
duplicated in model experiments and pointed out that rocket 
exhaust jets produce enthalpies much too low to simulate 
flight velocities in the range of 20,000 fps. 

An approximate analytical investigation of glassy materials 
in general was made by Bethe and Adams (16). This analysis 
accounts for the large viscosity variation through the liquid 
layer and treats stagnation ablation with and without evap- 
oration. A method is also outlined for calculations around a 
given body. By a second order analysis the theory is shown 
to be accurate to about 3 per cent in predicting the heat of 
ablation. In the absence of surface radiation, the result is ob- 
tained that the heat of ablation depends primarily upon en- 
thalpy or flight velocity, as found by Lees, but is inde- 
pendent of heat transfer rate or body size. Further, the 
theory predicts that the energy absorbed is quite insensitive 
to altitude or stagnation pressure. A simple formula was ob- 
tained for the stagnation point heat of ablation of Pyrex glass 
(for negligible radiation) in terms of free stream density and 
flight velocity 


‘is (he (Btu 
= 32 Pow ( ) (Fe 
Hett 320 ( ) RT, Ib [5] 


Interesting papers by Roberts (14, 19) consider stagnation 
point ice ablation with and without evaporation. For the 
case of melting only, Roberts obtains results similar to Lees’, 
indicating that some additional heat is absorbed by the 
liquid after melting. The total heat absorbed can be as much 
as twice the energy required to melt the ice. 

An important consideration for liquid layers is the stability. 
Though the liquid Reynolds numbers are, in general, low 
enough to insure laminar flow (15), it may be possible for some 
materials, under large deceleration conditions, to exhibit 
Taylor instability. Feldman (20) investigated this problem, 
and for the particular case of Pyrex glass he found the 
amplification notes to be quite small. 

Another effect of deceleration is the introduction of a body 
force in the direction of the liquid layer flow. References (8, 21 
and 22) show that the liquid flow will be retarded for 
some re-entry conditions, thus reducing the ablation rate. 
Under certain conditions the liquid flow direction may even 
reverse (21). 


Experimental Data 


Because of the lack of practical interest in purely melting 
ablation, only limited experimental data exist. The Chicago 
Midway Laboratories (3,23) obtained results on aluminum, 
copper and stainless steel, in both air and water arc jets. 
The excessive ablation rates make it difficult to interpret the 
data, but it appears that these materials absorb roughly 50 
per cent additional energy beyond melting, as might be ex- 
pected from the analyses of Roberts (14) and Lees (15). 


Conclusions 


Although studies of melting ablation are important from 
the point of view of understanding the overall ablation 
problem, materials which flow without evaporation are not 
of much practical interest. The heat absorbing capability of 
such materials represents only a minor improvement over 
the heat sink. This improvement may be considered im- 
portant if regarded as a safety factor for a material designed 
as a heat sink. 


Subliming Ablation 


_ Subliming ablation represents the other extreme from melt- 
ing, and the absence of liquid flow simplifies the problem in 
some respects. Sublimation differs from melting in that now 
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the interaction between the ablating material and the gas _ 
boundary layer is very important. The heat transfer is no 

longer independent of the ablation, as it was for melting; the < 
injected vapor leads to a reduction in heating by the well- 


known transpiration or heat blocking effect (24). = 


Consider an evaporating material in a steady-state ablating 


Simple Evaporation 


condition. The aerodynamic heat transfer to the solid 
boundary is absorbed by the heat capacity plus the heat of 

evaporation of the ablating material. In addition some heat a 
qr may be radiated away from the surface. In equation form 
the heat balance is given by 


mA(CpT; + hv) [6] 


The aerodynamic heat transfer q; is, in turn, related to the 
vapor injection rate m,. Here we draw upon the large body 
of work on transpiration cooling,‘ which gives a ratio of heat 
transfer with and without transpiration 


The quantity y is a function of mass injection rate, the 
molecular weight of the vapor, the external flow Mach number ~ 
and the nature of the boundary layer (i.e., whether laminar 
or turbulent). Utilizing available experimental and theoreti- 
cal transpiration results, several authors (7,24,25) have ob- — 
tained correlation formulas for the y-function. They found 

where (Ah)o defines the enthalpy difference across the bound- 
ary layer, without transpiration. The factor 8 will be dis- 


cussed later. 
Equations [6 and 7] then give 


CpTs + hy + B(Ah)o [8] 
1 — (€oTs4/qo) 

Thus the transpiration term can be interpreted as additional 
heat absorbed by the vapor as it diffuses through the gas 
boundary layer. From the references mentioned, 8 ~ 0.6 for 
laminar heating and roughly 0.2 in the turbulent case. Hence 
the transpiration term can account for a major contribution to 
Hs, especially for laminar heating. For satellite re-entry, 
for example, (Ah)o = 9000 Btu/Ib at peak heating. 

In the definition of Her as given by Equation [8], radiation 
emitted from the surface leads to an “apparent” additional 
energy absorbed by the material. Notice that in the limit of 
radiation equilibrium with the aerodynamic heating, Hes = 
co. In some cases it is more convenient to define a heat of 
ablation without radiation, i.e. 


(0) 
Het = C,T + he + B(Ah)o 19] 


This quantity is then independent of both heat transfer rate 
and radiation emission, and depends only on the enthalpy 
difference across the boundary layer. The former quantity 
(Eq. [8]) is usually determined from experiments, and conse- 
quently it will exhibit a heat transfer rate dependence, if the 
radiation is a significant fraction of the aerodynamic heating. 

The transpiration factor 8 is reasonably well-defined for a 
laminar boundary layer in air and is given approximately by 
the formula (7,24,25) 


= N(29/M,)* (10] 


where 0.67 < N < 0.72, 0.25 < a < 0.4, and M, is the 
molecular weight of the injected vapor. 
As usual, the turbulent case is not so clear, and only limited 


4 An excellent survey of transpiration cooling can be found in 
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theoretical and experimental results are available (26 to 29). _ 


air, then the data indicate that roughly 


Br ~ (1/3)Bz [11] 


ues of Br comparable with 8,. Ablating materials of concern — 
here, e.g., plastics and glasses, produce vapor molecular 
_ weights roughly comparable with air; hence Equation [11] 
seems to be the best estimate at present. 
___ The heat of evaporation h, should include the energy of dis- 
-- gociation of the injected vapor only if dissociated products 
appear at the surface (30,31). Lees (30) points out that the 
_ dissociation energy is beneficial only when the vapor dis- 
_ sociates right at the surface or when the dissociated vapor is 
able to diffuse back to a noncatalytic surface. 
' For some materials the surface temperature 7’; is not simply 
>, the evaporation temperature. For example, composite ma- 
terials often contain melamine or phenolic resins (2,5,8) 
_ which pyrolyze into a gas and leave a carbonaceous residue. 
This residue may then reach temperatures much greater 
than the pyrolysis temperature (8). There is no simple way 
to determine the surface temperature for these rather com- 
_ plicated materials, and experimental measurements are 
usually required. 


Effects of Combustion 


Several investigators (30 to 34) have studied the problem 
of combustion and other chemical reactions within a laminar 
_ boundary layer. A simple expression for combustion heating 
is derived by Hartnett and Eckert (34) under the assumption — 
that all of the oxygen diffusing toward the surface reacts with — 
material vapor. This condition leads to an upper bound 
estimate of the combustion effect, the increase in interface 
heat transfer rate being given by (34) 


0.21h, 
te.” 


- where h, equals the heat of combustion per unit mass of oxy- 
gen. Equation [8] is then modified (35) to read 
Hag = t + + 
1 — qr/qo1 + G) 
_ This same expression would be expected to apply as well for 
turbulent boundary layers, since the correspondence between 
conduction and diffusion, i.e., the Lewis number, should be 
4 similar to the laminar case. 


= [12] 


[13] 


Experimental Results 


Teflon 


Limited data exist for sublisaing materials, though results 
J are available for Teflon ablation at the stagnation point 
q (7,8). These data were obtained in air arc tunnels and are 
; : shown in Fig. 1, where Hes: is plotted vs. (Ah)o; the equiva- 

~ lent hypersonic flight velocities are also indicated. The 

Teflon ablation temperature was estimated (36,37) to be about 

800 F; hence radiation emission was negligible for the ex- 
ae perimental conditions. Two theoretical lines are also shown 
_ assuming either combustion or no combustion. The theoreti- 
eal lines were calculated as follows: 


1 No combustion 
(CaF"s) polymer vapor + hy 
h, = 750 Btulb (C.F,) (Refs. 36, 37) 
2 29 1/4 
B= 3 = 0.49 
CpT; = 0.25 X 800 F = 200 Btu/Ib 
From Equation [9] 
Here = 950 + 0.49(Ah)o, Btu/]b 


Combustion 
_ If the vapor molecular weight is comparable to or heavier than ere 


On the other hand, data for helium injection (27,29) give val- = an 


(CoF's) vapor + 2COF, + 
= 10,000 Btu/Ib(O.) (Refs. 36, 37) 
From Equation [13], with g, = 0 . 


Hest = + 0.49(Ah)o, Btu/Ib 


1 + [2100/(Ah)o] 


Either calculated curve gives reasonable agreement with the 
data. This is because the combustion effect, if present, is 
negligible in the range of enthalpies covered by the experi- 
ments. By fitting the theory to the data, the terms (C,7, + 
h,) are found to be 750 Btu/Ib instead of 950 Btu/Ib. This 
would indicate a reduced energy of evaporation meaning, per- 
haps, that the Teflon polymer is not completely broken down 
into monomers. 

The estimated Hs: for turbulent heating is also on Fig. 1. 
The turbulent estimate was obtained by multiplying the 
slope of the laminar data by one third, i.e., Br = (1/3) 8,. 
No turbulent data are available with which to compare. 


Phenolic-Nylon 


Considerable attention has been given a plastic material 
composed of a phenolic resin reinforced by nylon (see Steg, 
(8)). Although pure nylon will melt, the carbonized matrix 
residue of the phenolic traps the nylon melt causing it to 
vaporize within the matrix beneath the surface (8). Steg’s 
work indicates that the nylon phase change beneath the sur- 
face is very beneficial in reducing the effective thermal con- 
ductivity; hence the additional material required for insula- 
tion should be minimized by this effect. 

= theoretical formula derived (8) for phenolic-nylon gives 


Hen = 1500 + (1/2)he (14) 


which is 750 Btu/lb greater than for Teflon. Another ad- 
vantage of phenolic-nylon over Teflon is its capability of 
radiation emission. The observed carbonaceous residue 
should possess a high emissivity, in addition to yielding high 
surface temperatures (8). 

By use of an air arc wind tunnel, Steg measured Her: = 
4480 Btu/lb at the stagnation point for an enthalpy h, = 5100 
Btu/Ib. Equation [14] gives Hers = 4050, possibly indicat- 
ing some radiation benefit. 


Conclusions 


An important contribution to the effective heat of ablation 
of subliming materials arises from the well-known transpira- 
tion or heat blocking effect. This effect contributes an added 
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Fig. 1 Stagnation point heat of ablation measurements for 
Teflon, compared with theory 
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absorption energy over and above a material’s heat capacity 
plus energy of evaporation. The additional energy absorbed, 
per pound of ablated material, increases linearly with the en- 
thalpy difference across the gas boundary layer, and the 
proportionality constant is roughly three times larger for 
laminar heating than for turbulent heating. 


the Presentday air arc wind tunnels appear satisfactory for 
, is laminar ablation experiments; however, facilities to study 
eri- materials under turbulent heating are notably lacking. 

+ Surface radiation leads to a reduced net heat transfer and 
“his therefore an apparent increase in He. In this case Hes is 
Der heat transfer rate dependent. Data analysis is more difficult 
wn with radiation present, as now the surface temperature and 

emissivity are required. 

| 2 Mechanical erosion, though not discussed here, will reduce 
‘the Hs for some materials. Since surface erosion effects cannot 
Br be separated from ablation as such, the shear and pressure 


forces anticipated for application should be duplicated in 
laboratory ablation tests. 


Combined Melting and Sublimation 


teg, In this category, most attention has been given the glass- 
trix plastic combinations and high viscosity glasses alone. The 
t to glassy materials have been found attractive because of their 
eS high viscosity, low thermal conductivity and rather high heat 
sur- of evaporation. In addition, they exhibit good resistance to 


con- thermal stress. Plastics are combined with glass fibers or glass 
sula- cloth for several purposes, such as to simplify fabrication and 
, to reduce the thermal conductivity below that of the glass 
ves alone. 

Theoretical and experimental results are available and 


[14] comparisons will be made later, but first let us see how com- 

ail bined melting and sublimation influence Her. 

y of If a heat balance is made, as in the previous section, noting 

‘d that now only a fraction of the ablated material evaporates, 

idue 

high then it follows that 

5100 1 — (€oT4/qo) 

icat- The evaporation rate is denoted by m,, and the total ablation 
rate is m, as before. This formula differs from Equation [8] 
by the addition of another unknown quantity, the fraction of 
evaporation m,/m. Given a material with known physical 

—_ properties, the problem now is to determine the two unknown 

; quantities m,/m and 7. Theoretical analyses are available 

Oded which enable the calculation of these unknowns for ice (19) 
and glasses (9,16,38). Only the glasses, alone and in com- 
bination with plastics, will be discussed here as both experi- 
mental and theoretical results are available for these materi- 
als. 
Low Viscosity Glass in Composite With a Plastic 

<a An experimental and theoretical study of this composite, 


termed fiber glass, was carried out (35). It was observed that 
the glass (essentially the same as Pyrex) flowed without 
evaporation, as predicted by theory (16) and that the phenolic 
plastic component sublimed. It was further found that the 
surface temperature was determined by the flow of the glass 
(for the case of a composite in which glass constituted the 
—_ larger fraction of the mixture). A formula for Her can be 
easily derived for this case, if certain physical properties are 
assumed for the plastic. The molecular weight of the phenolic 
is comparable to that of Teflon; hence it seems reasonable to 
22 assume that their specific heats and vapor molecular weights 


ro. are also comparable. We assume further that their heats of 
evaporation are equal. The calculation of H.: will now be 
sed shown as an illustration of the procedure. 
It is convenient to calculate the heat of ablation without 
radiation, i.e., the numerator of Equation [15], and then 
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Fig. 2 Stagnation point heat of ablation results for quartz and p 
fiber surface radiation) 


the radiation correction factor can be applied for particular — 
heat transfer rates. In the example considered, the plastic — 
component is taken as 30 per cent of the total mass, the mix- | 
ture for which experimental data are available. 


Case a: Stagnation Point 


From Equation [5] 


Poo 1/13 hs 0.28 hs Btu 
C,T; = 320 = 300 { — 
(t=) Po RT, Ib 
My 


— [hy + B(Ah)o] = 0.3[750 + 0.49( Ah)o] < 


295 + 5 hi < hs 
RT. 


16 
RT, + 300 RT. [16] 


= 225 + 5 — 


Equation [16] is shown on Fig. 2,5 compared with data ob- 
tained in both an air arc wind tunnel (35) and a rocket ex- 
haust (2). This rather simple theory is seen to predict the 
data reasonably well. The close correlation between data — 
from these different facilities indicates that the chemical com- 
position of the test gas is not important for the fiber glass ma- 
terial. 


Case b: Estimates for Turbulent Heating 


Turbulent heating is considered here in a rather crude 
manner, only to give an indication of the differences to be 
expected. No data are available with which to compare the h 
calculations. 

As shown by Lees (15), the surface temperature for turbu- 
lent heating can be roughly 30 per cent greater than for 
laminar heating. Using this condition along with the as- © 
sumption that By = 8,/3 then 


1/13 0.2 
Po 


Equation [17] is shown on Fig. 3, but this result should be re- 
garded only as a crude estimate. 


5 The curve is drawn for l-atm stagnation pressure, and it is 
noted that the pressure dependence is very weak. 


629 


LJ 
“int 
| ‘@ 


ro 


¥ 


REF (18) 


QUARTZ 


(0) 


HEAT OF ABLATION, Hey ~(10° BTU/ #) 


37 
24 riper cLass 
J 
: FLIGHT VELOCITY ~ (10°F T/SEC) 
10 (2 14 16 18 20 22 


STAG, ENTHALPY, h, ~(10°BTU/#) 
Fig. 3 Estimates of the heats of ablation for quartz and fiber 
glass, subjected to turbulent heating — surface radia- 
tion) 
pe 


Results for Quartz 


Case a: Stagnation Point 


Of the many glasses, quartz (or fused silica) has the highest 
viscosity, and as a result it has been found to exhibit appre- 
_ ciable evaporation under laminar heating conditions (9). 

_ The theoretical calculations given in (9) can be approximately 
correlated to give the following result 

) —0.14 


h, Po 
= 1500 + 0.5 
1.54 14 
1500 + 0.55 ( (@) 

RT, Po 
The first term represents the heat capacity C,7,, which is 
roughly constant, and the second term is the contribution from 
the heat of evaporation plus the mass injection effect. Equa- 
tion [18], with ;4-atm pressure, is shown on Fig. 2. It is 
- compared with arc wind tunnel data of (9), and good agree- 
ment is observed. 

Rocket exhaust data exists (2) for a composite containing 
70 per cent quartz fibers in combination with 30 per cent 
phenolic resin. This material being predominantly quartz, its 
Hest can be reasonably compared with Equation [18]. The 
two data points given in (2) indicate Hf, ~ 5050 Btu/Ib as 
compared with a calculated value of 1790 for 10-atm pressure. 
The calculated temperature is approximately 4000 F com- 
pared with a value of 2790 F measured with an optical 
pyrometer. A possible explanation of these discrepancies 
could be that the optical pyrometer measurement is in error. 
If it is assumed that the true surface temperature is equal to 
the theoretical value of 4000 F, then the rocket data would be 


modified as indicated below, assuming a surface emissivity 


[18] 


qo 
Btu/ft?sec Btu/Ib —qr/qo Btu/Ib 
_ Pyrometer T = 2790 F 908 5390 0.06 5050 
Calculated T = 4000 F 500 3000 0.36 1900 


The close agreement between the modified (1900) and calcu- 
lated (1790) Hq values is probably fortuitous. 


Although this example may not be the proper explanation 
for the rocket exhaust data, it certainly serves to illustrate the 
difficulties encountered in data interpretation. The im- 
portance of knowing the surface temperature and radiative 
emissivity is apparent. 


Case b: Turbulent Heating 


The theory of (16) has been extended by Hidalgo (18) for 
quartz ablation around symmetrical bodies under both 
laminar and turbulent heating. His results show that with 
turbulent heating the fraction of evaporation is greater than 
for laminar heating; roughly this fraction is increased by a 
factor of two. An approximate theoretical Hgg curve can 
then be obtained by modifying the laminar results in Fig. 2. 
This modification involves doubling the fraction of vaporiza- 
tion and reducing the transpiration coefficient 6 by a factor 
of three as before. The results of this calculation are shown 
on Fig. 3; they should not be considered as accurate to better 
than 25 per cent. 


Effect of Radiation 


It should be remembered that the Hs; values shown on 
Figs. 2 and 3 were defined for the net heating, aerodynamic 
minus radiation. Quartz surface temperatures are, in general, 
large enough to emit significant radiation (under re-entry 
conditions the temperatures are calculated to be between 4500 
and 6000 F (9)). If the quartz is made opaque (e. g. by the 
addition of a carbonizing plastic (2)), then the magnitude of 
the radiant energy can be as large as 250 to 800 Btu/ft? sec. 
The heat of ablation, with radiation, can therefore be ap- 


proximated® by 
1 — 


> The importance of the radiation benefits thus depends upon 
the aerodynamic heat transfer rates. For satellite re-entry 
conditions, for example, the heating rates are relatively small, 
and H.s—> © provided the emissivity and nose radius are 
such that «VR > 0.35 (ft)! (see (9)). The surface will then 
be in radiation equilibrium, with negligible ablation. Radia- 
tion equilibrium has been experimentally observed at the 
Chicago Midway Laboratories (6). It was found that 
Foamsil (a form of opaque foam quartz) exhibited no ablation 
for aerodynamic heat transfer rates appropriate for satellite 
re-entry (100 Btu/ft? sec). 

Significant radiation emission can also be deduced from 
other tests (6) at low heat transfer rates. A quartz-phenolic 
composite gave heats of ablation in excess of 20,000 Btu/lb 
and glass-phenolic about 7500 Btu/Ib. This difference is ex- 
plained on the basis of a higher temperature for quartz, esti- 
mated at 4000 F , a8 compare 2d with 2400 F for the low viscosity 
glass. 


(0) 
Herr 


Her = [19] 


Conclusions 


It is evident from the results shown in Figs. 2 and 3 that the 
heat of ablation is strongly affected by the stagnation en- 
thalpy, as was the case with subliming materials. The heat of 
ablation is again heat transfer rate dependent only if the radiant 
energy emitted is a significant fraction of the aerodynamic 
heating. In general, the radiant energy will not be known a 
priori for a given material; hence to unravel experimental 
data it is necessary to measure or calculate the surface tem- 
perature and the emissivity; ideally, the net radiated energy 
should also be measured under test conditions. 

In contrast with Teflon, theory indicates that quartz will 
exhibit a larger ablation energy for turbulent heating than 


6 This approximation is subject to some error as H® is itself 
coupled with the radiation (9, 38). 
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for laminar heating. The larger fraction of evaporation, with 
turbulent flow, more than counteracts the reduced transpira- 
tion factor. Again, turbulent ablation experiments are 
needed before confidence can be placed in theoretical predic- 
tions. 


Examples of Applications 


Satellite Re-Entry 


As an illustration of the combined ablation plus insulation 
problem, let us consider the example of a recoverable satellite. 
We draw upon (8 and 35) in which heat shields for satellite 
re-entry are discussed. The satellite chosen for this example 
enters the atmosphere at 2 deg with the horizontal and is 
subject to a nearly triangular heat input vs. time. The heat- 
ing is assumed laminar with a peak heat transfer rate of 65 
Btu/ft? see and a heating time of 240 sec. The total inte- 
grated heating is then 7650 Btu/ft?. 

We will compare three materials: Teflon, fiber glass and 
opaque foam quartz. Because of the relatively low heat 
transfer rates in this example, only Teflon ablates appreciably. 
Foam quartz will radiate the aerodynamic heating, as pre- 
viously discussed. Fiber glass is also expected to radiate 
most of the aerodynamic heating (see (35)), although some 
mass loss will be experienced due to plastic decomposition be- 
neath the surface. In such a case it is convenient to define 
the insulation as a thickness required beneath the surface, 
and then the internal mass loss is charged as part of the insula- 
tion. 

To simplify the calculation of insulation weights, it was as- 
sumed that the conduction energy Q, had sufficient time, after 
the end of heating, to diffuse and heat the material to a uni- 
form temperature of 500 F. This assumption approximates 
the condition that the back face of the material, or the backup 
structure, will not experience a temperature rise in excess of 
500 F. The mass of material required to absorb the conduc- 
tion energy at 500 F is then termed the insulation weight. 

The results of this calculation are summarized in Table 1. 
See also (35 and 39). 

The long heating time associated with satellite re-entry 
leads to rather large insulation requirements. The insulation 
weight is proportional to V pK Tmax; hence foam quartz is 
better than fiber glass. The smaller insulation weight for 
Teflon arises from both a lower thermal conductivity and 
surface temperature. By comparison with the numbers in 
Table 1, a beryllium heat sink would require 22 lb/ft?. 

By contrast with satellite re-entry, a ballistic missile ex- 
periences much larger heat transfer rates and a greatly re- 
duced heating time. These differences change the magnitudes, 

as well as the relative importance, of the ablated and insula- 
tion weights. In the missile application, the ablated weight 
will in general be larger than the insulation. 


Rocket Nozzl 


In addition to —_ and other high speed flight applica- 
tions, ablating materials should also be useful for rocket 


oss 


nozzle walls. This application differs from re-entry in at least 
three important ways: 

1 The stagnation enthalpy is lower by about a factor of 
three. 

2 Surface radiation is not beneficial because of the in- 
ternal cylindrical geometry. 

3 The gas composition of the boundary layer is, of course, 
quite different. 

If it is assumed that the wall heating is not strongly in- 
fluenced by the gas composition, then the same heat transfer 
and ablation results, applicable for an air boundary, may be 
carried over to the nozzle problem. Some justification for this 
lies in the fact that ablation results obtained in a rocket ex- 
haust correlate very well with data in an air are wind tunnel 
(see Fig. 2). 

The boundary layers on the nozzle walls are undoubtedly 
turbulent; therefore in carrying over our ablation knowledge 
we must rely on theory. It is recalled that no turbulent abla- 
tion data are available at present. Theoretical estimates of 
Hem are given on Figs. 1 and 3; for rocket enthalpies (ap- 
proximately 2500 Btu/Ib) Teflon gives about 1000 Btu/lb, 
fiber glass 1300 Btu/lb and quartz about 2500 Btu/lb. These 
numbers are considerably smaller than those for re-entry con- 
ditions because of the reduced enthalpy. 

One important problem in relation to ablating materials is 
that of a changing nozzle area with time, the most serious 
location being at the throat. It can be shown that the frac- 
tional change in throat diameter is given by the approximate 
relation 

AD (Ah)o U 1 


D* pm (Re*)%? 


[20] 


where the starred(*) quantities refer to sonic conditions. 

Let us apply this formula to a particular case and determine 
the minimum throat diameter so that less than a 10 per cent 
increase will occur during 100 sec of operation. If we choose 
a specific impulse of 300 sec (h, = 2500 Btu/lb), a chamber 
pressure of 35 atm and a combustion temperature of 5500 F, 
the following result is obtained 


= 


for AD/D* <0.1. 

It is apparent from Equation [21] that a material having 
both a high ablating temperature and large Heq is required 
to obtain a reasonable minimum throat size. With quartz 
materials (Ah)o/H& ~ 1/2; hence the minimum throat 
diameter would be 3 ft (about 10° lb thrust). This size could, 
of course, be reduced if a larger change in throat diameter 
were tolerable. However, for throat dimensions less than 1 ft 
(less than 10° lb thrust), it appears that materials different 
from those discussed here are called for. For the same operat- 
ing conditions and requirements as above, a material with 

cx = 5000 Btu/Ib must have an ablation temperature 
greater than 5000 F to be satisfactory for a 1-ft throat diame- 
ter. 

Though the numbers quoted were obtained for specific 


Table 1 Heat shield weights for satellite recovery 


Avg. Hers Abl. wt Ins. wt. Total wt. 
Material Btu/lb- pK Q./Q lb/ft? lb/ft? lb/ft? 
fiber glass 0.125 ~3600 30.25 15.5 15.5 
opaque foam quartz © 0.068 3600 ies 0.14 0 8.5 8.5 
Teflon 4000 0.075 800 0.04 1.9 2.3 4.2 
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examples, it nevertheless would appear that new materials 
are required if ablation is to compete with internal cooling, 
g transpiration cooling or a heat sink (at least for the throat 
region of small rockets). Experimental work is certainly 
. needed to investigate materials for the rocket nozzle applica- 


Cp = specific heat 


= diameter 
h = enthalpy or energy per unit mass ptr) 
= enthalpy difference across boundary layer “4 


K = thermal conductivity 

k = thermal diffusivity, K/pC, 

H.+ = effective heat of ablation, including radiation emission 
Het = effective heat of ablation, excluding radiation emission 
= mass rate 

M = molecular weight 

q = heat transfer rate ie 
= time integrated heating 
= energy stored by conduction 
R. = Reynolds number, puD/p 
7 RT, = reference energy, 33.86 Btu/Ib 
= time 
temperature 
Tq = material ablation temperature 
Je = gas velocity 
7 Vy» = ablation velocity 
6; = thermal thickenss 
= transpiration factor 
= emissivity 
viscosity 
2 = ratio of heat transfer, with and without mass injection 
= density 
¢ = Stefan-Boltzmann constant 


Subscripts 
= combustion ; 

= free stream 

= solid-gas interface 
= laminar 

= material 

= condition for nonablating surface or sea level condition 


m 
0 

r = radiation 

T 
) 


= stagnation condition 
= turbulent 
= vapor 
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eliminated. 


NONLIFTING vehicle entering the Earth’s atmosphere 

at orbital velocity and grazing angles (@z < 1.5 deg) ex- 
periences & maximum longitudinal deceleration of about 
8.5 g (1,2,3).5 Peak G increases quite rapidly with increasing 
entry angle and reaches values of 11 g for #g = 3 deg and 16g 
for #g = 5deg. If entry velocity is near lunar or interplane- 
tary trajectory velocity (say 35,000 fps), not only does de- 
celeration build up rapidly with entry angle, but if the angle 
is too shallow there is danger of skipping out. These two 
limits allow an entry “corridor” less than 2 deg wide. Thus, 
there are obvious advantages in reducing peak G and pro- 
viding some degree of pilot control during entry. 

Eggers (1), Gazley (2) and Chapman (3) have shown that 
small lift-drag ratios permit large reductions in peak G during 
entry from decaying orbits (@¢ = 0 deg). But it is well- 
known that for such grazing trajectories the landing point 
is extremely sensitive to small errors in entry angle. The 
present paper adopts the point of view that the employment 
of aerodynamic lift eliminates the need for restricting the ve- 
hicle to grazing entry trajectories. In fact, we want to deter- 
mine the degree of flexibility in entry angle and peak G pro- 
vided by lift-drag ratios below the values required for an 
equilibrium glide path. At the same time we want to show 
that the “skip” phase can be entirely eliminated by means 
of relatively simple lift programming. 

Although there is a general agreement about the effects of 
lift on peak G, there seems to be a general impression that the 
use of lift always increases the duration of the heat pulse and 
the total aerodynamic heat transferred to the vehicle during 
entry. This statement is correct for a fixed entry angle, but 
is certainly invalid when we take advantage of the steeper 
entry made possible by lift. For example, according to 
Chapman’s calculations (Fig. 22b (3)) the total aerodynamic 
heat transferred and the maximum heating rate have about 
the same values for a nonlifting vehicle entering the Earth’s 
atmosphere from a decaying orbit, and a lifting vehicle with 
a constant L/D = 0.25 and 6g = 2 deg. Now the ratio of 
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Lift During Entry | 
Atmosphere 
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By employing aerodynamic lift during entry into the Earth’s atmosphere at either orbital disie 
cape velocity, the range of allowable entry angles for a prescribed peak deceleration is greatly in- 
creased, while the total heat energy transferred to the vehicle can be held to about the same value 
Only modest lift-drag ratios are required beyond peak G to prevent the _ 
deceleration from exceeding the peak value, or to prevent the vehicle from skipping out of the _ 
Thus, the difficult guidance and control problem is greatly alleviated; in par- 
ticular, for return from the moon or other planets, the necessity for multiple-pass drag braking is 


aerodynamic heating rate to the rate of dissipation of kinetic 
energy by drag is proportional to 
Vova 1 

so that this ratio is minimized by bringing the vehicle quickly 
to the lowest altitude consistent with the maintenance of 
laminar flow and limitations on peak G. The restriction to 
extremely small entry angles for the nonlifting vehicle favors 
a radiation-cooled design, while the employment of aerody- 
namic lift permits the designer to take further advantage of 
combined ablation and radiation, and to optimize the vehicle 
with respect to peak G, heat transfer and guidance require- 
ments. This general concept is illustrated in the present 
paper by means of a few typical numerical examples. 

Positive control over peak deceleration is even more im- 
portant for a manned vehicle returning to Earth from a lunar 
or interplanetary mission, i.e., at escape velocities. When 
drag braking is employed, the guidance and terminal control 
system must direct the vehicle to a corridor of the order of 10 
miles high. If the entry angle is less than about 4 deg the 
vehicle will skip out of the Earth’s atmosphere. On the other 
hand, for 6g > 4 deg, the peak deceleration experienced by a 
nonlifting constant-drag vehicle increases very rapidly with 
entry angle, reaching values of 9 g for 0g = 4.5 deg, 12 g for 
Oz = 5 deg and 20 g for 02 = 6 deg (Fig. 9, (4)). Eggers (1), 
Chapman (3) and others have studied the process of drag 
braking by means of successive passes through the Earth’s 
atmosphere, in order to limit peak G and reduce the velocity 
and eccentricity of the orbit to near-circular conditions. 
However, the guidance and control problem is hardly simpli- 
fied by this procedure. Drag modulation has been sug- 
gested (3,4) for direct entry, and Phillips and Cohen (4) have 
shown that the peak deceleration can be reduced to about 60 
per cent of its value for a vehicle with constant W/C DA, or, 
alternatively, the entry angle can be increased. However, 
the employment of aerodynamic lift should provide even 
more sensitive control over peak deceleration. Just as in 
the case of entry from orbital velocity, skipping can be 
eliminated by means of simple lift programming, provided 
92 > 4 deg, and the permissible margin of error in the guid- 
ance system is greatly increased. Again, these ideas are illus- 
trated in the present paper by means of some typical numer- 
ical examples. In a subsequent paper we hope to extend the 
analytical treatment of lift programming to the escape ve- 
locity case. 
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Equations of Motion 


A nonrotating two-dimensional inertial coordinate tee is 
employed with its origin at the center of the Earth. The 
gravitational field is assumed constant, independent of alti- 
tude. Because of the motion of the atmosphere, the magni- 
tude of the relative aerodynamic velocity V, does not equal 
the magnitude of the inertial velocity V. For the same 
reason V and V, are not co-linear in general, and the lift and 
drag will tend to curve the inertial flight path in the hori- 
zontal plane. This effect becomes very pronounced if the 
analysis is continued to low altitudes, where the velocity be- 
comes of the same order of magnitude as the rotational ve- 
locity of the atmosphere. 

The relation between V and J/, is illustrated in Fig. 1, while 
the coordinate system is shown in Fig. 2. From Fig. 1 it 
can be shown that 


where a is the course angle measured from east. 
In terms of the more commonly used angles, namely, those 
relative to the Earth’s surface 


V. 
= vl (1 — cos 6, cos — cos 0, a 


But V, = 1519 cos 6, where 6 is the latitude angle; V./V S 
0.14 over the portion of the trajectory where deceleration 
forces and heating are of most concern. Therefore in what 
follows, we replace V by V, in both magnitude and direction, 
so that the lift and drag force vectors lie in a vertical plane. 
The force components along the Cartesian inertial coordinates 
are 


F; = + Lain — mg sin = m(d/dt) (V cos 
[1] 

F, = Lcos@ + Dsin — mg cos = m(d/di(V sing) [2] 

By resolving these forces in the velocity direction, substituting 


for F; and F,, from Equations [1 and 2], and recognizing that 
6 = — y, we obtain 


gor —D + mg sin 6 = m(dV/dt) [3] 
Local Vertical 
‘ 


2 
) cos cos a | 


= 
In the same manner by recognizing that 


tude relation ¢ = a(h). 


dp _ dy _d0 
we obtain ; 
V cos 0 dé 


Assuming that V = V, i 
L = (1/2)pV°CLA 
= 
the equations of motion become > 
1dV _(1/2)pV?_ 
g dt W/CpA 
g dt W/CpdA \Cp g(R. + h) 


Consider first the case where entry velocity is the circular 
satellite velocity, or 26,000 fps. When lift is employed, the 
loss in vehicle velocity up to peak G is modest, as we shall 
show later. Therefore, in Equation [7] it is a good first ap- 
proximation to neglect the difference between gravitational 
foree and centrifugal force, compared with the lift force. Also, 
in the region of interest the contribution of the Earth’s 
gravity along the flight path is negligible compared to the 
aerodynamic drag, i.e., in Equation [6] - 


(1/2)pV* liane 
sin < A) 


With these approximations the reduced equations of motion 
are as follows 


V _ (1/2)pV? _ (1/2)p0 
ian v2 8 
(W/CpA) A ” 
where 4 = W/C)DA, andé 
—V6 _ (1/2)eV? Cx _ (1/2) (2) 
g (W/CpdA) Co Co 


Simple Analysis: Constant Lift-Drag Ratio to 


Peak G 
By dividing Equation [8] by Equation [9], we obtain 


[10] 


[11] 


independently of any statement regarding the density-alti- 
For constant Cz/Cp 


Co 
(Oz — 


In order to relate ¢ and @ some definite statement must be 


Ve 
— = ex 
V p 


we 


[12] 


® These equations are similar to Equation [14] employed in 
(5) to analyze the skip trajectory. However, the ap roximations 
employed in arriving at these equations are quite di erent in the 
two cases. (Note eo that @ in (5) denotes the local inclination 
of the velocity vector to a fixed horizontal direction, whereas 
here @ denotes the angle between ¥ and the local horizontal.) 
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made about o(h). For the present we take 


o = p/p. = 


where po is adjusted for best fit (1,5) 


Po 1.5Psea ievel 
Since 
h = —Vsin0 [13] 
go = —Boh = BoV sin 0 [14] 


By dividing Equation [9] by Equation [14], one obtains 


dé 1 / pg Cr 
sin = = -3(%)() 


For constant C,A/W 


—or=2 [16] 


which for small angles reduces to 7 7 


GE) 
Pog 


Now the magnitude of the deceleration vector as felt by a 
pilot or an accelerometer is given by 


ll 
— 

— 


or 


By utilizing Equations [12 and 17] the following expression is 
obtained for G along a trajectory with constant C,/Cp in 
terms of entry conditions - 


C. 
exp [-2 (08 ~ [18] 


For constant (C,/Cp) and constant (W/CpA) = A, peak G 
(denoted by G*) is experienced when oV? is a maximum, or 


when d/dt(oV?) = 0. But le 


d 

(oV?2) = + 2VV = o (sis 0 
according to Equations [8 and 14]. Therefore G = G* when 


= (poy/BA)o* 6* [19] 


sin 6* 


for small angles. 
By combining this relation with Equation [17] we obtain 


[20] 


which can be combined with Equation [18] to obtain G* ig. 
3). 

By neglecting og in comparison with o*, we obtain 


on. 


Lift "ce. g.'"' of Moving Vehicle 


Local Horizontal 


= 
Velocit 
= 
» de 


Fig. 2 Inertial coordinate system 
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Fig. 3 Illustration of the relationship between 6 and o 


and 


[22] 


Gt = 


where 6* is given by Equation [21]. 
When, in addition, 


<1 
then 
6* = - [23] 


(2) | - (2 | 
G 2 1+ Oz? exp 2 Oe 


In Fig. 4 the value for peak G with lift is compared with the 
maximum longitudinal deceleration experienced in a ballistic 
(nonlifting) trajectory for entry from a circular satellite orbit 
(Vz = 26,000 fps). When Cz/Cpy < 1, the approximation 
contained in the reduced Equation [9] of the present analysis 
is not applicable, and the peak deceleration for Cz/Cp = 0 
must be obtained by numerical integration. According to 
Fig. 4, peak deceleration can be reduced to values as low as 
2 to 3g over a practical range of entry angles by employing 
modest amounts of lift. Appreciable reduction in peak G is 
obtained with a value of C,/C as low as 0.25, which can be 
produced by such simple means as an offcenter shift in center 
of gravity location on a blunt nosed vehicle. 


V*? = Vg" exp [—2(Cp/Cz) 6s] 


and 
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Lift Cutoff 


In the preceding section we considered the simplest case of 
lifting entry trajectory from a satellite orbit, in which C,/Cp 
is held constant up to a predetermined peak value of G = G*. 
The question naturally arises as to the proper time at which 
the lift can safely be reduced to zero, with the assurance that 
Gz* < G* over the ballistic portion of the trajectory. In 
most cases 9* is quite small, and analytical solutions for the 
initial portion of a ballistic trajectory beyond G = G* are 
difficult to obtain. However, from the equations of motion 
it is not difficult to verify that: 

1 Ina ballistic trajectory, G decreases monotonically be- 
yond a local maximum or stationary value of G when 6* is 
sufficiently large. 

2. G initially increases when the lift is “turned off” at peak 
G if 0* is small. 

3 Only modest values of lift are required to prevent G from 
initially increasing beyond peak G. soit 

These conclusions will be demonstrated. 

In the last section we saw that 
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“Fig. 4 Peak G with constant (Vz = 26,000 fps) 
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By considering the equations of motion ([7 and 14]) we find 
that 


—(1/2)pV? CL ve 


do BoV? sin 0/g 
(oV?) = 0 


(-2)-(- 


For a ballistic trajectory, the inequality 


gR.} \ g 


is the necessary and sufficient condition for a monotonically 
decreasing G beyond a local maximum or level value in G. 
However, for interesting values of 6g and C,/Cp, the values 
of @* and (— V/g)* are almost always small enough so that 
(d?/dt?) (oV?) > 0 at peak G, and some lift is required. Ac- 
cording to Equation [26] when C,/Cp ¥ 0, (d?/dt?)(a V2) < Oat 
peak G as long as 


CL 


(V2/gR.) 


— 20 
—(V/g) 


For example, if V* = 20,000 fps, (— V/g)* = 3, and 6* = 
0.01 radian, then initially (d?/dt?) < Oif > 0.11. 
The fact that (d?/dt?) (oV?) < 0 initially does not guarantee 
that G will continue to decrease; therefore, we must consider 
the remainder of the trajectory. According to Equation {7] 
aerodynamic lift must provide an opposing normal acceleration 
of the order of g(1 — V?/gR.), but the longitudinal decelera- 
tion is of the order of several g, so that only modest values of 
C1/Cp should be required. Subsequently, lift can be turned 
off when sufficient velocity or kinetic energy has been dis- 
sipated by drag, and the values of @ and o are large enough 
so that the inequality given by Equation [27] is satisfied. 
Small values of C,/Cp are sufficient to keep G below (— V//g)* 
throughout the rest of the trajectory, as one can show by con- 
the following simple trajectory’: 1 constant C,/- 
_= (C1/Cp): up to peak G; 2 constant G = G*/- 

1+ (C1/Cp):? = (- V/g)* beyond peak G up to the point 
at which the required lift is zero. In this case, where ter- 
mina] lift i we small - G = constant, 


From Equation [26] these conditions yield 


Cu _ 1 — (V?/gR.) 
(—V/q) 


- (22)(5) 


7 This simple lift program example i is by no means optimum; 


— 206 


[28] 
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Two special points are of interest here: 
1 The condition for maximum (C,/C») and the 
(Ci Cob)m. 
2 The conditions for Cx,/Cp = 0. By differentiating 
Equation [28] with respect to o, one finds that C1/Cp is a 
maximum when 


A/m pg YR. 
and 

BR. 


independently of the value of (— V/g)*. 
value of 


1 4 
|  (—Vig)* VBR, 


For the planet Earth, 


value of 


maximum 


The 


[29] 
independently of the value of W/C)A. 
the numerical values are 
(o/A)m = 1.26 10-5 
6, = 0.033 rad = 


}:vidently the values of C,/C» required to maintain G con- 
stant beyond peak G are also small; for example, when 
(— V/g)* = 3, = 0.20. Also, when 

R. 


q 


1.9 deg 


— 0.132 


(—V/g)* 


(—V/g)* > 7.6 (Earth) 


the required lift is always negative.’ 


According to Equation [28], C1/Cp = 0 when 6 = 6 where 


BR. = — N?2) [30] 

and 
VBR. 


Thus the lift passes through zero when @ = 6,, the larger of - 
the two roots given by Equation [80]. From Fig. 6 one can 
see that 04 is generally large enough so that lift can be set 
equal to zero here without danger of (— V/g) peaking sub- 
sequently. 

Returning to the consideration of 6,, two cases can be dis- 
tinguished 


BR. 
« 
Viz 


In the first case, the curve of C,/Cp vs. o has a distinct maxi- 

mum beyond @ = 6* (Fig. 5), where the value of (Cx/Cp)m 

is given by Equation [29]. In the second case, Cz/Cp 

decreases monotonically with increasing o beyond peak G 

(Fig. 5), and Cz/Cp is always negative when 6* > 0+. A 

curve of 0:4 vs. (— V/g) is presented in Fig. 6. 
Now, 6* and 


case il 


( G* 

a V1 + 

: Of ¢ course, in the practical application the concept of nega- 
tive lift being required in this phase is meaningless, since zero 


lift will allow G to decrease. 
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Fig. 6 Angle at which lift becomes zero with constant G = 
(—V/g)*; N = [4(—V/g)]/VBRe 
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C,/c, Prior to Peak G's 


6 3 
Entry Angle, 8, (deg) at 400, 000 ft 


Fig. 7. Entry conditions for which no lift is required (shaded 
region) to hold G = (—V/g)* beyond peak G (Vz = 26,000 fps) 


are determined as functions of (Cz/Cp); and 0g for Vz =— 
26,000 fps in Equations [21 and 22], from which the regions 
in the (Cz/Cp); — plane associated with two cases 
discussed are readily ascertained. In Fig. 7, the area below — 
the curve 


- \ BR. 
represents all pairs of values of Cz/Cp and 6 for which 6* is | 
greater than V/1 BR., case ii. Similarly, the area below the — 

curve (— V/g)* = 7.6 in Fig. 7 represents conditions for which 

(— V/g)* > 7.6; here C_/Cp required to maintain constant 
G beyond peak G is always negative (Eq. [29]). One sees | 
from Fig. 7, when 6* > 1 V/ it is also true that (—V//g)*> 
7.6, except for a small region of the diagram corre-— 
sponding to low entry angles and small lift-drag ratio. We _ 
conclude that small lift-drag ratios fo rg to prevent 
G from exceeding (— V/g)* beyond peak G for (— V/g)* <'7:6. 

Our simple example used the case of G = (-V/g)* = 
constant beyond peak G. It is of much more in- 
terest to ask how much lift is required to prevent G from _ 
exceeding G*. Fig. 8 shows the amount of lift necessary to | 
prevent G after lift cutoff from exceeding the original G*. q 
These results were obtained from computer calculations. — 
For these calculations C,/C after peak G (terminal lift) was — 
held arbitrarily constant. 


is significant in peak G the 
values calculated for constant C,/Cp can be obtained by 
proper programming of the lift. 


= 
BRL 
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[28] Lift Programming 
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Fig. 8 Cz/Cp required after peak G to prevent the subsequent 
deceleration from exceeding peak G(Vz = 26,000 fps) 
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Fig. 9 Sketch of lift and G variation for step reduction in lift- 


drag ratio 


! 
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Simple Example: Step Reduction in Lift-Drag Ratio 


As a simple example of the gains that can be realized, con- 
sider first the case of a stepwise reduction in C,/Cp at o 
o;, and suppose one requires that peak G = G,* should be no 
larger than G,, the value of G at o = o (Fig. 9). We are 
particularly interested in the ratio G.*/G,*, where G,* is the 
value of peak & that is experienced by the vehicle if C,/C, = 
(C',/Cp), throughout phase 1. 

By utilizing relations like Equations [12 and 17] for the 1 
portions of the trajectory with constant Cz/Cp, one obt: 


where 

Cr 1 + Cr 2 CL 
Furth 
Furthermore, when 


then 


In this case, by utilizing Equations [18, 22 and 23] and re- 
quiring that G.* = G, (Fig. 9), one finds that 


— (612/08) [ (°°) 
xp | 2 6%) 
(Cr/Ci) 
V1 + 
* 


[34] 


6,2 
(: exp [2 (22) (0, 0; | 


Two numerical examples will give some idea of the quan- 
titative effects of stepwise lift reduction. 


G* 


Case 1 


= 0.10 radian = 5.73 deg; (Cx/Cp); = 1.0; = 
0.50. From Equation [33] one finds that 0,/@2 = 0.45, or 
6, = 0.045 radian = 2.6 deg, and from [34] G.*/G,* = 0.88. 
Since G,* = 5.4, G.* = 4.7. 


Case 2 


Oz = 0.15 radian = 8.6 deg; (Cz/Cp), = 2.0; (Cx/Co)2 = 
0.50. From Equation [33], 0,/0¢ = 0.555, or 0, = 4.75 deg, 
-and from Equation [34], G.*/G;* = 0.75. Thus the reduc- 
tion in peak G is small when (C,/Cp); = 1.0, as expected, 
but becomes appreciable (= 25 per cent) when (Cz/Cp); = 
2.0. Somewhat larger gains can be achieved in this last 
case by continuous variation of C,,/Cp up to peak G. 


Continuous Lift Programming 


We have seen how peak G can be reduced by a step re- 

duction in Cz/Cp. It follows that two steps will reduce the 

peak even more, and an infinite number of steps or continuous 

programming will produce the greatest reduction in peak @ 

(Fig. 10). 

7 There are two conditions to be satisfied: During the lift 

~ modulation period G will be constant; at the end of the period 
the lift will be zero. 

Despite the fact that the deceleration is at its largest value, 

the quantity which shows the least change during this perio: 
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If one assumes the velocity to be constant 
C 
Vet exp -2 (Or — ] 135] 


(1/2 C 2 
Co/i 


A 


4 4 


Cp o 


Substituting this last expression into Equation [7], one ob- 


tains 
,  (1/2)pog __, KE g ( v2 
— a} -) -1-=(1-—) [87] 
A o V gR. 


= BoV sin 0 = 


Carrying out the integration and employing Equation [36], 
we find that 


is the velocity. 


then from 


one obtains 


14 


Peak Deceleration, 


Fig. 11 Peak G using modulated lift (Vz = 26,000 fps) 


where V? can be approximated by 


» 
Cp 
6 J exp | -2 (2) 
But at o’ = 
dG 
a 
From Equations [19 and 36] 


Also 


C (oz — + Oz? [40] 


By combining Equations [38, 39 and 40] we obtain a quad- 
ratic equation for o;, the solution of which is 


tps) for W/CpA = 100. In Fig. 12 Cz/Cp is constant up 

to G*, and corresponds to Fig. 4, where Vz =26,000 fps. — 
Fig. 13 shows G* for the modulated case, and corresponds to 
Fig. 11. Fig. 14 shows the effectiveness of lift modulation | 
compared with constant lift. This figure applies for both — 
Ve = 26,000 fps and Vz = 35,000 fps, and is independent of © 


6x. Finally, Fig. 15 shows the effect of lift on heat rate and — 


total heating for Vz = 35,000 fps, holding G* = 10 g (modu- 
lated lift). 

The question naturally arises as to the lift modulation re- 
quired beyond peak G in order to prevent the vehicle from | 
skipping out of the Earth’s atmosphere. For steep angles — 
and low C,/Cp, no lift is required beyond peak @ (Fig. 16). 
For moderate angles, some negative lift is required beyond peak 
G, but naturally higher negative C,/Cp is required for the case — 
in which a constant positive Cz/C'p is used prior to peak G | 
than for the case of positive lift modulation prior to peak G. | 
However, in no case is the negative Cz/Cp required greater | 


* 


AB \Cd/; 


) on + + x] 


[41] 


| 2 


Thus utilizing Equations [35, 40 and 41] together with the 
entry conditions, one can determine G* obtainable with con- 
tinuous programming of C;/Cp. The G* using lift modula- 
tion is shown in Fig. 11; for a comparison with the values of 
G* for a constant lift see Fig. 4. 


Computer Results for Entry at 35,000 fps 


In this section, electronic computer calculations (IBM 704) 
were made for the case of lunar flight return (Vez = 35,000 
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than the positive C,/Cp used prior to peak G.® At first 
glance this result may seem surprising, but if one looks at the — 
expression for (— V/g) (d@/dt) one sees that the term (—1 + | 
V2/gR.) cannot be greater than 1, so the term (— V/g) x _ 
(C1/Cp) need only be equal to 1. Hence C,/Cp need never | 
be larger negatively than —1/(— V/g) to prevent skip. At 
this point (— V/g) is certainly greater than 4-5 (unlike re- 
entry conditions), and hence only modest values of negative 
C1/Cp are required. See Fig. 16. 


® For entry angles less than 5.5 deg, negative lift must be ap- 


plied immediately on entry to prevent skipping. However, for 
02 < 4.5 deg, the negative lift required is unreasonably large, as 
expected. 
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Fig. 12 Peak G (G*) with constant C,/Cpy (Ve = 35,000 fps) 
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Fig. 13 Peak G using modulated lift (Vz = 35,000 fps) 
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Fig. 14 The effectiveness of lift modulation (independent of 
62 or Ve) 


Conclusions 


By employing aerodynamic lift during entry into the 
Earth’s atmosphere the range of permissible entry angles for 
a prescribed peak deceleration is greatly increased, yet the 
total heat energy transferred to the vehicle can be held to 
about the same value as for a nonlifting vehicle. For entry 
from orbital velocity, analytical expressions are obtained {or 
peak G in the simple, illustrative case of constant lift-drag 
ratio up to peak G. For example, this analysis shows that the 
entry angle at 400,000-ft altitude can be as great as 9.5 deg 
with a lift-drag ratio of 2.0 and a peak deceleration of 10 g as 
compared with a maximum allowable entry angle of less than 
3.0 deg for a nonlifting vehicle (see Fig. 4).1° The analysis 
also shows that beyond peak G, only modest lift-drag ratios «re 
required to prevent the deceleration from exceeding the peak 
value, and the unnecessary skip phase is entirely eliminated 
(Figs. 7 and 8). 

By employing continuous lift programming or lift modul!a- 
tion, the maximum permissible entry angle is increased for a 
prescribed peak G, or, alternatively, peak G is reduced. For 
example, if peak G is held to 10g, an orbital vehicle can enter 
at an initial angle of 12.5 deg, if the lift-drag ratio is initially 
equal to 2.0. Alternatively, for an entry angle of 9.5 deg, peak 
deceleration is reduced to 6 g (Fig. 11). 

The effect of aerodynamic lift on the trajectory is even more 
pronounced for a vehicle returning to Earth from outer spice 
at the approach velocity of about 35,000 fps. By employing 
positive lift up to peak G and negative lift beyond peak G, 
the vehicle is prevented from skipping out of the Earth’s at- 
mosphere for an entry angle as low as 5.5 deg. By using 
negative lift initially, the minimum entry angle is reduced to 
about 4.5 deg (Fig. 16). When the maximum deceleration is 
held to about 10 g, the maximum allowable entry angle is 
10.5 deg for a lift-drag ratio of 2.0 (Fig. 13). Thus the diffi- 
cult guidance and control problem is greatly alleviated, and 
the necessity for multiple-pass drag braking is eliminated. 
Total heat energy transferred to such a vehicle with a 
W/CpA = 100 is about (40,000 Btu/ft?)/~/ Ro, and the maxi- 
mum stagnation point heat transfer rate is about (800 Btu/ft?/ 
sec) in/ Ro, where Rp is the nose radius in feet (Fig. 15). 


Nomenclature 

Cp = drag coefficient 

Cri = lift coefficient 

D = drag 

G = deceleration experienced by an accelerometer, g 
= altitude above Earth’s surface 


_ 1° The reader is reminded that the entry angle definition must 
include the entry altitude taken as 400,000 ft in this paper. 
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Fig. 15 Re-entry heating for 1-ft radius hemisphere with surface 
temperature = 2100 F. The trajectory is flown so as to keep 
G* = 10g (Ve = 35,000 fps, W/CnA = 100) 
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Fig. 16 Entry conditions which do not skip for: Negative lift at 
entry, no lift after G*, and programmed lift (negative) after G* 


eee 


lift 

radius of the Earth 
inertial velocity 
velocity relative to the atmosphere — 
local velocity of the Earth’s atmosphere in inertial space 
vehicle weight 


angle between east and the horizontal projection of V 
angle between east and the horizontal projection of V, 7 4 


g/RT ~ (1/23,000) ft-! = (1/h) log (po/p) 
latitude 
(W/CpA) 


ordinate of inertial coordinates 


6 = angle between local horizontal and V 
9, = angle between local horizontal and V; 
— = abscissa of inertial coordinates 
p = density ~ 
po = 1.5 X sea level reference density ~ 0.0034 slug/ft® 
= p/p = eh 
¢ =6+y 
= angle between local horizontal and 
Subscripts 
B = during the ballistic or no-lift phase S=_ 
E = entry condition 
1 = at the beginning of the lift modulation period 
2 = during the C1/Cp step phase 
m = maximum ot 


‘ = end of lift modulation period and beginning of ballistic phase 
* = condition of maximum deceleration . 
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URRENT interest in satellites, space travel and re-entry 

has resulted in numerous analyses of aerodynamics at 
hypervelocities, that is, at vehicle speeds of up to the escape 
velocity. Experimental studies at these extreme conditions 
are required to augment the theory. Research at these very 
high velocities is often performed in a shock tube, but higher 
velocities and temperatures than can be obtained in pressure- 
driven shock tubes are necessary. One new method of ob- 
taining shock waves at higher velocities is by use of the elec- 
tromagnetically driven shock tube. 

Various types of electromagnetic shock tubes have been 
used previously for generating strong shock waves to obtain 
very high gas temperatures, but their use to produce high 
velocity gas flow for aerodynamic studies is new. Before 
any one configuration can be considered for aerodynamic re- 
search, an investigation must be made of the suitability of 
the flow conditions existing in the shock tube. The first 
part of this paper is a brief report of the successful perform- 
ance of a tapered tube-type electromagnetic shock tube. 

The remainder of the paper describes the use of this shock 
tube for an experimental investigation in magneto-aerody- 
namics, which results from the interest in the application of 
magneto-aerodynamics to very high speed vehicles. When 
a vehicle travels through the atmosphere at speeds above 
approximately Mach 15, the air behind the shock wave in 
the region near the nose becomes thermally ionized and, hence, 
electrically conducting. Theoretical studies have shown 
that the force produced by the interaction of the electrically 
conducting fluid with an applied magnetic field can produce 
appreciable changes in the flow pattern about the body. 

Although considerable theoretical work has been done on 
this subject, very little experimentation has been reported. 
_ Therefore, an experimental investigation was made of the 
effects of an applied magnetic field on the bow shock standoff 
_ distance of a blunt body. A report of this study and a com- 
parison of the results with theoretical predictions are given. 
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Experimental Investigation in 
Magneto- Aerodynamics 


ao A 3-inch diameter electromagnetic shock tube is studied as a means of producing the hyper- | 
velocity and ionized air flow required for magneto-aerodynamic experimental research. The simple 
construction and operation of the shock tube is described, and an evaluation is made of its perform- 
ance. It is shown that although usable test times are only about 20 microseconds, velocities up to 
12,000 meters per second, and stagnation temperatures up to 25,000 K are readily produced, well into 
the range suitable for magneto-aerodynamic studies. The electromagnetic shock tube is used to 
study magneto-aerodynamic interaction in the ionized flow about a blunt body. The test body isa 
hemispherical cylinder containing a coaxial pulsed magnet coil in the nose, which produces field 
strengths at the stagnation point of up to 40 kilogauss. Quantitative measurements are made of 
the change in the bow shock standoff distance upon application of the magnetic field. The mag- 
netic field displaces the bow shock upstream and the standoff distance increases by a factor of 
7.5 for the condition (0; ;*rp/p,U2) = 69. The experimental results are com pared with theory, and 
good agreement is obtained within the range of experimental conditions. ; a 
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Electromagnetic Shock Tube 


In the tapered tube-type electromagnetic shock tube (1)* 
which has been studied, the driving energy is suddenly added 
to a portion of gas at one end of the tube by a high voltage 
electrical discharge. In addition to the rapid expansion of 
the heated gas, a magnetic pinch is created simultaneously 
with the discharge which ‘‘squirts” the expanding hot gas 
down the shock tube. Together, these effects generate a very 
strong shock wave with a Mach number as high as 100 or 
more, with accompanying high velocity flow and high temper- 
ature conditions behind the shock front. 


The shock tube and associated equipment consist mainly 
of the shock tube assembly, spark gap switch, high voltage 
power supply, condenser bank for energy storage, trigger 
pulse circuit and vacuum system. These are shown sche- 
matically in Fig. 1. The shock tube proper, shown in Fig. 2, 
is a commercial 3-in. diameter double tough Pyrex glass pipe, 
and the tapered tube section is formed from Pyrex glass. 
The shock tube is evacuated to the operating range of 25 to 
300 u absolute pressure, and the pressure is controlled by a 
needle valve bleed which admits air or bottled gases to the 
shock tube. The bleed also provides flushing of the tube 
between shots to remove contaminants. 

The electrical energy is supplied from six 1-uf capacitors 
which are charged to 25 kv, and each is coupled to the dis- 
charge assembly through RG-8/U coaxial cable. The elec- 
trodes consist of a central electrode at the bottom of the 
tapered section and a ring electrode between the cone and 
tube sections. The energy in the capacitors is transferred 
through a spark gap switch consisting of two slightly convex 
copper plates, separated by a ceramic insulator. The upper 
plate is made integral with the central electrode. The spark 
gap is pressurized with nitrogen to about 10 cm Hg gage pres- 
sure to reduce corrosion of the plates, and the entire switch 
is submerged in transformer oil to prevent high voltage leak- 
age. 

The spark gap prevents current flow until a 15-kv pulse 
is applied to the trigger wire in the. lower switch plate. A 
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discharge occurs between the wire and the lower plate, ionizing 
the gas in the gap and making the gas conducting. Conse- 
quently, the capacitors discharge across the gap and then 
across the electrodes, with a ringing frequency of 125,000 eps. 
The current flows from the center electrode to the ring elec- 
trode and then returns through six copper straps placed 
evenly around the conical tube. In addition to the sudden 
energy release by the joule heating, the magnetic field created 
by the current in the ionized gas pinches the heated gas. 
The heating and the self-pinching together, generate a strong 
shock front which emerges from the discharge section. The 
shock wave then propagates along the tube to the test region 
near the end of the tube, accelerating the gas to very high 
velocity. 
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Fig. 1. Schematic diagram of the electromagnetic shock tube 
and allied equipment 


Fig. 2. Electromagnetic shock tube in operation. Test body is 
installed on shock tube end plate 
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_ which provides extremely short exposure times and accurate 
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Photographic Techniques of Analysis ; 


In the study of shock tube performance, useful and reliable _ 
data have been obtained with direct, streak and image con-_ 
verter photographic techniques, each one serving a specific: 
role. Direct photography has been used to obtain good 
qualitative pictures of stationary shock wave formations. 
The camera shutter is opened just prior to the firing of the 
shock tube, and the exposure made by the luminosity of the _ 
gas in the shock tube. Although the fluid is generally lumi- — 
nous after passage of the incident shock, the luminosity behind — 
the shock waves attached to a model is much greater and eas-_ 
ily distinguishable. 

For the observation and measurement of dynamic phe- 
nomena, such as the travel of the incident shock wave along 

_ the tube, a rotating mirror-type streak camera (2) has been 
used. The streak picture is essentially a time vs. distance _ 
- of the shock wave travel, and the velocity of the shock 
front can be measured from the slope of the image. In par- 
_ ticular, the streak camera was used to measure the shock front — 
velocity and attenuation along the tube. Typical streak 
pictures are shown in Fig. 3. The vertical streaks are light 
reflection off the model nose and, later, illuminosity of the 
} bow shock wave. 
The image converter camera (3) is an electronic camera — 


_ synchronization with external events. The camera is actu- 
cite by a voltage pulse, and takes one exposure or a series of 
_ three exposures with an exposure time of between 5 to 200 
-millimicrosec on either Polaroid or conventional cut film. A 
pulse from a pickup coil, which occurs simultaneously with 
the gap breakdown, is used as the synchronizing event with a 
variable delay connected in series. Hence, the exact time 
of the camera exposure can be preset to within a microsecond. 
This camera has been used to photograph the traveling 


Fig. 3 Streak photographs of traveling shock waves in air at 
M, = 38, 18 and 15, left to right, respectively. Flow is from left 
to right, time increasing downward 
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Fig. 4 Measured effect of shock tube initial pressure on shock 
Mach number in air 
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wave front as it arrives at the model and the resulting transi- 
tory phenomena during the establishment of equilibrium shock 
conditions. It also has been used to measure, quantitatively, 
bow shock standoff distance during the steady flow test 
period. : 


Performance Evaluation 


An investigation was made of the performance character- 
istics of the taper tube-type shock tube together with an 
evaluation of its potentialities and the limitations. In par- 
ticular, the shock wave velocity as a function of initial pres- 
sure was determined, and the uniformity and duration of flow 
in the test section were analyzed. Also, the thermodynamic 
properties of the gas behind the shock front were determined 
for the range of conditions encountered in the shock tube. 

The shock velocity was measured as a function of the ini- 
tial pressure in the shock tube for a constant charging voltage 
of 25 kv. From the data shown in Fig. 4, it can be seen that 
the shock Mach number is inversely variant with the initial 
pressure. Shock Mach numbers in the range of 15 to 40 can 
be readily obtained with air. From streak pictures, such as 
in Fig. 3, it was found that from the midpoint to the end of 
the tube there is little attenuation of the primary wave 
velocity. However, secondary or successive shocks were 
evident in the streak pictures. These weak shocks, probably 
resulting from discharge oscillations, will cause small irregular 
variations of the density and temperature in the free stream. 
Experiment does show, however, that in spite of the presence 
of these fluctuations, the flow Mach number is reasonably 
uniform, as evidenced by the uniformity of stationary shock 
waves on the test models. 

In order to determine the uniformity and duration of the 
high velocity flow, a study was made of the formation and 
constancy of the bow shock on a 2-em diam hemispherical 
body. Image converter pictures were taken which showed 
the shock layer thickness and symmetry at various times 
after the incident shock arrival. From several series of photo- 
graphs at different initial pressures, it was concluded that 
the establishment of the bow shock wave requires about 3 or 
4 microsec, and that after this, steady conditions are main- 
tained for 10 to 35 microsec. Longer durations of steady 
flow were obtained as the initial pressure was lowered, and 
this variation is shown in Fig. 5. 

Next, the state of the shocked gas is to be considered. 
The thermodynamic properties of air, including density, pres- 
sure and temperature, behind a traveling shock front have 
been computed by Gilmore (4) for shock Mach numbers of up 
to about 50. The data, based on real gas properties, assume 
that the shock front is a single, normal and plane discon- 
tinuity and that thermal equilibrium exists. The initial gas 
temperature is taken to be 273 K. From the theoretical 
density ratio across the shock front and the measured shock 
velocity, the velocity of the gas behind the traveling shock 
is readily computed. The free stream velocity, shown in 
Fig. 6, is almost independent of the initial density, but is 
approximately linearly variant with the shock Mach num- 
ber. The theoretical free stream temperature behind an ideal 
shock is given by Gilmore and is shown in Fig. 7. The fact 
that the actual initial temperature in the shock tube is ap- 
proximately 300 K introduces only a small error when using 
the data of Gilmore. 

Feldman (5) has also computed the properties of air behind 
a traveling shock and also behind a normal standing shock, 
for an initial temperature of 300 K. Unfortunately, the data 
are limited to shock Mach numbers of up to about 22. 

As mentioned previously, it is clear from streak pictures 
that successive and multiple weak shocks occur behind the 
primary shock front, being more pronounced at higher initial 
pressures. Some of these waves are observed to reflect par- 
tially from both the shock front and other traveling waves. 
‘Hence, the state of the gas at some distance behind the inci- 


ARS JourNAL 


rr) 


» 
4 
- 
- 
- 
4 
*t 
‘et 
ad 
10 
| 
|_| 
~ 


Tansi- 
shock 
‘ively, 
y test 


acter- 
th an 
1 par- 
pres- 
f flow 
namic 
nined 

e ini- 
1 that 
nitial 
0 can 
ch as 
nd of 
wave 
were 
bably 
gular 
ream. 
sence 
nably 


shock 


f the 
1 and 
prical 
owed 
times 
hoto- 
that 
3 Or 
nain- 
eady 
and 


ered. 
pres- 
have 
of up 
sume 
3con- 
] gas 
tical 
hock 
hock 
n in 
ut is 
ideal 
fact 
ap- 
Ising 


hind 
ock, 
data 


ures 
the 
itial 
par- 
Wes. 
inci- 


NAL 


15,000 
| 


T 
o? 
w 


3 
N 


TEMPERATURE 


3 
\ 
X 


3 
N 


3,000 
10 15 20 25 30 35 40 


SHOCK MACH NUMBER 


Fig. 7 Theoretical free stream temperature vs. shock Mach 
number for air at several initial density ratios 


dent shock will be different from that in the case of an ideal 
normal shock, due to the additional energy input. 

In order to determine the effect of the successive weak 
shocks on the state of the fluid behind the shock, the electrical 
conductivity of the shocked air was measured according to 
the magnetic technique described by Lin, Resler and Kan- 
trowitz (6). Once the electrical conductivity is measured, 
the temperature of the gas can be estimated from the con- 
ductivity data of Bush (7) and Meyer (8), which are in good 
agreement with the experimental conductivity data of Lamb 
and Lin (9). The experimentally determined temperatures 
are found to be from 400 to 800 K higher than those computed 
for an ideal shock wave. The difference is larger at lower 
shock velocities, as would be expected from study of the 
streak pictures. The free stream temperatures, reaching an 
estimated value of 12,000 K for a Mach 40 shock wave, are 
sufficient for appreciably dissociating and ionizing the gas 
molecules. The stagnation temperatures will be considerably 
higher and may be as high as 25,000 K. The temperature 
and electrical conductivity in the shock layer will be considered 
further in a subsequent section of this paper. 

As a result of these investigations, it is concluded that this 
type of shock tube can produce very high velocities and tem- 
peratures for times sufficiently long to perform certain kinds 
of aerodynamic tests. The principal limitations are that the 
operation of the present assembly is somewhat restricted as 
to the range of initial density, and that the short test times 
make the use of instrumentation by other than optical means 
rather difficult. 

The performance of the shock tube reported in this paper 
is summarized in Table 1. The free stream velocities and 
temperatures far exceed those obtainable in pressure-driven 


shock tubes, and hence, the electromagnetic shock tube ex- - 


tends the regions of possible laboratory exploration. 
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q Table 1 Shock tube performance characteristics 
for air 
shock Mach number, WV; 15 to 40 
initial density ratio, p:/po 10-*to3 10-5 
free stream density, p2/po 3.6 X 10-% to 6.0 
free stream Mach number, M2, 2.9 to 3.9 
free stream velocity, V2 4500 to 12,000 m/sec 
free stream temperature, 7's 6000 to 12,000 K 
test duration 10 to 35 microsec 


A theoretical analysis has been made by Bush (10) to find — 
the effects of applying magnetic fields, on the flow of an ion- 
ized gas about a blunt body of revolution. It was predicted 
that the application of a magnetic field in the region of the , 
stagnation point would increase the standoff distance of the © 
bow shock wave. In addition, the theory predicted that the — 
heat transfer rate to the body at the stagnation point would = 
be reduced when the field was applied. 

The increase in the bow shock standoff distance shiniictbiadg 
correlates with the parameter 


Q* = pol, 


Magneto-Aerodynamic Study 


where 


o; = electrical conductivity of the ionized gas in the shock — 
layer 
By = field strength at the stagnation point on the body 
ry = nose radius 
po, U2, = density and velocity in the free stream, respectively 


The change in standoff distance is given as the ratio of the 
standoff distance with the magnetic field applied to that 
with no magnetic field, i.e., A/Ayy,, and this ratio increases 
with increasing Q*. Proof of the existence of magneto- 
aerodynamic interaction and the quantitative verification of 
the theoretical relationship were the primary objectives of 
the following experimental study. 


Test Model 


The test body, shown in Fig. 8, was a 2-em diam cylinder 
with a hemispherical nose. A magnet coil was located co- 
axially in the nose section and covered by a 1-mm thick Pyrex 
glass nose piece. The energy for the magnet coil was sup- 
plied by a high current pulser which consisted of a high volt- 
age power supply, two 100-uf capacitors and an ignitron 
switch. The field produced by a capacitor discharge is, ad- 
mittedly, nonsteady, but since the test period is very short, 
of the order of 10 microsec, the field will be essentially con- 
stant if the frequency of oscillation is below 5000 cps and if 
the peak of the cycle is used. 

The magnet coil was made of three layers of 10 turns each 
of 20-gage quadruple Formvar copper magnet wire wound 
on a 6.3-mm core, and had a length of 0.95 cm. The coil was 
wound on a plastic core, wrapped with fiber glass tape and 
impregnated with Epon resin for mechanical strength. The 
measured inductance was 6.1 whenries, and since the 
circuit is underdamped, using 200-uf of capacitance, the 
ringing frequency is, theoretically, 4560 cps, and the time to 
the quarter period (very nearly the time of maximum current) 
is 55 microsec. 

Measurements were made of the peak flux density at the 
stagnation point Bo with a single turn pickup coil placed on 
the nose of the model. 8, which was found to be a linear 
function of the capacitor voltage, was 40 kilogauss at 1500 v. 
Operation above 1500 v was avoided to prevent damage to 
the magnet coil. 

In the magneto-aerodynamic analysis (10) which is to be 
verified, it is assumed that the magnetic field distribution 
in the free stream is that of a dipole at the origin (i.e., at the 
center of curvature of the nose radius). This distribution i is 


mathematically expressed by La 
- _ 
B= + 


= ——— & 
r3 2r8 


forr> c. This is similar to the field distribution of a sol- 
enoid. On the axis of symmetry, the field variation for the 
dipole will be, referred to the stagnation point 


B = B(r./r) 


Steady-state measurements were made with a gaussmeter of 
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Fig. 8 Diagram of the test body and agent coil 


pom variation of B along the axis of the solenoid, and found 
to be 

B = Bo(r./r) 3.61 

Hence, the external field along the axis of the magnet coil 

varies approximately as the field along the axis of the dipole, 

and so the magnetic field distribution in the experiment is 

similar to that assumed in the theoretical model. _ 


| 
Shock Tube Tests 


The test body was installed in the electromagnetic shock 
tube, and the apparatus connected as shown schematically in 
Fig. 1, but with the addition of the high current pulser for the 
magnet coil. The magneto-aerodynamic tests consisted of 
measuring the standoff distance of the bow shock wave on 
the test body for variations of the applied magnetic field 
strength and of the free stream flow conditions. The model 
was exposed to steady hypervelocity ionized air flow for pe- 
riods of 15 to 30 microsec, during which the bow shock was 
photographed with an image converter camera. The free 
stream flow conditions were varied by altering the shock tube 
initial pressure, which, in turn, varied the shock velocity, 
free stream velocity, density and temperature. 

The photographic exposure time was 110 millimicrosec, 
and synchronization with the shock tube and magnetic flux 
peak was maintained at each pressure by means of the vari- 
able time delays controlling the events. Thus, at about 10 
microsec after the arrival of the shock front at the model, and 
within the duration of steady flow, the camera exposure 
occurred, coincident with the magnetic flux maximum. The 
fact that the shock standoff distance did not measurably vary 
during the period from the establishment of the bow shock 
to the end of the steady flow period, which was 10 to 20 micro- 
sec after the camera exposure, indicates that the driver gas 
did not yet arrive or that it had essentially the same gas 
properties as the driven gas. 

A series of tests, without the magnetic field and then with 
incrementally increasing field strengths up to 40 kilogauss, 
were made at shock tube initial pressures of 70 and 200 yu. 
The standoff distance was then measured from the photo- 
graphs with a possible error of +8 per cent due to the diffuse 
nature of the shock boundary. 


Table 2 Values of various parameters at the test conditions 


Pr, Uz, T;, 03, oatp/p2U2 
OM; p2/po m/sec K mbho/em  (gauss)~? 
70 21.5 11.9 X 10-4 6710 12,000 55 6.12 K 1078 
200 17.5 28.0 X 10-4 5370 10,300 28 1.66 X 1078 


Fig.9 Typical increase of the bow shock standoff distance upon 
application of a magnetic field 


Experimental Results 


The results of the shock tube tests showed that as the 
strength of the applied magnetic field increased, the bow 
shock standoff distance also increased, proving the existence 
of magnetic and aerodynamic interactions. The effect was 
very pronounced, and at the higher field strengths the shock 
layer thickness increased manyfold, as shown in Fig. 9. 

In order to correlate the experimental data with the theory, 
the values of the factors appearing in the Q* parameter had to 
be evaluated for the test flow conditions. The temperature 
in the free stream was determined from measurements of the 
gas electrical conductivity, and the other thermodynamic 
properties of the air in the free stream were determined 
as described previously. To obtain the conditions in the 
shock layer, the energy equation was applied across the shock 
discontinuity. Since U, > Us, the internal energy in the 
shock layer becomes approximately i 


BE, +(1/2U2 


At the high temperature occurring in the shock layer, the in- 
ternal energy of the air is a function of the temperature and 
density. Therefore, the density in the shock layer must be 
known, and this is found from the relation 

Ps _ ¥3' + 1 


pe 


in which y’ is an effective specific heat ratio and is tabulated 
as y’ = f(T, p) in (4). 

Knowledge of the temperature and density in the shock 
layer permits determination of the electrical conductivity 
from the curves of conductivity for air computed by Meyer 
(8). In this manner, all the thermodynamic factors in the Q* 
parameter are found. Table 2 shows the variation of these 
factors for initial shock tube pressures of 70 and 200 yu. 

The measured standoff distance was normalized as A/ Ay y, 
the ratio of the magnetic to the nonmagnetic case, and is 
plotted vs. the correlating parameter Q* in Fig. 10 superim- 
posed upon the theoretical curve. The data exhibit moderate 
scatter, but the relation between A/Ayy and Q* is evident 


- and not far different from the theoretical prediction. The 


measured standoff distance ratio is greater than predicted 
for both initial pressures over the range tested, and for Q* = 
69, becomes as large as 7.5. 

In addition to inaccuracies in the experimental measure- 
ments and calibrations, and in the values of the computed 
thermodynamic properties of the shocked air, other expla- 
nations exist for the deviation of experiment from theory. 
Specifically, the magnetic field distribution of the solenoid is 
not exactly that of a dipole. Also, the electrical conductivity 
in the free stream, instead of being negligible as assumed in 
the theoretical analysis, is as high as one quarter of the con- 
ductivity in the shock layer. Moreover, whereas theory 
assumed hypersonic flow in the free stream and an incompres- 
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sible flow in the shock layer, in the experiment, hypervelocity 
flow existed, but due to the high free stream temperature the 
upstream Mach number did not exceed 4. As a result, the 
subsonic Mach number in the shock layer was considerably 
higher than theory assumed, and hence compressibility in 
the shock layer was more pronounced. In addition, the con- 
ditions of constant electrical conductivity and negligible heat 
conduction and viscosity in the shock layer are obviously not 
realized. However, the generally good agreement between 
experiment and theory indicates that these discrepancies are 
not severe. 


Conclusions 


The results of the experimental investigation led to the 
following important conclusions. The electromagnetic shock 
tube produces the high velocity and high temperature condi- 
tions required for magneto-aerodynamic research. The bow 
shock standoff distance increases with the application of a 
magnetic field, and the standoff distance ratio correlates 
with the parameter 


a By Properties of Air to 24,000 K,’’ Rand Corp. Rep. RM-1543, Aug. 24, 1955. 
5 Feldman, S., ‘‘Hypersonic Gas Dynamie Charts for Equilibrium 


at least up to Q* = 69, the limit of the experimental data. 
Finally, the experimental data for the standoff distanc 
ratio agree with the theoretical predictions of Bush in spite 
of several discrepancies between the theoretical and labora-— 
tory models. 

The success of the correlation between theory and experi-— 
ment for the bow shock standoff distance suggests a similar 
validity for Bush’s theoretical prediction of a reduction in the | 
heat transfer rate to the body upon the application of a suit- | 
able magnetic field. For the limiting experimental condition 
corresponding to Q* = 69 and A/Ayy = 7.5, theory predicts 
that the stagnation point heat transfer rate should be reduced — 
by about 65 per cent. Experimental measurements of the 
heat transfer rate have not yet been made, but efforts in this — 
direction are being conducted by the author. 


Nomenclature 
B = magnetic flux vector 
c = radius of the bow shock 
ér, eg = unit directional vectors in r and 6 directions, respectively — 
= internal energy 
M; = shock Mach number, based on sonic speed of 332 m/sec — 
p = pressure 
Q* = o3Bo?ry/p2U2, nondimensional parameter 
r = radial distance from center of curvature of nose 7 
7 = temperature 
U velocity 
y' = effective specific heat ratio 
= standoff distance 
p = density 
6 = angle between axis of symmetry and radius vector from _ 
center of curvature of nose 
= unit of pressure, m Hg 
a = electrical conductivity 
Subscripts 
b = body 


= bowshock 

0 = stagnation point 

= ahead of the traveling shock wave 
2 = behind the traveling shock wave 
3 = bow shock layer 
\ nonmagnetic case 
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The Arc Wind Tunnel] 
Tool for Atmospheric 


Entry Research THOMAS R. BROGAN! 


Aveo Research Laboratory 
Everett, Mass. 
— 
2 Because the available test time is limited, questions involving the exposure of materials to the 
= high energy re-entry environment cannot be answered in contemporary quasi-steady sources of 
dpe high temperature air. This paper describes the development and calibration of an arc powered 
wind tunnel designed to provide steady flow simulation of the conditions encountered during 
re-entry. Flight velocities between 16,800 and 21,800 fps can be simulated with this device. The 
flow properties in the tunnel test section can be determined to a degree of accuracy sufficient for 
quantitative experiment. Further work now in progress will improve the degree of simulation 
which can be achieved with are powered wind tunnels. 


= 


HEN aa satellite or ballistic missile re-enters the at- to establish steady ablation can be estimated using a simple 
mosphere, its surface is exposed to air at temperatures analysis. As shown in Fig. 1, an infinite half plane of material 
which may be as high as 8000 K, with stagnation pressures as with specific heat C,, density p, and thermal conductivity K 
high as 50 atm. The simulation of such an aerodynamic en- is receiving a heat input q at its open face. The material has 
vironment is an interesting task for the designer of research a heat of ablation Q*. The surface of the ablating material is 
- tools. In the first place, the methods by which air at the in- stationary, so that the material as a whole moves toward the 
teresting temperatures can be produced are limited in number, surface with a velocity v., given by 
and, second, when the requirement of a very high gas tem- a 
perature is combined with the demand for some sort of scale vo = g/Q [1] 
simulation, one finds that the necessary power inputs to a In steady ablation, the temperature distribution in the ma- 
steady flow test device can easily become prohibitive. terial is described by wow © 
This difficulty with the power input has led to the develop- a -_ 
ment of quasi-steady methods for producing the high energy ay’ ay = : 
air at the required flow rates. With these methods, a large 
. amount of stored energy is suddenly released to produce high - 
enthalpy gas for a relatively short period. In spite of the w hich i the solution a) 
limited test time, it is possible to establish steady flow for a 7 a es = a 
sufficient period to obtain quantitative data. exp ( = 5 = [3] 
The shock tube is a classic example of such a quasi-steady 7 Ts 6 T. PlaCp - 


research tool. It is simple, but extremely versatile, and can 
be operated over a large range of all the interesting parame- 
ters. Shock tube studies have contributed immensely to an 
understanding of the re-entry problem, particularly in the 
field of high temperature gasdynamics. Other devices em- 
ployed with varying degrees of success include shock tunnels, 
ballistic ranges and the “hot shot” tunnel. The pebble bed 
heater has been used for high speed flight research, but it can- 
not produce gas at the enthalpies appropriate to missile or 
satellite re-entry. 

Although all of the devices mentioned have added to our 
knowledge of hypersonic flight, they are not adequate because 
of their limited test time when questions involving the ex- 
_ posure of materials to the high energy environment must be 
answered. Such problems must be studied over relatively 
long test times, and this becomes particularly important when 
the re-entry body has an ablating surface. For this applica- 
tion, a research device with a relatively long test time is 

needed. 


The order of magnitude of the minimum required test time 


where 7’, is the surface temperature, and 7) the temperature 
atz = ©,W hich may be assumed small. 


Presented at the ARS 13th Annual Meeting, Nov. 17-21, 1958, 

1 Researc ientist. 


Fig. 1 Temperature distribution in steady ablation 
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Table 1 Time to reach steady ablation 


Ablating material 


Pyrex 


to, sec 
0.001 
Teflon 0.0005 
carbon-pure graphite brick 0.8 
carbon-amorphous tthe 0.03 


Now to establish steady ablation, the thermal diffusion 


depth Vv Kto/ pC, must be of the order of 6, where fo is a 
me:sure of the time to reach steady ablation, or 


K 
4 
pC, ( q | 


For times less than to, the material behaves primarily as a 
heat sink, and not as an ablating body. In order for an 
ablation test to be meaningful, the test duration should be 
several times the value of fo. 

The maximum test time in a shock tube or shock tunnel is 
a few hundred microseconds, while in the ballistic range or hot 
shot tunnel, the test time is about 1 millisee (1,2).2 Assume 
that an 0.5-in. diameter sphere of the material to be tested 
is fired into atmospheric air at a speed equivalent to a Mach 
number of 15. At the stagnation point (which may be shifting 
due to rotation) the heat transfer rate is about 27,000 w/em? 
(6). A similar heat transfer rate is to be expected in the hot 
shot tunnel (1). Assume then that the quasi-steady test de- 
vices deliver a heat transfer rate of 27,000 w/cm? for a period 
of 1 millisec. 

Under this test condition, consider an ablation test on three 
typical materials: Pyrex, Teflon and carbon. For Pyrex, a 
typical value of Q* is 610 cal/gm (3), and for Teflon (which 
sublimes), 2800 cal/gm (4). For carbon, the value of Q* may 
be estimated as being equal to the heat of vaporization or 
about 5600 cal/gm. The values of to for these materials are 
given in Table 1. 

Since the test duration should be several times the value of 
lo, it is evident from the data in Table 1 that the quasi-steady 
test devices with a test time of 1 millisec-are unsatisfactory, or 
at best very marginal for ablation ‘studies. On the other 
hand, while the heat transfer rates which can be produced in 
an electric are wind tunnel are much less than those available 
in the quasi-steady devices, the test time in the are wind tun- 
nel is unlimited. 

The common electric arc is used for countless applications 
where high gas temperatures are required. As an example, 
in the ordinary welding arc in air at l-atm pressure, the gas 
temperature in the core of the arc is in the vicinity of 6500 K, a 
temperature which would exist at the stagnation point of a 
blunt body flying at 22,000 fps. The are is an attractive 
source of high temperature gas, and in the past few years 
attention has been centered on its application in a hypersonic 
aerodynamic research device. 

Much of the interest was stimulated originally by the work 
of the German group under Maecker (5), who, in connection 
with their extensive research into the mechanics of ares, built 
the first examples of the type of apparatus which has become 
known as a “plasma jet.’’ A diagram of the original apparatus 
is shown in Fig. 2. Water is swirled about a carbon cathode. 
This rotation ‘‘stabilizes’’ the arc, that is, confines the arc to 
strike only inside the carbon anode into which a nozzle is 
machined. The hot plasma (water and carbon vapor) is 
ejected through the nozzle. The arc was operated at pres- 
sures as high as 50 atm. Using spectroscopic methods, the 
gas temperature at the nozzle exit was found to be 8000 K. 
By similar methods employing Doppler shift, the average 
gas velocity was measured at 6500 mps. The extensive arc 
research of this German group included detailed spectroscopic 


? Numbers in parentheses indicate References at end of paper. 
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analysis, determination of gas velocities in the arc, and 
studies of the effect of a magnetic field on the arc stability. 
They were not primarily concerned with the application of 
the are as a gas source for a hypersonic research device. 

This paper, however, is concerned specifically with just 
such an application. Indeed, we will neglect, insofar as is 
possible, the mechanics of the are itself and concentrate on 
using the arc as a source of high enthalpy air for a wind tunnel. 
Since the arc powered device is of primary use in the study of 
ablation, we will be concerned mainly with the development 
of facilities, which provide adequate simulation for phenomena 
associated with ablation. It is sufficient to remark here that, 
in any ablation test, the flight enthalpy is the most important 
simulation parameter and must be duplicated in the test de- 
vice (3,4). 

The paper begins with a description of a low density arc 
wind tunnel built at the Avco Research Laboratory in order — 
to study the feasibility and operation of arc powered re-entry _ 
test devices. Next, the calibration of this tunnel is described, 
with particular attention to the determination of the gas en- 
thalpy. Finally, there is a short discussion of further arc — 
wind tunnel development which is now in progress. 

Low Density Arc Wind Tunnel _ 

In order to study the feasibility and operation of arc 
powered test devices, a low density are wind tunnel was con-— 
structed and put into operation at the Avco Research Labora- 
tory. Development of this low density are wind tunnel 
was begun using a configuration similar to that shown in Fig. 
2, but with two important exceptions. First, air was used as 
the working fluid instead of water, and, second, the carbon _ 
anode was replaced by a forced convection, water cooled — 
copper anode, which forms a nozzle for the gas to be drawn 
from the arc. Since the major fraction of the carbon impurity 
in the German arcs was due to erosion of the carbon anode, 
we expected that substitution of a water cooled copper anode _ 
would appreciably reduce the impurity in the gas produced by 
the arc. This has proved to be the case. 

In an attempt to reduce the impurity still further, a water 
cooled metal cathode was substituted for the carbon cathode. — 
This was unsuccessful; the metal cathodes led to very large 
impurity levels. Turning again to carbon cathodes, we at- 
tempted to remove the carbon impurity by applying vacuum 
back through the cathode itself, but this suction had a | 
drastically unfavorable effect on the arc stability and again — 
was unsuccessful. Finally, we settled on 0.5-in. diameter _ 
National Carbon “Intensare” welding electrodes as giving the | 
minimum are contamination. At the highest power input to 
the tunnel (130 kw), there is a 2 per cent carbon impurity in 
the gas produced by the are. At lower power inputs, the con- 
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tamination is much less, about 0.4 per cent at the minimum 
power input (50 kw). 

After establishing a stable are configuration, we proceeded 
to a preliminary examination of the gas produced by the arc. 
Two serious defects were noted. First, the intensity of the 
light emitted from the gas flow varied in time by a factor of 
25 or more, and at a rate corresponding to the frequency of 
oscillations in the power input to the are. Second, the test 
section gas properties were obviously nonuniform in space; 
the presence of a “hot core’ ’ composed of gas from the center 
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Fig. 4 The low density arc tunnel 


_ Fig. 5 Five-sixteenth-in. diameter hemisphere cylinder in 
1.033-in. nozzle 


of the are was clearly evident. The nonuniformity in both 
time and space rendered the initial arc configuration un- 
suitable as a gas source for a high enthalpy wind tunnel. 

To reduce the nonuniformity, it was decided to use a wind 
tunnel with a two-throat system to isolate the arc from the 
test section. A sketch of our final arc tunnel configuration is 
shown in Fig. 3, and a photograph in Fig. 4. The are burns 
down inside the first nozzle and is stabilized by air injected 
tangentially at the nozzle entrance. The gas passes from the 
are into a settling or plenum chamber, which serves to damp 
both the space and time variations in the gas properties. [In 
order to damp the time variations, the volume of the plenum 
chamber must be large in comparison with the volume of 
gas which passes through the first throat during one oscillation 
in the power input to the are. To reduce the space varia- 
tions, there must be enough dissipation in the plenum chamber 
to provide good mixing. However, the dissipation must not 
be so great that the heat loss to the wall becomes prohibitive. 

The amount of dissipation in the plenum chamber is regu- 
lated by the ratio of the first and second throat areas. After 
some experimentation with the plenum chamber size and 
shape, the final configuration consisted of a spherical plenum 
chamber of 3-in. diameter, with a first throat 0.6 in. in diameter 
and a second (sonic) throat 0.3 in. in diameter. With this 
configuration, it appears that about 10 to 15 per cent of the 
energy imparted to the gas in the arc is lost in the plenum 
chamber. 

Although we have not carried out a complete study of the 
effects of the plenum chamber on the uniformity of the gas 
properties in the test section, it is obvious that the uniformity 
is much improved. For instance, examination of the light 
emitted from the gas in the test section shows that the time 
variations are reduced from a factor of 25 to a factor of two or 
three. Also, the light intensity variation is concentrated in 
the CN impurity bands, so that it is not certain that the time 
variations in light intensity represent actual differences in 
the gas properties as a whole. With regard to the space 
variations, we have conducted only visual examinations. 
Using the plenum chamber, there is no evidence of a “hot 
core” in the test section, as was clearly evident before the 
plenum chamber was incorporated. 

At the second throat (Fig. 3) the flow is sonic, and the 
Mach number increases as the flow passes out the second 
nozzle. These nozzles are designed to produce a slow expan- 
sion in the vicinity of the throat, but no attempt was made to 
use characteristics in the nozzle design. Nozzles up to 2-in. 
exit diameter can be used free of shock waves; with the 
maximum exit diameter, the flow Mach number at the nozzle 
exit is about 4. 

The nozzle discharges into a test section 6 in. in diameter. 
Tests are conducted } in. downstream of the nozzle exit. The 
tests may be observed through quartz windows installed in 
the test section. Models are mounted on a water cooled 
strut. 

The gas from the test section passes into a water cooled, 
staggered, cross-bank heat exchanger, where it is cooled very 
nearly to ambient temperature. It is then fed to a vacuum 
pump of sufficient speed (500 ft?/min) to maintain the low pres- 
sure required in the test section. 

Fig. 5 illustrates the supersonic flow about a 3-in. diameter 
copper hemisphere cylinder in the are tunnel. The photo- 
graph is taken with a nozzle of 1.033-in. exit diameter, using 
the light emitted by the gas itself. - 


Calibration of the Low Density Arc Tunnel 


This section is concerned with the determination of the 
flow properties of the gas produced in the arc tunnel. The 
most important quantity to be determined is the stagnation 
enthalpy of the gas in the test section. This parameter 
is of highest importance in all ablation studies (3). Our 
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efforts have been directed toward determination of this 
quantity to a degree of accuracy sufficient to permit quantita- 
tive interpretation of data obtained in the tunnel. 

At the exit of the supersonic nozzle (Fig. 3), two quantities 
are measured. These are the stagnation pressure p, and the 
average heat transfer rate to blunt, axisymmetric bodies. 
The stagnation pressure is measured with a water cooled 
pitot tube. The heat transfer is measured by placing a metal 
model in the nozzle exhaust and noting the time which elapses 
before the model begins to melt. Applying the theory of 
laminar heat transfer to blunt bodies of revolution, the gas 
enthalpy can be determined from these two measurements. 

Assume a constant, uniform heat input g into an infinite 
half plane of material of density p, constant specific heat C, 
and thermal conductivity K. The material is initially at a 
uniform temperature 7, and melts at a temperature 7’,,. 
Then the heat transfer rate is given by 


1 
q= (T, — To) 5] 


where é is the time required to melt the surface. A long 
cylindrical rod with insulated sides subjected to nonuniform 
heating at the free end will behave as if it is receiving an 
average heat input given by 


where 
= local heating rate 


q 
R = radius of the rod 
{ = surface area at the free end 


Equation [6] is valid if, and only if, the thermal diffusion 
depth in the melting time ¢; is much greater than the radius 
R, that is 


>R 


For the copper and aluminum models used in the calibrations, 
the thermal diffusion depth is always greater than 10 times 
the body radius. These remarks also apply to the case of a 
long hemisphere cylinder with insulated sides. Since the 
temperature gradients in the hemispherical portion of the 
body are small, the body behaves as if it were flat-faced and 
receiving an average heat input given by Equation [6]. In 
this case, the integration is performed over the hemispherical 
nose. 

A sketch of a typical heat transfer calibration experiment 
is shown in Fig. 6. The supersonic nozzle has an exit diameter 
of 1.033 in. and a throat diameter of 0.300 in.; the nozzle 
area ratio is 10.9, and the flow Mach number at the exit is 
about 3. The nozzle contour is designed to produce a slow 
expansion in the vicinity of the throat; the exit is a point of 
inflection where the half angle is 8 deg. The models are 
placed + in. downstream from the nozzle exit. A flat-faced 
$-in. diameter model is shown, but other diameters and hemi- 
spherical models are also used. In all cases a fiber glass or 
Astrolite shield prevents heat transfer to the side of the body. 
To avoid heat loss from the model to the shield, there is a 
0.015-in. air gap between the shield and the body. The lead- 
ing edge of the shield is swept back so as not to influence the 
apparent nose radius of the body. 

Depending on the enthalpy at which the tunnel is operated, 
either copper or aluminum models are used for the heat trans- 
fer determination. For both of these materials, the assump- 
tion of constant thermal conductivity over the whole tem- 
perature range is reasonable, since the variation in the con- 
ductivity is small, i.e., 25 per cent for copper. An average 
value of K is used in computations. 

The melting time is measured either from motion pictures, 
or with a stopwatch. Melting is very sharp and uniform, with 
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Fig. 7 Stagnation point heat transfer to a hemisphere cylinder — 


the entire model surface melting at almost the same time. — 
This attests to the fact that the calibration experiment does, | 
indeed, average the heat input over the face of the model. 

The stagnation point heat transfer to blunt bodies of revo- 
lution in dissociated air has been studied theoretically by 
Fay and Riddell (6), and experimentally by Rose and Stark | 
(7). Theory and experiment are in excellent agreement. For 
present purposes, the results are summarized in Fig. 7; 9,0 
is the heat transfer at the stagnation point of a hemisphere 
cylinder; F is the nose radius; H is the stagnation og oe 
® the gas constant, and p, the stagnation pressure. T» is 
room temperature (293 K). Note that for given tunnel con- — 
ditions, the product 9g,0V R is a constant. 

In the case of the flat-faced bodies, the stagnation point — 
heat transfer cannot be predicted theoretically since the nose — 
radius is infinite. Measurements by Stoney and Markley (8) — 
show that g, = (0.55 + 0.05)q,0, where q, is the stagnation _ 
point heat transfer on the flat-faced cylinder of radius R. 

The laminar heat transfer distribution over a hemisphere 
cylinder has been studied by Kemp, Rose and Detra (9). 
Using their distribution and performing the integration indi-— 
cated by Equation [6], we find that for the hemisphere cylin- | 
der g = 0.92 go. For the flat-faced bodies, from the studies _ 
of Lees (10) and the experiments of Stoney and Markley (8), 
g = 0.669¢,0. Remember that 9,0 is the stagnation point heat — 
transfer to a hemisphere cylinder of radius R. 

Since g is proportional to 1/V R, and the melting time is" 
proportional to (1/g)?, the melting time should be directly — 
proportional to the model size with tunnel conditions fixed. 
Further, for the same model size, the ratio of the melting times — 
for the flat-faced cylinders to that of the hemisphere cylinders — 
should be (0.66/0.92)? or 0.52. In Fig. 8 the melting time is 
plotted vs. body diameter for both shapes from a typical 
calibration. A straight line is drawn through the data for 
each shape, and a value of gGVR is determined using Equation 
[5]. The value of GoVR is then known, and the enthalpy 
can be determined from Fig. 7. 

The measured melting times are in the ratio 0.57, as com- 
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pared with the predicted value of 0.52. This discrepancy 
corresponds to a 5 per cent difference in enthalpy. The scat- 
ter in the melting time data represents a further 8 per cent 
uncertainty in the results. By averaging the results between 
the two shapes, we arrive finally at an overall uncertainty of 
10 per cent in the determination of the stagnation enthalpy. 

At the low gas densities in the tunnel test section, the gas 
may be out of thermodynamic equilibrium. However, as 
shown by Fay and Riddell (6), the heat transfer is changed 
only a negligible amount if the gas is out of equilibrium, pro- 
vided that H is the total enthalpy, however distributed, and 
the model surface is catalytic. Since the surface temperature 
is never above 1350 K, it is very likely that the wall is cata- 
lytic. Hence, lack of equilibrium does not influence the 
calibration. 

It is interesting to consider the stream area over which the 
calibration averages the gas properties. With the model sizes 
and gas densities in the arc tunnel, about one tenth of the 
stream energy incident on the body is transferred to the body, 
but the heat which goes into the body comes not from the gas 
incident on the body as a whole, but only from that gas which 
enters the boundary layer. The boundary layer is composed 
of that gas which passes through the bow shock in the im- 
mediate vicinity of the stagnation point. Therefore, the 
heat transfer measurements represent an average determina- 
tion of the flow properties over a stream area of approximately 
one tenth the frontal area of the body used in the calibration. 
For the largest models (3g in.), the gas properties are averaged 
over a stream diameter of 0.1 in. 

For practical purposes, the tunnel has been calibrated at 
three different operating conditions corresponding to stagna- 
tion enthalpies H/RT, = 178, 217 and 263. These values of 
enthalpy correspond in turn to flight velocities between 16,800 
and 21,800 fps. 


Further Development of Arc Powered Wind 24 
Tunnels 


Although the are wind tunnel described accurately simu- 
lates important re-entry enthalpies, other parameters, such 
as scale, i.e., Reynolds number are not simulated. Further 
work directed toward providing scale as well as enthalpy 
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Fig. 8 Melting time vs. diameter, copper models 1 . 


simulation is now in progress at the Avco Research Labora- 
tory. This work involves a large increase in are power level. 

The first step in the direction of increased power has been 
the construction of a 2000-kw tunnel which is now in opera- 
tion. It is similar in basic design to the tunnel previously 
described, and operates at stagnation pressures up to 20 atm. 

The heat transfer rate at the anode of a high intensity arc 
is extremely large. The water cooled copper anodes used in 
the wind tunnels can be made to handle average heat transfer 
rates of up to about 5 kw/cm?. It is our experience that a 
power of about 2000 kw is the maximum that can be handl--d 
by our are configuration without exceeding the critical heat 
transfer rate to the arc anode. Any increase in power level 
above 2000 kw must involve an extensive modification of thie 
tunnel design. 

There is a simple way to get around the power input limit 
to a single arc. The method involves the simultaneous oper:- 
tion of several arcs, each at a power level of 2000 kw, and all 
discharging into a single plenum chamber. This configuration 
has other advantages: Excellent plenum chamber mixing, 
lower relative plenum chamber heat loss and the opportunity 
to remove the vorticity due to the air stabilization by proper 
geometrical arrangement of the individual ares. Such an arc 
wind tunnel is now under development at the Avco Research 
Laboratory. The total power input of 10,000 kw (supplied 
from a battery bank) is divided equally among five separate 
arcs. The tunnel is designed to operate at a pressure of 20 
atm. 


Conclusions 


Questions involving the exposure of materials to the high 
energy re-entry environment cannot be answered in the present 
quasi-steady test devices. The electric are wind tunnel pro- 
vides a steady flow of high enthalpy air, and can be calibrated 
to a degree of accuracy sufficient for quantitative re-entry ex- 
periments. Continued development of are powered wind 
tunnels will provide simulation of re-entry scale parameters in 
addition to simulation of flight enthalpy. 
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Nomenclature 

C, = specific heat, cal/gm K 

H = total gas enthalpy, cal/gm 

K = thermal conductivity, cal/em sec K i 

= stagnation pressure on a blunt body, psia 

= atmospheric pressure, 14.7 psia 

q = heat transfer rate, w/cm? 

dso = heat transfer rate at the stagnation point of a hemisphere 

cylinder 
qs = heat transfer rate at the stagnation point of a flat-faced 
cylinder 

G = average heat transfer rate to a blunt body - Vv = 

Q* = heat of ablation, cal/gm 

R = body nose or cylinder radius, em pret 

T = temperature, K + * 

T. = body surface temperature 

Tm = material melting point 

t) = time to reach steady ablation, sec 

i; = time to melt the body surface 

v. = ablation velocity, cm/sec 

xz = coordinate perpendicular to the body surface, em 

G = the gas constant for air, 0.0686 cal/gm K 

p = material density, gm/cm’ 
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“Heat Conduction in a Circular 
Cylinder Rotating With Its Axis 


Normal to the Direction of Flight 


ANTHONY DEMETRIADES! | 


California Institute of Technology 


Pasadena, Calif. 


The time-temperature history of an infinite solid cylinder of circular cross section is examined 
when heat is transferred to it in a manner resembling heat transfer in high speed flight, and — 
the cylinder is simultaneously rotating about its axis. The heat conduction equation is solved for 

_ boundary conditions peculiar to such heat transfer, and it is shown that the pertinent parameter _ 
in the solution is the ratio of the thermal wave length to the cylinder radius. For small values of — 
this parameter (corresponding to angular speeds above a few radians per second) it is found that x 


the cylinder rotation about its axis. 


HE PROBLEM of heat transfer to a solid body moving 

at hypersonic speeds through a frictional atmosphere has 
been the subject of extensive investigations. Attention is 
usually concentrated on the forward stagnation point of the 
body, since it is there that the most severe heat transfer rates 
occur. It is commonly known that these stagnation point 
heat transfer rates are often too severe to be withstood by 
current structural materials. 

This paper investigates the phenomenon occurring when 
the body is given a rotation about an axis normal to the direc- 
tion of flight. More specifically, we will investigate the time- 
temperature history of an infinitely long circular cylinder with 
its axis normal to that of flight, the cylinder also rotating 
about the same axis. This rotation will normally result in a 
number of different effects; for example, it will induce an 
aerodynamic circulation, may produce some lift, and it may 
produce centrifugal stresses in the body material of a severity 
related to the angular speed of rotation. We will ignore these 
effects in favor of concentrating on the problem of heat con- 
duction within the material itself. 


General Considerations 


The equation of heat conduction 


= DV?7' ® 
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the cylinder maximum surface temperatures are reduced by a considerable amount by virtue of 


since the surface (or boundary) heat transfer to a body is 
- proportional to a lower power of its dimension # rather than 


Transfer to a Rotating Cylinder 


- eylinder of radius R and infinite length rotates about its axis 


r.. 


can be nondimensionalized with the aid of a characteristic — 
temperature 7o, length R and speed of rotation w. Thus— 


Awt) wk? To 


The quantity D/w has the dimensions of length squared and — 
is the square of the thermal wave length d aan 


(7) (7 
We therefore see that if the thermal wave length is small com- — 
pared to the characteristic body dimension, the time rate of 
temperature change is small. 

Since the thermal diffusivity © is a characteristic of the | 
body material only, the above considerations lead to the con- 
clusion that for a given angular speed w(which appears in the | 
boundary conditions) the time rate of change of temperature — 
depends on the absolute magnitude of the scale factor R. 
This is, of course, directly evident from physical reasoning. 


overall heat capacity. 


We shall now solve the following problem: A circular 


which is normal to the direction of flight. Heat is transferred 
to the cylinder from the flow about it. We seek the cylinder 
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First, write =... 
+ Ts 


so that if we rewrite pawn 


( a+ e + 3 e 


we have 


9 


The solution is carried out with the aid of Laplace transform~.* 
Thus, 7; is 


2 


Fig. 1 Above: Cross section of cylinder showing the body ~ 4 
motion, coordinate system and angular rotation. Below: As- where aj, Q, a, ete., are given by the solution of = 


of Qs) 


sumed distribution of surface heat transfer aD 


Jia.) = 0 


The expressions for 7’, and 7; are complicated by the fact 
that, unlike 7,, they depend on the coordinate @ 


—iw — B, 2% 2 
J; r Ji Bs *Dt/R? 
2 d B2D\[ d d [— 
dk D R? p= —6sD/R? 
= 
= = [7] 
d iw s-1 dd 
dR dp dR p= —Bs*D/R? 


The denominators of these expressions are given by 


temperature asa function of the radius 0 < r< R, the angle @ d fico lio 1 ‘ica 
and the time as in Fig. 1. J; r) Vo Jo r) J; 


The equation of heat conduction is 


Having coordinates chosen on the cylinder we can now imagine — De ees = 


[6] 


Nin 


that the stagnation line rotates around it with an angular ' 4 
speed w. Further, for the distribution of the heat transfer er Pee va D 
rate around the cylinder circumference, we can use the laminar ; 

heat transfer rate distribution given, for instance, by in _ Since 7 is the complex conjugate of T2, we write 


If we approximate this distribution by a cosine function, we 
can finally write T2 +7: =b | 


e(ac + bd) + fiad — bc) tse lag + 4 


(OT /dr);-r = a + b cos (6 + wt) [2] 
where a and 6 are defined in Fig. 1. We can also write where * a £40 
T = To = constant [3] = cos fo, , [9a 
at t = 0. These two relations serve as boundary conditions 1 ri 


for the problem. Condition [2] is somewhat conservative in 
that for a finite radius R, its right-hand side should include a Red; Vi)eos wt + ImJ, (- wt 


nonharmonic time factor; this omission, however, serves to — 


simplify the problem without concealing its salient features. [9e] 

2 Lees, L., ‘‘Laminar Heat Transfer over Blunt-Nosed Bodies ’ The author is indebted to Dr. E. Hellund for his cooperation 
at Hypersonic Flight Speeds,” Jer Proputston, April 1956, pp. in discussing the mathematical aspects of this and similar prob- 
259-269. ems. 
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Fig. 2 Contribution 7’; to the cylinder stagnation point tempera- 
orms.* ture. a is defined in Fig. 1 
= ( Vi) sina +ImJ; Vi os wt 
(9d) 


4 =~ — 1) [9g] 
| = B2/R? [9h] 
(9i] 

| where Bs for s = 1,2,3,...are given by 
Jo(Bs) = Ji(Bs)/Be [10] 
| We can now proceed to interpret the contributions 7), T» 


and 7’; to the body temperature. The term 7), as given by 
Equation [4], possesses angular symmetry and is thus inde- 
pendent of 6, hence also of w. Of the contribution of 7, + 

e) T;, the first term in the bracket of Equation [8] is purely 
harmonic with wt; the coefficients of the trigonometric func- 

& tions are functions of the ratios \/r and \/R as shown by Equa- 
— tions [9] and as expected, of course, from the preliminary 
remarks in this paper. The second term is also a function 

of the ratio X/R. 

. A numerical example has been computed using the analyti- 
| cal results obtained. We considered the case of a solid steel 
circular cylinder (D = 0.07 em?/sec) with a 20-cm radius 

rotating with angular speeds ranging up to 100 radians/sec. 
Although more generalized calculations using \/R could be 
also made, the numerical example as computed is of direct 
illustrative value. The results are shown in the accompany- 

(8) ing graphs. In these, the computed temperatures are at the 
forward stagnation point (r = R, 6 = 0) of the body in 
terms of its own coordinate system. The temperatures them- 
[9a] selves are given in terms of a and b as taken from Fig. 1. 
Direct substitution of the numerical values of a and b (in 


ad cgs units) for specific flight conditions will yield the 
numerical value of the temperatures. 

wl lig. 2 gives the time history of the term 7, and thus shows 
the lowest stagnation point temperature attainable notwith- 

[9e] stunding the magnitude of w. Fig. 3 shows the term 7, + 7; 

ion for different values of w. It should be noted that each 

= curve in this figure includes the harmonic envelope of the 


first term in the brackets in [8], and thereby denotes the 
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[5] = 25 Fal (2 (2 vi) 
f=-- b —her|R ImJ,(R Of 


= 
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Fig. 3 Maximum values of the contribution 7. + 7’; for various — 
values of the speed of rotation 


maximum value of 7, + 7’; for that particular value of w. 
It is evident that for angular speeds of the order of 1 


radian/sec or more, the term 7, + 73 becomes negligible; 
this condition corresponds to having a thermal wave length 
smaller than about one-hundredth the cylinder radius. — 
Under these circumstances, the stagnation point (or, more 
accurately, the stagnation line) temperature of the cylinder 
will be reduced to as much as half the value which an identical 
cylinder would exhibit under identical (high speed flight) con- _ 
ditions, but without rotation. Such results, however, are 
applicable only as long as the heat transfer rate relation [2], 
based on a quasi-steady boundary layer around the cylinder, 
is not violated. Roughly speaking, this requires that the — 
local flow velocity in the inviscid field around the cylinder be 
much larger than the surface velocity wR produced by the 
rotation, 

It would be of interest to extend these considerations to the 
case of a slender body of revolution rotating about its axis 
while also moving through the air at high speeds and at high 
angles of attack. The problem of the rotating sphere is also 
under consideration by the author. 


Nomenclature 

= temperature 

t = time 

D = diffusivity = k/pc 


k = thermal conductivity of material 
p = density of material 

c = specific heat of material 

To = initial or reference temperature 
R = cylinder radius 

w = angular speed of rotation 

r = thermal wave length = (D/w)'/? 
r = radial distance from cylinder axis 
6 = meridian angle 

a,b = constants defined by Equation [2] and Fig. L 
as = solution of Equation [5] 

Jo, J: = Bessel functions 

3 = solution of Equation [10] 

Re... .) real part of terms following 


Im(...) = imaginary part of terms following ' 
ber, bei = the Bessel functions defined by Jo(zVi) = 
i ber x 


ber x — 


655 


| 
| 
| | 
| 4 | 
| u)= 
| ii 0.01§ 
| 
| j 
ad | 
o.! 
i | 
e 
| 


phasized. 


N SPACE, man will be completely cubis on 1 the space 
vehicle and its systems to sustain life. He cannot long 
exist outside the sealed cabin nor can he return quickly to 
the life-supporting regions of Earth’s atmosphere in case of 
malfunction or failure. As the limits of outer space are ap- 
proached, the probability of successful escape would appear 
to become smaller and smaller; the man becomes more and 


more dependent upon the machine for survival until at some - 


point the fate of the two are, to all intents and purposes, one 
and the same. 

The prime consideration in this matter is the overall reli- 
ability of the manned, sealed cabin system. It is evident 
that the system must be so conceived that the probability 
of failure of both the human and the machine approaches 
zero. The physiological and psychological stresses on the 
individual must be reduced to levels which are well below 
maximum tolerance when applied either individually or 
collectively over extended periods of time. The sealed cabin 
and its internal systems must be executed in accordance with 
the concept that no single failure can be catastrophic and, 
further, that the effects of any predictable hazardous occur- 
rence shall be mitigated to survivance levels. The most 
reliable sealed cabin is one which provides: 

1 Optimum environment to permit maximum human 
performance. 

Maximum structural reliability. 
Maximum protection for occupants. 
Maximum reliability of internal systems. 
Minimum fire hazard. 


Orde 


Human Performance 


During the Symposium on the Physics and Medicine of the 
Upper Atmosphere, given at the USAF School of Aviation 
Medicine in September 1951, Ulrich Luft (1)* stated that, 
“From a physiological point of view, it would appear desir- 
able to onpely the aviator with a cabin environment which 

Petal at the ARS 13th Annual Meeting, New York, 
N. Y., Nov. 17-21, 1958. Revised April 1959. 

3 Chief, Interiors Design Section. 
? Assistant Chief, Strength and Dynamics Section. 
3 Numbers in parentheses indicate References at end of paper. 
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This paper ounsiders the “ie antages to the crew of an optimum internal atmosphere and pres- 
sure and the effects on a sealed cabin structure. 
fast fracture are discussed and shown to be essentially constant for various internal pressures with 
proper design. Certain safety advantages of a cabin designed for high pressure are pointed out, 
for example, a greater resistance to meteoroid penetration and a longer decompression time should 
a penetration occur. Means of increasing the time of decompression to permit effective emer- 
gency action are discussed, and the need for complete reliability of the sealed cabin system is em- 


The reliability of the cabin and its resistance to 


resembles the conditions on the earth’s surface as closely us 
possible. This would certainly ensure an optimum range of 
performance and endurance. However, the chances are 
that it may not be possible to adhere to this ideal situation 
for various, mainly technical, reasons. It is necessary, 
therefore, to anticipate certain limitations of the environment 
which may have to be imposed for the duration of flight and 
also to consider marginal conditions arising in emergencies.” 

The “limitations of the environment” mentioned by Luft 
refer primarily to conditions of the cabin atmosphere, i.e., 
the total pressure in the cabin and the partial pressures of 
the constituent gases. The technical reasons for these en- 
vironmental limitations were very real in 1951; indeed they 
exist today in conventional pressurized cabins which rely 
upon the external atmosphere as the principal gas supply. 
A sealed cabin, however, is entirely independent of this at- 
mosphere. It offers complete isolation from the external 
environment and has the potential capability of maintaining 
an internal sea-level atmosphere which will satisfy all of 
man’s known physiological requirements and permit him to 
function at maximum efficiency for extended periods. Shirt- 
sleeves can become the normal flight clothing in sealed cabins 
under such conditions. In terms of human performance, the 
advantages of a sea-level atmosphere have been clearly 
demonstrated by the experiences of Ross and Lewis (16) 
during the recent Strato-Lab High 2 and 3 flights. 


Structural Reliability 


The greatest possible assurance of the reliability of a 
structure lies in the ability to predict accurately the behavior 
of the structure and its materials under thé actual conditions 
of loading. 


Resistance to Fast Fracture and Rapid Decompression 


The state of the art for pressure cabin design has progressed 
rapidly in the past few years. The cause of uncontrolled 
fast fracture of the cabin shell under cyclic loadings and the 
resultant effects of rapid decompression are reasonably well- 
known. Recent investigations by Christensen (2,3), Irwin 
(4), Sorensen (5) and others into the mechanism of cracking 
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and the fracturing characteristics of certain materials have 
established a firm basis for predicting and preventing fast 
fracture. 

The reliability of a pressurized structure can be assured 
for any design pressure differential through proper design 
techniques, the utilization of suitable fatigue and tear re- 
sistant materials operating at safe stress levels, and a suitable 
test program. For example, a representative section of a 
jet transport fuselage has been subjected recently to simulated 
service loading tests in which over 140,000 complex loading 
eycles, including a 9.3-psi pressure differential, were imposed 
on the specimen (6). This is the equivalent of 40 to 120 
years of endurance, depending upon the type of operations 
encountered by a given aircraft. 

The tear resistance of a material is a measure of its ability 
to resist the propagation of a crack through the material. 
Under cyclic stress, fatigue cracks are believed to form very 
early in the cycling process. Safe, low stress levels may 
cause little or no propagation of the crack after its formation. 
Under higher, unsafe stresses, however, these cracks propagate 
slowly through the material and, under certain conditions, can 
attain critical length. At this point, the elastic energy being 
stored within the shell becomes equal to the critical tear re- 
sistance of the material; the stored energy is released, and the 
crack suddenly progresses through the material in a very 
rapid and uncontrolled manner. Irwin (4) has referred to 
the elastic energy accumulated and released at fast fracture 
as the “crack extension force.” 

Critical tear resistance varies over a wide range for different 
structural materials and frequently bears little or no relation 
to the ultimate tensile strength of the material. Fig. 1 indi- 
cates typical values of critical tear resistance and of critical 
crack lengths vs. mean separational stress for a few represent- 
ative materials. Critical crack lengths can be either meas- 
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Fig. 1 


Critical crack length as a function 


ured experimentally or derived theoretically according to the 
formula 


cr 

where 

lee = critical crack length 

F., = critical tear resistance 

E = modulus of elasticity 

o, = P/A = mean tensile or separating stress, in the 


material 


From the curves, critical working stress levels can be de- 
termined. If the actual working stress is below the critical 
level, a crack of length equal to or shorter than the corre- 
sponding critical length will not propagate in fast fracture. 
If the following conditions are satisfied, the possibility of 
fast fracture of the structural shel] will be completely avoided, 
regardless of the design pressure differential: 

1 The critical crack length shall be greater than the 
maximum boundary dimension of the stressed panel. 

2 The maximum working stress in the panel shall be below 
the critical stress level. 

3 The boundary conditions of the panel shall prevent 
propagation of fatigue cracks across the boundary. 

When discussing the likelihood of rapid decompression of 
a sealed space cabin, one is likely to think in terms of the reli- 
ability of the pressurized compartments on current military 
aircraft. The structural concepts for the two applications 
are quite different, as are their potential reliabilities. Aero- 
dynamic considerations for space vehicles will permit the 
sealed cabin to be of a shape more suitable for pressure vessels 
than are the flat-sided, minimum area cross sections of today’s 
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military cockpits. Of particular importance in this regard 
will be the likely disappearance from the sealed cabin of the 
vast areas of brittle, transparent plastic material currently 
required for direct vision. They will undoubtedly be replaced 
by more modest structural glass openings whose reliability 
will be at least as great as the metal structure around them. 

Zaller (6) indicates that the measure of fatigue damage in 
a pressurized commercial aircraft structure is primarily re- 
lated to the number of takeoffs, flights and landings; i.e., 
the number of cabin pressurizations, rather than hours or 
years of service, with some adjustment for various oper- 
ational conditions, such as air turbulence experienced, pilot 
skill, structural maintenance, corrosiveness of the environ- 
ment, towing and ground handling techniques, etc. It 
appears evident that anticipated endurance of a sealed space 
cabin will be far longer than that of the relatively short-haul 
commercial aircraft of today. 

Reliable sealed cabin structure can be constructed from 
known, available materials. In the developing state of the 
art, however, it is anticipated that new materials and tech- 
niques will provide structures having greatly improved tear 
resistance properties. These properties are dependent upon 
a number of factors, including: 

1 Ductility of the material. 

The nature of the crack or penetration. 

The direction of crack relative to grain direction. 
The temperature of the material. 

The dimensions of the stressed skin panel. 

The location of the crack within the panel. 

The effects of restraining or stiffening members. 

For example, sandwich-type construction and techniques, 
such as coining (6), or the pre-stressing of materials in com- 
pression, have been successfully used to provide increased 
fatigue resistance. Of particular interest are the suggestions 
made by Hoffman (7) and Topping (8) regarding the use 
of high strength, nonhomogeneous materials. Properly com- 
pounded, such materials would be extremely resistant to 
fatigue, heat and impact. 


SID Ot GO bo 


Cabin Leakage 


Even though the propagation of fatigue cracks can be 
controlled and critical crack length never attained, fatigue 
crack formation probably can never be completely elimi- 
nated. Leakage rates from such cracks and minute meteoroid 
penetrations will be quite small and exceedingly difficult to 
detect and locate from the inside of the cabin during flight. 
This suggests the vital necessity of a plastic sealant layer on 
the cabin wall to seal such small leaks which would other- 
wise remain undetected until a substantial loss of cabin gas 
had occurred. Leakage of this nature would be significant 
primarily because of this Joss of gas which must necessarily 
be replenished from reserve supplies of oxygen and nitrogen. 


Protection for the Occupants 


Probability of Meteoroid Penetration 


Except for the yet little-known perils of cosmic radiation, 
the possibility of collision with meteoroids in space probably 
offers the greatest external hazard to the astronaut through 
penetration of the sealed cabin. Only penetrations which 
may be survived need be considered, where the mass of the 
meteoroid is small enough so that the occupants survive both 
the impact and probable explosion and are able to initiate 
emergency procedures and effect repairs to the vehicle. It 
is obviously impossible to assign an upper limit to the size o/ 
a penetration which may be survived; however, it is highly 
improbable that structural failure, involving fast fracture of 
the skin as a result of the penetration, would occur unless 
the diameter of the penetration exceeded the critical crack 
length for the structure. If this crack length is assumed to 
be in the neighborhood of 20 in., the probability of collis‘on 
with a meteoroid of such proportions is fantastically small, as 
are the chances of crew survival in the unlikely event of such 
a collision. 

Based on Whipple’s data (9) published in October 1957, 
Table 1 shows the probability of meteoroid penetration of a 
hypothetical 7-ft diameter sphere. Although the indicated 
frequency of collisions and penetrations may be grossly in 
error, considering the scant information at our disposal, the 
significant reduction in probable frequency of penetrations as 
the wall thickness is increased is of great importance. For 
example, each doubling of skin thickness will increase the 
probable time between penetrations by a factor of approxi- 
mately 6.3. Whipple predicts that the use of his proposed 
bumper wall on the space side of the primary sealed cabin wall 
will reduce the number of penetrations of the primary wall by 
a factor of 10 to 100 times. 

Thus, through the combination of a doubled thickness of 
sealed cabin skin and the utilization of Whipple’s external 
bumper wall, which is 10 per cent of the primary skin in 
thickness, we can safely double the cabin pressure, up to a 
maximum of 14.7 psi, and at the same time increase the prob- 
able time between penetrations of the sealed cabin by a factor 
of from 63 to 630. Long duration flights will require corre- 
spondingly long periods between probable penetrations and 
consequent increases in skin thickness. In such cases, a sea- 
level cabin would become more important and structurally 
more feasible at the same time. 


Decompression Due to Meteoroid Penetration 


The physiological effects of a decompression from a given 
cabin pressure to a complete vacuum are something of an 
unknown quantity. Reliance on the RGE factor as an 
indicator is of httle value, since the RGE approaches in- 
finity at 63,000 ft (47 mm Hg). Armstrong (10) reports 


? Based on data by Dr. F. L. Whipple published in Oct. 1957. 


Table 1 Probable meteoroid penetrations of sealed cabin” 


Penetration Operating 
Pressure of skin pressure ———With buffer wall —— 
differential, thickness, stress level,” Meteoroid ————._ No buffer wall,’ days/penetration 

psi in. psi Mass, gm in. days/penetration Factor = 10 Factor = 100 
5.5 0.016 7220 5.67 X 107 0.011 .15 1.5 15 
7.0 0.020 7350 9.56 X 107 0.013 0.26 2.6 26 
10.9 | 7150 3.48 X 10-6 0.020 0.95 9.5 95 
14.7 7700 6.11 X 10-6 0.024 1.67 16.7 167 
14.7 4830 2.14 X 107° 0.037 5.96 59.6 596 
14.7 2 3420 5.67 X 107% 0.051 16.0 160 1600 
14.7 2415 1.38 X 10-4 0.068 37.6 3760 


0.129 


> Assumed sealed cabin: 7-ft diameter sphere (150 ft? surface area), aluminum alloy skins. 
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two decompression experiments which indicate that the 
immediate and critical symptoms of decompression to low 


where 


eae ambient pressures are those of hypoxia rather than of severe —_e Co = 41.43 (n)!? WT. + 459.72 
bodily damage caused by excessive gas expansion. In the = —p, = flow rate (Ib/see) [2] 
is 4 first of these experiments, Luft and Hornberger explosively — ; 
hich decompressed themselves to 62,300 ft while breathing oxygen, where 7 7 : 
the and remained at that altitude for 10 sec before descending. 
bot No symptoms of boiling of the body fluids were evidenced, ine & Pe = Pov re) ; 
tia : and no harmful effects, other than a period of unconsciousness 
It from hypoxia during the descent, agi observed. In the Decompression times for no augmentation are calculated 
ze 0! second experiment, Adelman and Hitchcock explosively de- by direct substitution in Equation [1]. An iterated solution 
ghly compressed 10 dogs from 750 to 30 mm Hg in 0.08 si Some of the two equations is used for the decompression times with 
e of of the dogs were breathing 90 to 100 per cent oxygen during continuous augmentation throughout the decompression. 
nless the experiment while others — breathing air. All sur- Based on the above, Table 2 and Fig. 2 indicate decom- 
rack vived ee to 1 min at 30 mm Hg. Although they pression data and pressure-time relationships both with and 
d to collapsed in 6 to 8 wets and suffered extreme distension of the without augmentation for various hole sizes in a 150-ft* 
ison abdomen, they regained — within 10 min after cabin at both 14.7 and 7.0 psi initial pressure differentials. 
], as vies to ged at All of the animals appeared essen- Also shown are the quantities of augmenting gas required to 
such - | —— ter bed 30 min. ' maintain both a constant pO. during the decompression and 
Time is therefore critical. In ean decompression of the a minimum decompression time of 600 sec to reach 4.36 ps. 
057. sealed cabin should occur, the shirt-sleeved astronaut must For the problem, the following arbitrary conditions were 
a. have time to safely effect certain emergency procedures: iii 
sted To quickly plug the penetration, or to don a suitable pressure 1 The assumed orifice coefficient (0.6) is believed to be 
y in see and rey eles ape which would permit him we reasonably conservative considering the probable quality of 
the either continue his journey or return to base. If 10 min are the edges of an actual penetration. In the decompression 
nie necessary and sufficient to don the pressure garment, the testa reported by Konecei (12), where conventional sharp- 
mac. astronaut must be protected for that length of time against edged orifice plates were used, an orifice coefficient of 0.43 
the loss of useful meres through hypoxia and against brought the theoretical calculations and the measured data 
“a, the disabling effects of decompression sickness or severe into close agreement. 
ind bodily distension from expansion of gas in the gastro-intestinal 2 A constant pO; = 160 mm Hg for both initial cabin 
vall tract. , ; pressures: however, any pOz less than 160 mm Hg can be pro- 
| by The cabin pressure at the end of 10 min should probably vided during the decompression by substituting appropriate 
; not be less than 4.36 pei (30,000 ft), the critical altitude for quantities of nitrogen for equal amounts of oxygen in the 
5 of decompression sickness. To avoid hypoxia, the pO: avail- augmenting gas. 
far able to the man should not fall below 110 mm Hg (11,12) 3 A constant relationship between the weight of nitrogen 
" and should preferably be at 160 to 180 mm Hg to maintain added (WNA) and the total gas flow during any given time 
oa approximately 95 per cent blood saturation. If either of increment. Through a series of preliminary IBM runs, the 
C these minima is exceeded, some augmentation of the cabin relationship dWNA = 0.266 dM, was found to cause the 
‘tor atmosphere during the decompression would be required. pressure-time curves to converge at or very close to 4.36 psi. 
rre- : 1 he pressure-time relationship during a cabin decompres- Except when arbitrarily fixed at a minimum 600 sec, the 
and tee WF the event of . meteoroid or netration can be com- times of decompression are shown to be essentially propor- 
ia puted using the following basic equation. tional to the initial pressure differential and to the cabin 
ally ( 1 ( 2n- volume, and inversely proportional to the area of the pene- 
=< i die tration; similarly, the total amount of oxygen added to main- 
P. = Po \[ @AC 2 \*t!_ /n-1 = 
— (2, tain a constant pO» during the decompression is proportional 
Ve ee n } to the cabin pressure differential and cabin volume, although 
compartment pressure (PSF) [1] the amount remains constant for the various hole sizes. 
ven 
an 
= Table 2 Decompression times with and without augmentation, compartment volume = 150 ft® ’ 
rts Initial compartment pressure 
14.7 psi — 7.0 psi. 
Hole, WGA WOA WN. WGA WOA WNA 
a in. A/V Ratio T p, sec T pa, sec (Ib) (Ib) = (Ib) Tp, sec T pa, sec (Ib) (Ib) (Ib) 
0.02 1.45 10-8 1.16 105 2.18 10° 6.87 6.87 4.27 X 1.10 10° 3.11 3.11 
0.04 5.81 10-8 2.88 X 101 5.44 X 104 6.87 6.87 1.06 X 10 2.74 3.11 3.1) 
0.06 1.31 X 1077 1.28 X 104 2.42 X 10‘ 6.87 6.87 4.74 10 1.22% 10¢ 
0.10 3.63 1077 4630 8730 6.87 6.87 10" 4:38 10" 
0.50 9.06 x 1076 187 6007 iy 11 6 68.4 600° 14 10 4 
| 1.0 3.63 X 10-6 46.3 600 61 43 18 17 600 55 40 15 
2.0 1.45 X 10-* 11.6 =. 237 172 65 4.27 600 216 158 58 
4.0 5.81 X 10-4 2.88 600 — 940 680 260 1.06 600° 859 630 229 
T'p = time of decompression to 4.36 psi without augmentation - aa 
Tpa = time of decompression to 4.36 psi with augmentation except as noted 7 *, r,s a 
WGA = total weight of augmenting gas added during 7T'p4 (WGA = WOA + WNA) eared - SF 
WOA = weight of O, added to maintain pO, at 160 mm Hg aloe dD ot pi 
WNA = weight of Nz added to maintain 7p, at 600 sec minimum (dWNA = 0.266 dMc) 7 
Orifice coefficient 0.6 
Initial cabin temperature 70 F 
* 4.97 psi at 600 sec 
> 4.52 psi at 600 sec : 
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The augmentation with nitrogen and oxygen to provide a 
600-sec decompression time to 4.36 psi appears quite feasible 
for penetrations up to 1 in. in diameter or slightly larger, but 
somewhat less practical for holes approaching 2-in. diameter 
and larger because of the rather considerable weight of aug- 
menting gas required. A 1l-in. diameter hole is not very 
large by many standards; but, if the sealed cabin is out of range 
of all turbines and other rotating machinery, the most likely 
cause for such a penetration is a meteoroid. According to 
Whipple (9), a meteoroid of the fifth visual magnitude is of 
the same proportion (0.834-in. diameter) as the 1-in. pene- 
tration, and probably would strike a 7-ft diameter sphere 
(154 ft?) only once in 91,000 days (250 years). Thus it 
would appear that the suggested augmentation offers 
adequate safety to the astronaut. 


Fire Hazard 


Sealed cabin atmospheres in which ground level partial 
pressure oxygen is maintained at cabin equivalent altitudes 
up to 25,000 to 27,000 ft have been proposed as being satis- 
factory from a physiological standpoint. From an engineer- 
ing point of view, however, this condition must be viewed as 
a serious potential fire hazard. Although little quantitative 
information is available on the effects of a low pressure at- 
mosphere of high oxygen content on the likelihood of ignition 
of various materials or on their burning rate, Clamann, and 
more recently Simons and Archibald (11), have reported 
limited experimental data on burning rates. The latter 
indicate “that anything that will burn at ground level will 
burn much more rapidly at high altitudes at the same pOo.”’ 

Fire has long been considered as one of the greatest of in- 
flight hazards, and both the armed services and the aircraft 
industry have worked ceaselessly toward its elimination. 
Any factor, such as a cabin atmosphere of low pressure and 
high oxygen content, which tends to increase either the 
hazard or the consequences of in-flight fire during normal 
operations must be viewed with grave concern. 


DECOMPRESSION TIME (SECS.) 
Fig. 2 Sealed cabin decompression oe 


Internal Systems 


A sealed cabin system of the highest possible overall effi- 
ciency will be required if very long duration flights in space 
are to be successfully undertaken. Unlike conventional pres- 
surized cabins where system efficiency is generally measured 
in terms of maximum productivity, the sealed cabin and its 
occupants will constitute a rather delicately balanced closed 
system in which efficiency will be measured in terms of ‘‘just- 
enough” productivity, i.e., the degree of equilibrium estab- 
lished and maintained withia the system. 

Once equilibrium is established within the cabin, it is quite 
probable that the weight, the compatibility and the reli- 
ability of the internal cabin systems are essentially inde- 
pendent of the cabin pressure in which they must operate; 
provided, however, that the pressure remains constant. This 
will probably not be the case for a typical sealed cabin unless 
the cabin structure is designed for a 14.7-psi pressure differ- 
ential and ground level pressure is maintained throughout 
the flight. 

If a lower pressure differential is required, the internal 
pressure must change rapidly during the exit and re-entry 
phases so that neither the positive nor the negative pressure 
allowances of the structure are exceeded. Thus, the func- 
tional balance between the cabin systems and the human re- 
quirements which existed at ground level must also be main- 
tained at reduced pressure despite the fact that the physical 
properties of the cabin atmosphere which are dependent on 
the atmospheric pressure will have changed, e.g., air density, 
the convective heat transfer, the boiling point of liquids, etc., 
while other factors, including the human requirements, will 
remain essentially constant. 

Most of the above variables will affect the thermodynamic 
balance within the sealed cabin and, in particular, the systems 
for control of pressure, temperature and humidity. These and 
other systems could be seriously unbalanced as a result of 
rapidly changing conditions within the sealed cabin, and a 
programmed change to a new level of control and operation 
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under each new condition may be required for each system 
affected. 

The requirement that interna] equilibrium be established 
at more than one cabin pressure during any one flight will 
certainly introduce complexity into the sealed cabin system 
with the consequent possibility of a reduction in reliability 
an an increase in weight. 


Weight Considerations 


The selection of an efficient pressure vessel shape for the 
sealed cabin and the use of highly tear resistant materials 
in the cabin structure are of prime importance in providing a 
minimum weight cabin. Although the structural weight 
of 1 static pressure vessel can generally be considered to be 
roughly proportional to the pressure differential, this is not 
the case with a pressurized aircraft or spacecraft structure, 
where all types and combinations of loads are encountered 
during takeoff, flight and landing, which do not necessarily 
add directly to the pressure load or to each other. 

\ structure properly designed to handle all these various 
louds can in many cases permit the application of significantly 
higher internal pressures without a corresponding increase 
in weight. For example, conservative design practice will 
dictate that a manned space structure be designed to the full 
dynamie pressure load to be experienced during exit and re- 
entry without benefit of any compensating internal pressure. 
While in space, a spherical shell designed to withstand this 
external crushing load will accommodate an internal cabin 
pressure at least several times greater than the above dynamic 
pressure at the same margin of safety and with no increase 
in weight. Similarly, it appears that thermal and radiation 
considerations, in addition to high loads from acceleration, 
braking and landing impact, may well require reinforcement 
of the capsule structure and thus permit a differential as 
high as 14.7 psi to be carried without significant weight 
increase, if any. Any weight increase from such cause must 
be balanced against the advantages gained from a sea-level 
cabin system. 


Summary 


In debating the optimum internal atmosphere for a sealed 
space cabin, first consideration must be directed toward 
the overall reliability of the manned sealed cabin system. 
It is concluded that a sea-level atmosphere witbin a sealed 
cabin designed for a 14.7-psi pressure differential offers the 
greatest promise for maximum reliability of the complete 
system. Considerations leading to this conclusion are as 
follows: 

1 The sea-level environment will satisfy all of man’s 
physiological requirements and will ensure an optimum range 
of human performance and endurance. 

2 At the present state of the art, the reliability of pressur- 
ized structure, in terms of its resistance to fast fracture as a 
cause of rapid decompression, is independent of the design 
pressure differential. Although limited fatigue cracking of 
the structural shell probably cannot be entirely eliminated, a 
layer of plastic sealant will prevent the loss of gas through 
such cracks and other minute penetrations. 

3 A thicker cabin shell which will accommodate an in- 
creased pressure differential will also provide a significant 
increase in the probable time between meteoroid penetra- 
tions. Should a meteoroid penetration occur, the higher 
cabin pressure will provide a proportionally longer decom- 
pression time. In addition, augmentation of the cabin 
atmosphere with nitrogen and/or oxygen can protect the 
shirt-sleeved astronaut from the effects of hypoxia and de- 
compression sickness for at least 10 min while he initiates 
emergency procedures and, if necessary, dons a suitable 
pressure garment. 

4 The sea-level atmosphere, as compared with a low 
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pressure atmosphere of high oxygen content, will involve no 
increase in either the hazard or the consequences of in-flight 
fire. 

5 A constant cabin pressure (14.7 psi) would simplify 
the contro] of those internal systems whose performance will 
vary with a change in cabin pressure. The equilibrium 
achieved in the cabin at ground level could thus be maintained 
throughout the flight without the introduction of complexity 
in the form of programmed changes required in system oper- 
ation to maintain cabin equilibrium under a rapidly changing 
cabin pressure. 

6 A typical spacecraft structure will be subjected to a 
multiplicity of severe loading conditions during exit, flight, 
re-entry and, perhaps, impact landing. Considering the 
potential structural reinforcement which may be required 
to withstand such loadings, it is entirely possible that a sea- 
level cabin pressure could be allowed without significant in- 
crease in weight, if any. Should a weight increase occur 
from such a pressure requirement, the increase must be 
offset against the advantages of a sea-level cabin. A weight 
increase in a system which is weight-critical can be and has 
been justified on the basis that the net result improves, rather 
than penalizes, the system. 


Nomenclature 


n = polytropic exponent, assume n = 1.4 for adiabatic ex- 
pansion, dimensionless 

Po = initial pressure in cabin for any dt, lb/ft? abs 

P, = instantaneous pressure in cabin, lb/ft? abs 

6 = orifice coefficient, dimensionless 

A = orifice area, ft? — 

V. = cabin volume, ft* 

To) = initial cabin temperature for any dt, F 

dt = time increment between successive calculations, sec 

Peo = specific weight of gas in cabin for any dt, lb/ft? 

Cy) = speed of sound in air for any dt, fps 

M. = weight of compartment atmosphere, lb e 
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Technical Notes 


_ Preliminary Evaluation of Helium as 
Electric Are Heated Propellant’ 


'EXHE DEVELOPMENT and calibration of an electric arc 

wind tunnel which simulates the gas enthalpies encountered 
by a vehicle re-entering the atmosphere are described else- 
where in this issue.* A preliminary study of the properties of 
helium as an arc heated propellant has been carried out using 
this tunnel. No discussion of the arc heater as a propulsion 
device can really be meaningful without considering a specific 
vehicle and its particular mission. However, in general, the 
gas heater (or accelerator) must produce the value of specific 
impulse which minimizes the total vehicle weight. Moreover, 
the conversion of the electrical power input to the propulsion 
device into exhaust kinetic energy should be accomplished 
efficiently. 

The conversion efficiency limits the maximum specific im- 
pulse which is practical for an electric arc gas heater. The 
reason for this limit is as follows: The voltage across an arc 
is divided into two parts—the drop in the are proper (the 
positive column drop), and the electrode drop. Only the 
positive column drop contributes to useful gas heating; the 
electrode drop is wasted. As the specific impulse is increased, 
the are operates at successively higher temperatures, increas- 
ing the electrical conductivity in the positive column. Conse- 
quently, the electrode voltage becomes a larger fraction of the 
total are potential, reducing the efficiency. 

The are delivers a certain enthalpy to the propellant gas, 
and this enthalpy must be converted to kinetic energy in a 
supersonic nozzle. In the case of a slightly ionized mona- 
tomic gas, the expansion poses no problems; all of the en- 
thalpy can be converted to kinetic energy in a nozzle of 
reasonable area ratio. For higher specific impulses where an 
appreciable fraction of the total enthalpy is invested in ioniza- 
tion, there may be a considerable loss in conversion efficiency 
due to the fact that the ions and electrons will not recombine 
during the expansion. If a polyatomic gas is used as the pro- 
pellant, it will probably be completely dissociated in the arc. 
In the case of hydrogen heated to 10,000 K, about half the 
total enthalpy is due to dissociation. Due to the slow recom- 
bination rate, it is unlikely that this dissociation energy could 
be recovered as exhaust kinetic energy for pressures less than 
about 100 atm; that is to say, the propulsive device can have 
an overall conversion efficiency no better than 50 per cent. 

Although the very low heat of vaporization of helium poses 
difficult storage problems, it does have features which make 
it attractive as an are heated propellant. It is light and 
monatomic. Among the elements, it has the highest ioniza- 
tion potential and can be heated to a high temperature before 
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excessive electrical conductivity reduces the arc efficiency. 
It has a low electrode drop and is nonreactive, simplifying the 
design of durable arc electrodes. 

As previously mentioned, the helium tests were conducted in 
an arc tunnel.* To reduce heat losses, the plenum chamber 
volume was reduced to 3 in.* (one third its normal value). 
Some changes in the gas injectors were necessary in order to 
obtain the desired helium mass flow rates. The supersonic 
(second) nozzle has an area ratio of 10.9; the half angle at 
the nozzle exit is 8 deg. For a monatomic gas, the Mach 
number WV at the nozzle exit is 5.2 ,;where 92 per cent of the 
stagnation enthalpy is converted to directed exhaust kinetic 
energy. This nozzle is nearly equivalent to expansion into : 
vacuum. 

The measurements include gas mass flow m, are power input 
Wo, arc chamber pressure po, and the pitot pressure (tot:! 
head) at the nozzle exit p;. After the first few runs, it was 
found that the ratio ps/po remained steady at a value of 
0.133 + 0.005 and the pitot tube measurements were omitte« 
in later runs. The value 0.133 agrees closely with the ex- 
pected value for the expansion of a monatomic gas in a nozzle 
with an area ratio of 10.9. The mass flow was varied between 
1.07 and 4.70 gm;sec, while the power input was varied be- 
tween 40 and 210 kw. 

The thrust F is given by 


F= 


where 


u 
A nozzle exit area 

~Pr static pressure at the nozzle exit 

Pe ambient pressure into which the nozzle discharges 


For M = 5.2, the second term in Equation [1] can be neg- 
lected. The mass flow is ¥ 


m = puA 
where p is the gas density at the nozzle exit. Therefore 
F = 
But, again for M > 1, ps = pu? sothat | 
F = Ap, = 0.133A po a 


[3] 


The specific impulse J 5p is 
Isp = (0.1334 /m)po 
and the directed exhaust kinetic energy W; is approximately 
W, = F2/2m 


The overall conversion efficiency of electrical power input into 

directed exhaust kinetic energy 7 is then given by 

We 2mW, | 

The test results are summarized in Fig. 1, where the overall 
conversion efficiency is plotted vs. specific impulse. The 
lowest measured value of J;, was 400, and the corresponding 
conversion efficiency was 73 per cent. As expected, the ef- 
ficiency tends to decrease as J,, is increased, but at the 
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maximum measured value of Js, 1120, the efficiency has de- 
creased only to 52 per cent. 


Conclusions 


Klectrode and recombination losses limit the maximum 
specific impulses which can be produced efficiently by an are 
driven propulsion device. Nevertheless, the results presented 
here show that the helium arc can be operated so as to pro- 
duce interesting values of specific impulse at reasonable con- 
version efficiencies. 


Ablation Characteristics of a Subliming 
_ Material Using Arc Heated Air' 
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A procedure is described indicating the use of arc heated 
air in the determination of a material’s ablation charac- 
teristics. Experimental results are reported for Teflon. 
The results are interpreted in terms of a simplified ablation 
theory. 


bbe RECENT development of the electric are air heater 
(1)* has provided a new tool for the study of the behavior 
of ablating systems in dissociated air. The are heater is 
capable of continuously heating air to the temperature and 
enthalpy levels previously attained only on an intermittent 
basis in the shock tube. It is the purpose of this note to report 
on the use of are heated air in the determination of material 
steady-state ablation characteristics. 


Experimental Procedure and Equipment 


In analyzing the steady-state ablation characteristics of a 
wide number of materials, it has been found convenient to 
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= air flow, Ib/sec 


where 
H = total air enthalpy, Btu/lb 
+= are efficiency 
P = power toare, Btu/sec 


define a steady-state thermochemical heat of ablation 


(H Hy) 1 rad 
-& (1) 
H pV pl 
where 
q* = steady-state thermochemical heat of ablation, 
Btu/lb {> 
; de = cold wall heat input, Btu/ft? sec 
H = total free stream gas enthalpy, Btu/Ib 7 
: Hy = total wall gas enthalpy, Btu/Ib : 
(A — H,,)/H = hot wall correction 
= material density, lb/it® 
V steady-state ablation velocity, fps 
rad = net radiation from ablating surface, 


Btu/ft? sec 


The basic experimental setup for the determination of the 
thermochemical heat of ablation is shown in Fig. 1. The 
Aveo electric are air heater has been described previously (1). 
Flat-faced samples of the material to be tested are placed in 
the subsonic exit jet of the arc. By suitable empirical choice 
of exit jet and sample diameter, it is possible to obtain a 
match such that the sample retains its initial flat shape during 
ablation. The steady-state ablation velocity (V) is obtained 
from 35-mm Arriflex film records (24 frames/sec) of the ablat- 
ing sample surface. The net radiation from the ablating 
sample surface is measured with an Eppley thermopile 
and mirror system. The sample surface temperature, and 
thus the hot wall correction (H — H,)/H, is estimated by 
comparison of the spectral radiation characteristic of the 
ablating surface with that characteristic of a black body. The 
spectral radiation distribution is obtained with a Perkin-Elmer 
fast recording infrared spectrometer. 

For each run, measurements are made of power input to the 
arc, power loss to the cooling water, air flow and are chamber 
pressure. The mean air enthalpy of the exit jet is given by 
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thermochemical heat of ablation 
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Test conditions 


Table 1 


50-150 Btu/sec 
0.005-0.025 Ib/sec 
2500-10,000 Btu/Ib 
500-1500 Btu/ft? sec 
1.00-1.75 atm 
1.00 atm 


Power to air, 7P 

Air flow, M.... 

Enthalpy, H........ 

Cold wall heat flux, ge... 

Are chamber pressure, P-. . 
Exit pressure, Po............. 


Fig. 2 35-mm film sequence—ablating Tefion 

T T T 
HEAT FLUX,q¢ = 1500 BTU/ FT°SEC 

AIR ENTHALPY,H = 6600 BTU/LB 

040} ABLATION VELOCITY, V = 0035 IN/SEC 

DENSITY, P = 137.3 LB/FT® 


F ABLATION, q* = 3630 BTU/LB 


T 


HEAT O 


0 20 40 60 8.0 
TIME (SECONDS) 


Fig. 3 Distance-time history of sample surface 
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Fig.4 Thermochemical heat of ablation (q*) vs. air enthalpy (H) 
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On the basis of laminar stagnation point heat transfer 
theory (1,2), the cold wall heat flux to a laminar stagnation 
region from a subsonic jet of dissociated air can be approxi- 
mated by an expression of the form 


KnP (P.\'7 

where 

K = heat transfer correlation factor, (lb/ft sec)'/? 

D = jet diameter, ft 

P, = are chamber pressure, atm 

Py) = exit pressure, atm 


For the determination of the heat transfer correlation factor A” 
the cold wall heat flux was measured with a flat-faced shielde:| 
water calorimeter (1). The correlation factor K exhibits « 
slight dependence on the exit jet enthalpy, pressure and Mach 


number. Details will be published shortly. om | 


A number of 2-in. diameter flat-faced samples of Teflon wer: 
subjected to a 3-in. diameter jet of arc heated air. The 
samples were located 1 in. from the exit nozzle of the arc 
heater. Experiments were carried out over the range of test 
conditions described in Table 1. 

During the course of a typical run (10 to 15 sec) the samples 
were observed to lose between } and 3 in. of their length due 
to ablation. By placing the water calorimeter shield combina- 
tion at different positions within this range along the arc axis, 
it was determined that there was no appreciable falloff in heat 
flux during an ablation test. 

A typical 35-mm film sequence for ablating Teflon is 
shown in Fig. 2. Of particular interest is the manner in which 
the sample retains its initial flat shape. A distance-time his- 
tory of an ablating Teflon surface is presented in Fig. 3. It 
is noted that the sample reaches its steady-state ablation 
velocity almost instantaneously. 

Referring to Equation [1], in the case of Teflon, the hot 
wall correction (H — H,)/H and the net surface radiation 
(qraa) are negligible. Therefore, the thermochemical heat of 
ablation for Teflon is given by 


Experimental Results 


q* = 4-/pV if] 


This, of course, is not true for ablating materials with charac- 
teristic high surface temperatures. The experimentally de- 
termined thermochemical heat of ablation for Teflon as a 
function of air enthalpy is shown in Fig. 4. It is noted that 
the heat of ablation is, within experimental scatter, linear with 
enthalpy. 


Discussion and Summary 


The heat input to an ablating sample is blocked and/or 
absorbed by four mechanisms: Heat conduction to the sample 
interior; surface sublimation, vaporization or reaction; 
transpiration due to the gases evolved at the sample surface; 
and surface radiation. An energy balance at the ablating sur- 
_ face takes the form 


dinput = cond Qvap + trans + 


_ The individual terms are discussed below. : 
The heat input to the ablating sample will be defined as the 
aerodynamic heat input to a nonablating surface which has 
the same temperature as the ablating surface. It is assumed 
that the radiation from the gases incident on the body is 
negligible (3) 


Ginput = (H — Hw)q-/H [6] 

The heat conduction to the sample interior is given by 
qeona = VpCATx — = VeC,AT [7] 
ARS JourNAL 
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where 


C, = material heat capacity, Btu/lb R _ 
T, = temperature of ablating surface, R eee 
[_.. = temperature of unheated sample interior, R ae 


It is assumed that the sample has reached a condition of 
steady-state ablation and has an effective semi-infinite thick- 


ness. 
The energy absorbed by surface vaporization is given by 
Qvap = fpVH, 

where 


fraction of ablated material which is vaporized 
effective heat of vaporization, Btu/lb 


H, 


As used here, surface vaporization is defined as any process 
by which energy is absorbed or released due to phase or 
chemical changes at the sample surface. 

‘The movement of the relatively cool gas from the ablating 
sample surface blocks heat transfer to the surface by both 
thickening the boundary layer and distorting the boundary 
layer temperature distribution (4,5). As a first approxima- 
tion for either laminar or turbulent flow, the heat blocked due 
to transpiration is directly proportional to the mass flux intro- 
duced into the boundary layer at the ablating surface 


Qtrans = a fpV(H Hy) [9] 


where @ is the transpiration factor. > 

The transpiration factor tends to be greater for laminar than 
turbulent flow (6) and decreases as the molecular weight of 
the injected gases is increased (7, 8). 

Combining Equations [5 through 9] the energy balance at 
the ablating surface can be expressed as 


pV|C,AT + fa(H — Hy) + fHy] = 
— Hw)/H) qe — qraa [10] 


Comparing Equations [1 and 10], the thermochemical heat of 


ablation is given by 
q 


* = C,AT + fH. + fa(H — Hw) [11] 


In the case of Teflon, visual observations suggest that all 
the ablated material vaporizes, i.e., f = 1. Thus, the experi- 
mentally observed linear relation between the thermochemical 
heat of ablation q* and air enthalpy H shown in Fig. 4 can be 
predicted on the basis of elementary theory. The slope of the 
curve is a direct measure of the laminar transpiration factor a. 
Similarly, the value of the thermochemical heat of ablation 
at the wall enthalpy H, is a measure of the material heat 
capacity and heat of vaporization, i... CpAT + H,. For 
Teflon, the temperature of the subliming surface is estimated 
to be less than 1000 F; thus, H, ~ 0. Referring to Fig. 4 
for Teflon, the sum of the material heat capacity and heat of 
vaporization is about 600 Btu/lb; the laminar stagnation 
point transpiration factor is approximately 0.48. 

The preceding experimental results obtained at 1 atm are 
in substantial agreement with results for Teflon obtained in a 
supersonic plasma jet at a stagnation pressure of ;/5 atm (8). 
Thus, it is probable that pressure does not have a significant 
effect on the steady-state laminar ablation behavior of 
Teflon. 

The experimental results reported above are preliminary. 
However, they do permit some confidence in the present 
understanding of the ablation characteristics of low tem- 
perature subliming materials. In particular, the main process 
of energy absorption for Teflon is that of transpiration. In 
the case of turbulent flow based on information presently 
available in the literature (6), it is expected that the value of 
the transpiration factor, and thus the thermochemical heat of 
ablation, will be significantly reduced. 
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val Method of Stabilizing an 


Astronomieal Satellite! 


ROBERT R. NEWTON? 


The Johns Hopkins Univers Silver Spring, Md. 


N THE operation of a satellite being used for astonomical 
observations, it is clearly advantageous for the satellite 
to have zero angular momentum with respect to the fixed 
stars. However, this condition, even if achieved at some 
epoch, cannot be maintained indefinitely; there are always 
small torques arising from meteorite impacts, gravitation 
and radiation pressure, acting to alter the angular momentum. 
There will therefore be a need for eliminating angular momen- 
tum in small amounts over a long period of time. 

Let us take as a requirement that the stabilization system 
used must reduce to zero, within one period, the angular 
momentum corresponding to an angular velocity of 1 deg per 
period. We suggest the use of a system which unbalances, in 
a controlled way, the mass distribution of a satellite and thus 
produces gravitational torques on the satellite. We wish to 
show that such a system is feasible. 

Let Ry be the vector’ from the center of the Earth to the 
center of mass of the satellite, and let go be the acceleration 
due to the gravitation of the Earth at the distance Ro. Take 
the center of mass as the origin of a coordinate system X YZ 
with Z positive upward, and let J J K be unit vectors along 
XYZ, respectively. Let p be the vector from the origin to 
any mass element dm of the satellite. Then the gravitational 
force dF on the mass dm has the magnitude 
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go [Ro?/|Ro + P|2]dm, 
and the direction 


— (Ro + ?)/|Ro + P| 
Therefore, the vector torque on the satellite is 
T=fpx dF 
= X (Ro + p)|/\Ro + 
Since p = 1X + JY + KZand Ry 
the torque readily reduces to 
T = + JI xz) [1] 


through first order in the dimensions of the satellite, for a 
spherical earth. Iyz and Jxz denote products of inertia.‘ 

The angular momentum vector has some direction in 
space, and we shall suppose that there is some method of 
sensing this direction. Since the X-axis is not yet defined, we 
can take it to lie in the vertical plane containing the angular 
momentum vector. Since there is no Y component of angular 
momentum, we want xz = 0. In order to obtain maximum 
torque about the X-axis, we want /yz as large as possible. 
The value of Jxy is immaterial. 

Although we cannot affect the component of angular 
momentum which is vertical at some particular instant, the 
direction of Z is constantly changing, and a component which 
is vertical at some instant will be horizontal one quarter of a 
period later. Thus, all components of angular momentum, 
taken with respect to the fixed stars, can be affected by the 
gravitational torques. 

In order to control the products of inertia Jyz and [yz, 
suppose that the satellite contains six equal masses m/2, 
which can move along three orthogonal axes xyz, with one pair 
of masses along each axis. The motions of a given pair will 
be equal and opposite, so that the center of mass is not affected 
by the motion. (Using only one mass along each axis would 
probably be acceptable; this would introduce a small amount 
of coupling between the various axes.) Further, assume that 
for some sct of positions of the masses, the moments of inertia 
of the satellite are all equal, and the products vanish. The 
directions xyz will then be the principal axes of inertia for any 
other configuration. Let 29 yo 2 be the positions of the masses 
in the balanced configuration. 

The conditions that Jyz = 0 and that Jyz be a maximum 
are sufficient to determine xyz in most cases. Let the direction 
cosines connecting xyz and X YZ be denoted by the scheme 


K Ro, the expression for 


Z 
z ny Ne n3 


The products of inertia Jyz and [yz are (p. 123, (2))5 
Iyz = + I,ymem; + 


Since = Jy, = when x = 4%, ete., and Il; + mam, + 
nny = 0 for i ¥ j, Equations [3] become, using Az? for 
x? — 27, ete. 


m|lel;( Az? — Ax?) + mom;( Az? — Ay?) ] 


‘The potential energy and torque on a nonspherical satellite 
have been derived in several places, for example, in Ref. (1). 
The writer has not seen this form for the torque obtained else- 
where. In deriving this form, the writer originally made the 
mistake of taking the infinitesimal force dF to have the direction 
of Ro, rather than of Ro + »; this error leads to a factor of 2 
rather than 3 in 7. The importance of this error was pointed out 
to the writer by Dr. E. Shotland of this laboratory: Of the total 
torque, two thirds comes from the variation in the size of gravity, 
and one third from the variation in direction. 

5 Numbers in parentheses indicate References at end of paper. 
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Most of the time, none of the direction cosines vanish. Even 
if some of them should vanish at some instant, they will not 
a short time later. Assuming that Js, and msm, do not vanish, 
we can use the condition Jyz = 0 to obtain 


Az? — Az? = — (Az? 
= — lelymgm)( Az? — Ay?) /Isl, 


Ti 3," 


The direction cosines depend upon the orientation of the 
satellite with respect to gravity and to the total angular 
momentum, and are presumed known. Suppose J3lm2m; — 
llymym;> 0 and I3l,> 0. Then, in order to make [yz a maxi- 
mum in the second of Equations [5], we try to make Az?, and 
hence z, positive and as large as possible, and y as small :is 
possible. If this cannot be done within the limits of tlic 
geometry without violating the first of Equations [5], we 
make Az? — Az? as large as possible with the correct sign, 
and give the corresponding value to Az? — Ay?. 

These conditions lead to a “bang-bang” type of operation, 
with two of the pairs of masses usually at their extreme 
positions. Since an interval of several minutes is a short tinie 
for this application, this does not constitute a severe require- 
ment; however, other modes of operation, such as making J, z 
proportional to the angular momentum, might be preferable. 

To estimate the values of m and Az, etc., needed, let \f 
be the mass of the satellite and & its radius of gyration, and 
suppose that the average value of Jyz obtainable for a full 
period is the maximum of }m(Az? — Ay?). For the torque 
requirement stated in the second paragraph, ignoring the 
sign 


m( Az? — Ay?) = goT? 


where 7’ is the satellite period. For a circular orbit, go = 


giving 
m( Az? — Ay?) = Mk?/10807r 


Suppose further that 2 = yo = k, and that each mass can be 
diplaced a distance D from the balance position, and use 
Az? = —Ay? = (k + D)? — k? = QkD. Then 


mD = Mk/4320r [6] 
It is reasonable to allow D = k/20, requiring a total movable 
mass 


3m 


M/724r 


or a requirement of less than 3 per cent of the total mass for 
this control function. 
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Adiabatic Wall Temperature Due to 
Mass Transfer Cooling With a 
Combustible Gas 

at 

D. B. SPALDING! 
Imperial College, London, England 


RECENT technical note by Sutton (1),? with the 

above title, discusses the influence of the burning of a 

transpiration coolant on the quantity of coolant necessary to 
Received June 17, 1959. 
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maintain a given wall temperature. The present note dis- 
cusses the same problem in a way which has been found useful 
in calculating the burning rates of solid and liquid fuels (2). 

Consider the transpiration cooling of a porous surface in a 
gas stream. Then a simple modification of the general mass 
transfer equations of (3) shows that, for unity Lewis number, 
we have to solve the usual equation 


GVb — V(Dpvb) = 0 [1] 
and boundary condition Pu, ; 


= |Dp(ob/dy) |w [2] 

where the conserved property 6 is defined by — 
b =h/(hy — hi) [3] 

Here h is the enthalpy of the gas (chemical terms being in- 
cluded), h,, is that at the wall, and h; is that of the coolant in 
the reservoir from which it is drawn. : 
This problem has solution 


where 
B= (he — hu)/(hw — hi) [5 
and the function F( ) depends on the geometry of the sys- 
tem. Once obtained, this solution permits calculation of the 
coolant rate m” for a value of B which is specified by stream, 
wall and reservoir conditions. m” always increases with B, 
but not proportionally. 

The present contribution concerns the determination of B 
by means of an enthalpy-composition diagram, without the 
assumption of constant specific heats and taking full account 
of chemical reaction. Fig. 1 shows an enthalpy-composition 
(h-f) diagram for gaseous hydrogen-air mixtures, based on 
data from (4). The kinks appearing in the isotherms are at 
the stoichiometric fuel fraction fstoicn = 0.0284. 

Suppose, as an example, that the gas stream has a stagna- 
tion enthalpy of 8000 Btu/lbm, corresponding to a relative 
velocity of approximately 20,000 fps; the appropriate state- 
point on the A-f chart is marked G. Suppose that the hydro- 
gen in the coolant reservoir has temperature 500 F abs; the 
appropriate state-point on the chart is marked F. Suppose 
that the surface has to be cooled to 2000 F abs; then the 
corresponding state-point S lies somewhere on the 2000 F abs 
isotherm. Now it is a consequence of the theory of mass 
transfer (see, for example, (5)) that, for unity Lewis number, 
the state-point S lies on the line joining G and F; it therefore 
lies at the intersection of this line with the 2000 F abs iso- 
therm. 

Inspection of Equation [5] shows that we can therefore 
write 


B = GS/SF [6] 
So a measurement of this length ratio on the diagram enables 
the ‘‘transfer number” B to be evaluated with ease; B = 1.5 
in the example shown. 

If combustion did not occur, all the isotherms would be 
“pulled tight” through the same points on the pure-air and 
pure-hydrogen lines. In our example, the surface condition 
would therefore be represented by the new state-point S’; 
the corresponding value of B is 1.28. The coolant rate m” for 
the new condition will be less than that for the old; the reduc- 
tion will be smaller than proportionate, but the solution 
}quation [4] for the geometry in question must be examined 
to find out precisely how big it is. 

The use of the h-f chart for such problems offers the follow- 
ing advantages: 

1 There is no necessity to assume constancy of specific 
heats. 
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2 No knowledge of the stream temperature is needed. 
This is just as well, since its calculation under these highly 
dissociated conditions is very difficult. 

3 In constructing the chart, dissociation of the gases near 
the stoichiometric mixture ratio can be ignored, for this merely 
rounds the sharp kinks in the isotherms in this region, but 
does not alter the positions of the isotherms at the moderately 
large values of f at which S is likely to be located. 

4 Inspection of the graph makes it easy to see how the 
transfer number B will alter in magnitude with changes in the 
specification of the stream and reservoir conditions and of the 
wall temperature. Furthermore, it shows clearly how, with 
specified temperatures, B is bound to increase greatly if a 
coolant is used which has a specific heat (spacing of isotherms) 
nearer to that of air than that of hydrogen. 

5 All these conclusions are independent of the geometry 
of the flow (e.g., stagnation point, flat plate, rocket nozzle, 
laminar or turbulent, etc.). It is only when translating B into 
terms of m” that the geometry must be considered. 

6 An exact analysis of the flow in particular conditions 
requires consideration of property variations in the stream. 
The chart facilitates this, for all fluid states appearing in the 
boundary layer (assuming unity Lewis number and equi- 
librium) lie on the line GS. 

The particular geometry considered by Sutton (1) was the 
axisymmetrical stagnation point with a laminar boundary 
layer. He neglected the pressure gradient and assumed unity 
Prandtl number, using the flat plate solutions of Emmons and 
Leigh (6). The inaccuracies resulting from these assumptions 
may be deduced from Fig. 2, which contains various m” vs. B 
relations in dimensionless form; curve (a) is the one used by 
Sutton. 

Curve (b) also holds for the flat plate; it is obtained from 
(7) which is based on the collection of boundary solutions in 
(8); it holds for Pr = 0.7, which is the value for air. Com- 
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B 
Fig. 2 Transfer (coolant) rate vs. transfer number B for 
various laminar boundary layers 


parison with curve (a) shows that the Pr = 1 assumption 
somewhat underestimates the coolant rate m” as is to be ex- 
pected. 

Curves (d) and (e) hold for the axisymmetrical stagnation 
point; the lower curve has been plotted by interpolation in 
the data of (7) followed by use of the Mangler transformation, 
while the upper represents the older K4rm4n-Pohlhausen 
solution of (9). Both hold for Pr = 0.7, but curve (d) is the 
more accurate. Comparison of curves (b) and (d) shows that 
neglect of the pressure gradient causes the necessary coolant 
rate to be greatly underestimated. Roughly speaking, the 
coolant rate must be three times as great as the values pre- 
dicted in (1). 

Curve (c) holds for the two-dimensional stagnation point. 
It is plotted for interest only. All the solutions hold for con- 
stant properties. In the absence of solutions for varying 
properties, it is probably wisest to use the corrections which 
Fay and Riddell (10) found to be valid for axisymmetrical 
heat transfer in the absence of mass transfer. 

In addition to the quantitative difference between the 
stagnation point and flat plate solutions, there is also a 
qualitative one: Whereas, for the flat plate, the ordinate of 
Fig. 2 approaches the asymptotic value of 0.6193 as B tends 
to infinity, there is no such upper limit when wu, increases with 
x. In the former case, the boundary layer is “blown away” 
at the limiting mass transfer rate; i.e., the velocity gradient 
at the surface tends to zero. In the latter case, however, con- 
sideration of the equation of motion of the boundary layer 
shows that the surface velocity gradient tends to the positive 
value (peue/m”) du,/dx; so the boundary layer thickness re- 
mains finite for all finite values of m”. 


Nomenclature 

b dimensionless conserved property 

B be — by, “transfer number’’ or driving force for mass 
transfer. See Equation [5] 


= typical dimension of body 
diffusion coefficient j 
mass of coolant per unit mass of local mixture a. 


= mass velocity (vectorial) of gas _ 
enthalpy of gas mixture 
mass flow rate of coolant per unit area of surface 
velocity in zx direction 
distance along surface from leading edge or stagnation 
point 
= normal distance from surface toward gas stream 
dynamic viscosity of gas 
kinematic viscosity of gas 
density of gas 


Subscripts 


$ ter, 
gas stream conditions 
= coolant gas reservoir 7 


w surface of porous skin exposed to gas stream 
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Hypersonic Conical Flow 
HARRY WILLIAMS! 


California Institute of Technology, Pasadena, Calif. 


HE FOLLOWING analysis presents an approximate 
solution for hypersonic flow about a semi-infinite cone at 
zero angle of attack. Simple expressions for surface velocity, 
pressure and temperature ratios are derived and compared 


with those calculated by Kopal. 


The flow of an inviscid, nonconducting perfect gas ro a 
semi-infinite cone at zero angle of attack with an attached 
shock wave can be formulated in terms of the equations 


General Theory 


avg = 
VXq= 
y¥-1 
+ Ma? (1 — 


The boundary conditions are specified on the shock wave and 


Received June 17, 1959. 
1 Senior Research Engineer, Prelimary Design and Research, 
Jet Propulsion Laboratory. Member ARS. 


ARS JourNAL 


— 
° 
mx 
Ab 
4 
*6 
Biss 
= 
; 
- 
» 
bet | | 


ion 


on the body surface as 


= Wy) = cos w 


On introducing conical coordinates, Equations [1, 2 and 3] 
with the tangency boundary condition lead to 


1 dua 1[a%n) (= 
Ke) - af, ] u(n)dn [4] 


Approximate Solution 


As (a) has been observed to be small at high Mach num- 
bers, it follows from Equation [4] that from an order of mag- 


nitude standpoint 


O(a) di/dh = 0(a?) 


Thus a(w) is a slowly varying variable and can be approxi- 
mately replaced in the integrand of Equation [4] by its mean 
value (w%m) where (Um) is arbitrarily taken to be 

Um = = w-) + Ww = wx)] 
In addition, since é is small, the speed of sound a(w) will also 
be a slowly varying variable. Hence, the integrand of Equa- 
tion [4] is approximately stationary and leads to the approxi- 
mate solution 


d(x) = —2umx 
W(x) = Uc — Umx? 


On applying the boundary condition, the shock wave angle is 
determined from the solution of 


2/(y — 1) 2a 
ti 


and the cone surface velocity follows as .. 


( 1+ a?/2 

ue =(- cos w 

1 — a2/2 

The parameters of greatest interest are the pressure and 


temperature ratios at the surface of the cone These are ob- 
tained from Equation [3] and the condition of isentropy as 


Pw Tx 8 1)/2y]/M? 


The limiting behavior of the wave angle (w,) and the pressure 
and temperature ratios as M.—> © is clear from these rela- 
tions. 

The limits of applicability of this solution can be approxi- 
mately fixed by substituting this solution into Equation [4] 
and noting the terms that were neglected. In this way, it 
follows that the approximate solution is appropriate, pro- 
vided 


a cot Wy < 1 


Thus, this solution would not apply to very slender cones ex- 
cept in the case of very large Mach numbers. A comparison 
of this solution with that of Kopal is presented in Figs. 1, 2 
angle. 


and 3 for several values of cone 


Nomenclature 
« = speed of sound 
h = (w — we)/a 
\f = Mach number 
M, = Mo sin Wy 
p = static pressure 
¢ = velocity 
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q/U 

temperature 

radial velocity 
u/U 

free stream velocity 
angular velocity component 
v/U 

(w — we) 

dummy variable 

standoff angle 

ratio of specific heats 

(z = we) /a 

polar angle 


= 


Subscripts 
condition on cone surface 
condition immediately behind shock wave 
condition ahead of shock wave 


Thermal Protection of a Re-Entry 
‘Gr ©, Satellite’ 
Na 


SINCLAIRE M. SCALA? 
General Electric Co., Philadelphia, Pa. 


HE DESIGN of a satellite for successful re-entry from 

orbit requires consideration of the heat flux generated by 
the vehicle as it enters the Earth’s atmosphere. This problem 
has been considered by a number of investigators (1 to 5),° and 
the heat transfer can be shown to be a function of such quan- 
tities as the ballistic parameter of the vehicle W/CpA, the 
entry velocity, the entry angle, and aerodynamic lift and drag 
forces, if employed. Once the type of entry vehicle is selected, 
the entry trajectory and the heat flux are also established, and 
then attention must be focused on the choice of heat protec- 
tion scheme to be used. 

Three basic types of heat alleviation schemes are the heat 
sink, the regulated mass transfer system (e.g., transpiration 
cooling) and the self-regulating mass transfer system (e.g., 
ablation cooling). Analyses are therefore presented for these 
mass transfer systems, and for purposes of comparison, a heat 
sink calculation is also performed. 

Results are presented for the ballistic re-entry of a satellite 
for three different values of the ballistic parameter W/Cp4A, 
i.e., 50, 100 and 200. For all three cases (6), re-entry is from 
an initial satellite orbit of 900,000 ft. For comparison’s sake, 
the velocity at initiation of re-entry is approximately 24,000 


fps at a path angle of 92} deg from the local vertical. This 


is achieved by firing retrorockets at time zero. _ 


Heat Sink Considerations 


Upon applying the theory of (7) to the three re-entry tra- 
jectories, one obtains the results shown in Fig. 1. One ob- 
serves that there are three sets of curves representing the 
aerodynamic heat transfer at each value of the ballistic param- 
eter W/CpA. This follows because each of the thermal 
protection systems to be considered here operates at a dif- 
ferent characteristic surface temperature. The effects of gas 
cap radiation are not included in these curves and are con- 
sidered separately. For all heat protection systems, it has 
been assumed that the surface emissivity is €, = 0.1. 
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A one-dimensional transient heat conduction calculation 
was first performed (8), and the typical results obtained are 
shown in Fig. 2, where time zero corresponds to the initiation 
of significant heating. 

In the case of a heat sink, the two major criteria which de- 
termine the thickness of the wall (and hence the mass of solid 
required) are that: (a) the surface temperature should re- 
main below the softening point of the material, and (b) the 
inner wall temperature must remain below the temperature 
which destroys the payload. 

Note that the outer wall temperature reaches a maximuin 
and then drops, while the inner wall temperature rises steadily . 


7 Transpiration Cooling 


In determining the mass of coolant required to maintain 
the surface of the satellite at a constant temperature, onc 
must resort to a mass and energy balance at the surface. On: 
obtains as the required flux of coolant 


(hyp — hres) + (AQ/Ati)w 


_ where the “effectiveness quotient”’ has been obtained (9) by 


the author as 


AQ 30 Btu 
—) [240 0.48 (he —h 


where the enthalpy h is in Btu/Ib. 

Except at very low flight velocities, the major contribution 
to the denominator is the “effectiveness quotient.” By em- 
ploying an oxidation resistant porous wall which can tolerate 
a surface temperature of 3500 R, sufficient heat will pass 
through the wall so as to cause the phase transformation from 
liquid air to gaseous air to occur between the reservoir and th: 
porous wall. Upon assuming a reservoir temperature of 130 R 
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and employing the data of (10 and 11) and the correlation 
Equation [2], one may calculate the instantaneous coolant 
flux required. 

The calculations were performed for a satellite having a 
nose radius of 6 ft and a surface emissivity of €. = 0.10, and 
are shown in Fig. 3. Note that due to high altitude decelera- 
tion, the shock layer radiation was negligible, and hence the 
results are also applicable at the stagnation point of a satellite 
having a different nose radius, provided €y+/R, is maintained 
equal to 0.245 ft”*. It is, of course, clear that the area under 
the curves shown in Fig. 3 represents the total coolant flux 
required during re-entry. 


Vaporization Cooling 

Under conditions of severe heating, a refractory material 
behaves initially as a heat sink, and then because of the ma- 
terial’s inherent low thermal conductivity, the surface tem- 
perature rapidly rises to a level sufficiently high so as to cause 
melting and vaporization to occur at the surface. 

When melting and vaporization occur, the total heat 
blocked or absorbed during the ablation process includes three 
separate contributions. These are the heat blocking action 
due to the thickening of the gaseous boundary produced by 
mass transfer, the heat absorbing action due to phase changes, 
and the heat blocking action due to convection in the liquid 
phase. If the symbol T is used to represent the fractional 
part of the solid which gasifies, then analysis shows that Q* 
the “effective heat of vaporization” may be written (12) 


AQ 
=T [ (42), + (hw — h(s)) + - no) | [3] 


which includes the three contributions discussed. 

If certain minimum physico-chemical data are known con- 
cerning a refractory material, then the dependence of Q* on 
environmental conditions can be determined. Thus, for the 
typical refractory material characterized by the chemical sym- 
bol ZOs, it was found (12) that Q* depends primarily upon the 
stagnation enthalpy and rises almost linearly with this quan- 
tity (or as the square of the flight speed V..”). 

If sublimation of ZO2 occurs (' = 1.0), one finds that Q*, 
which for this quartz-like material includes a large heat of 
vaporization, correlates within 5 per cent error in the form 


(Q*)r=1.0 = 8580 + 5.00X10*V.2 Btu/lb [4] 


where V., is expressed in fps for the entire hypersonic flight 
regime. 

It has also been found that (12), for a given flight speed, 
Q* is practically a linear function of TI’, hence it is convenient 
to write 


(Q*)r 410 = (Q*)rai0 — (1 — T)AQ*/AT [5] 


The first term on the right-hand side is given by Equation [4], 
and the reduction in Q* per unit decrease in T has been cor- 
related by 


* 
= 7430 +5.00x10¢V,2 [6] 


Ib 
where V. is expressed in fps, a form which differs from 
Equation [4] only by a constant (heat of fusion). 

The gasification ratio I is critically dependent upon the pre- 
cise temperature dependence of the vapor pressure of the 
vaporizing species and the viscosity of the molten layer, and 
consequently will vary with environmental conditions, and 
can be obtained by solving and matching the gas phase and 
liquid phase boundary layer equations (13). 

Qualitatively, however, one may anticipate that the actual 
(' of a material which exhibits a liquid phase will drop as the 
stagnation pressure rises and/or as the stagnation enthalpy 
drops. For purposes of discussion, however, it was assumed 
here that T’ is constant along a given trajectory and equal to 
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0.25 for one crystalline form of ZO. and equal to unity for the 
other. 
An energy balance yields as the total instantaneous mass of 
solid required for a quasi-steady ablation process 


Mot = Q* i 


It should be noted that this expression is quite similar to 
Equation [1] derived for transpiration cooling. 

Although it was found that reradiation can be conveniently 
neglected in the case of the heat sink, and its neglect will not 
cause an appreciable error for the transpiration cooling case 
(T,, = 3500 R), one cannot neglect the effects of surface re- 
radiation when considering a material which melts at a very 
high surface temperature (7, = 5000 R). 

Upon utilizing Equations [4 to 7] and an assumed surface 
emissivity parameter of €y\/R, = 0.245 ft”?, the rate of 
mass loss for the three ballistic parameters was determined and 
is shown in Fig. 4 for the two aforementioned values of the 
gasification ratio (I' = 0.25, 1.0). 


Conclusions 

Theoretical analyses have been presented for several dif- 
ferent heat protection schemes, where the effectiveness of 
the various systems is compared in terms of the total mass 
of material required for each at the forward stagnation point. 
The comparison is shown in Fig. 5. This study shows that 
the intrinsic mass of coolant shield (less appurtenances) re- 
quired for the thermal protection of a re-entry satellite is a 
sensitive function of the cooling scheme selected, and, further, 


the heat sink does not appear to be competitive on a weight 
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basis, with the mass transfer cooling systems. (The heat sink 
results do not appear because they lie off the graph.) 

The regulated cooling system which utilizes air is not op- 
timal from the point of view of the choice of the coolant gas 
(a lighter gas, such as helium would undoubtedly result in a 
smaller total mass requirement); however, the relative per- 
formance of air cooling under the conditions of hypersonic 
satellite re-entry is established, and, further, an excellent in- 
dication is given of the enormous blocking action of the mass 
injection process. It is also noted that the value of (AQ/Am),, 
. given in Equation [2] is nearly equal to the value of the 
“effective heat of ablation’ of many plastic materials which 
do not have an appreciable heat of combustion; hence the 
curves marked “‘transpiration cooling” are also indicative of 
the behavior of Teflon-like plastics. 

(It should be noted that only the intrinsic total mass re- 
quired for thermal protection is shown in Fig. 5, and hence in 
the case of transpiration cooling the weight of the storage 
tanks and valve system is not included.) 

The self-regulated vaporization cooling system which 
utilizes a quartz-like ablating material, having T = 0.25, 
represents a 50 per cent increase in efficiency over air cooling. 
Sublimation of a quartz-like material represents a 400 per cent 
increase in effectiveness over air injection (see Fig. 5). Note 
that the class of ablating refractories considered herein is 
attractive, because the operating surface temperature is so 
high that over certain portions of the trajectory the material 
actually provides a radiation cooling system, and the mass 
loss by vaporization is somewhat incidental to the overall 
process. With regard to sublimation, it must be understood, 
however, that it is unlikely that a value of I equal to unity 
will ever be obtained for a noncombusting refractory material, 
and hence the assumption that !' = 1.0 is optimistic for quartz- 
like materials. In the case of vaporization cooling, the as- 
sumption of quasi-steady ablation implies that the thermal 
diffusivity is sufficiently small so that conduction within the 
refractory material proceeds no faster than the ablation rate. 
If the process is truly quasi-steady, no insulating material is 
required behind the vaporizing material. On the other hand, 
if the ablation process is transient due to the long re-entry 
time, additional shield weight will be required in the form of a 
composite wall consisting of the high temperature refractory 
backed up by a low weight insulator. 

Finally, it is clear that experimental verification of the pre- 
dicted interaction between a material and its environment is 
also a necessity. Therefore, related experimental studies are 
currently being performed (14,15). 
Nomenclature - 

area 
static enthalpy of gaseous mixture 
mass transfer rate of mixture 
mean molecular weight of mixture 
heat transfer rate 
= “effective heat of vaporization” 


her ARB. 


“effectiveness quotient’’ of mass transfer process 
reduction in Q* with decrease in gasification 
nose radius of satellite 

solid phase 
temperature 

flight speed of satellite mn 

ballistic parameter 

pure element 

density 

emissivity 

Stefan-Boltzmann constant 

gasification ratio 


Subscripts 


aero aerodynamic 
edge of boundary layer 
gas cap 
liquid phase 
reservoir 
total 
wall, surface of condensed phase 


Acknowledgment 


The author would like to acknowledge the assistance o/ 
Leon Gilbert and Miss Eileen Foxhill in performing the calcu- 
lations and preparing the graphs, and Paul Gordon and J. A 
Royer who programmed the analysis for the IBM 704 digita 
computer. The author also wishes to acknowledge th: 
courtesy of H. Bloom, E. Nolan and C. Allen of the Genera! 
Electric Co., M. 8. V. Department, Philadelphia. 


References 

1 Allen, H. J. and Eggers, A. J., Jr., ‘‘A Study of the Motion and Aero- 
dynamic Heating of Missiles Entering the Earth’s Atmosphere at Hig! 
Supersonic Speeds,’’ NACA TN 4047, Oct. 1957. 

2 Gazley, C., Jr. and Masson, D. J, ‘‘Recovery of a Circumlunar In- 
strumented Carrier,’’ ARS preprint no. 488-57. 

3 Riddell, F. R. and Detra, R. W., ‘‘Returning Alive from Space,”’ 
AAS preprint no. 57-4. 

4 Hoglund, R. and Thale, J., ‘‘Recovery from Satellite Orbit,’’ ARS 
preprint no. 650-58. 

5 Paige, H. W., ‘‘The Technology of Manned Return from Outer 
Space,’’ Redding Lecture, Franklin Institute, Nov. 1958. 

6 Bloom, H., private communication. 

7 Seala, S. M. and Baulknight, C. W., ‘‘Transport and Thermodynamic 
Properties in a Hypersonic Laminar Boundary Layer—Part II Applica 
tions,’’ General Electric Co. TIS Rep. no. R59SD306, Feb. 1959. 

8 Nolan, E. and Allen, C., private communication. 

9 Scala, S. M., ‘‘Transpiration Cooling in the Hypersonic Laminar 
Boundary Layer,’’ to be published. 

10 White, H., ‘‘Modern College Physics,.’’ Van Nostrand Co., N.Y., 1951. 

11 Nat. Bur. Standards Staff, ‘‘Tables of Thermal Properties of Gases,’ 
Circular 564, Nov. 1955. 

12 Seala, S. M., ‘‘Sublimation in a Hypersonic Environment,’’ General 
Electric Co. Rep. no. TIS R58SD289, Oct. 1958. 

13. Seala, S. M. and Sutton, G. W., ‘‘The Two-Phase Hypersonic Lamina: 
Boundary Layer—A Study of Surface Melting,’’ 1958 Heat Transfer and 
Fluid Mechanics Institute, Berkeley, Calif., June 1958. 

14 Gruntfest, I. J. and Shenker, L. H., ‘‘Behavior of Reinforced Plastics 
at Very High Temperatures,’’ Modern Plastics, vol. 35, no. 10, June 1958, 
pp. 155-166. 

15 Dank, M., Nelson, R. A., Sutton, W. H. and Sheridan, W. R.., 
‘‘Water-Stabilized Arc Tests on Non-Metallic Materials,’’ J. Electrochem. 
Soc., vol. 106, no. 4, April 1959, pp. 317-321. 


Changes in the Density and Temperature 
of Liquid Propellants in Rocket Engine 
Supply Systems 


S. L. BRAGG! 
Rolls-Royce Ltd., Derby, England 


The basic thermodynamic relations are used to derive 


equations for the changes in propellant density and tem- 
_ perature in rocket engine supply systems. It is shown that 
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in a pump with pressure rise AP, these changes are given 


by 
Ap AV ) ] 
—-— = — = AP | — +aT) 
Br 
VAP(1— 17 
AT = —(—— +aT 
C, 


where the symbols have the conventional meanings. 
In a pipe in which friction causes a pressure drop — AP 


Ap AV aV 
aT A (1 — aT) 


Pp 


Numerical values for typical propellants show that these 
changes, although small, are not insignificant. ; 


N ANALYZING the performance of a liquid rocket engine, 
t is necessary to know how the density and temperature 
of the propellants vary through the supply system. The 
pump delivery head, for example, must be based on the fluid 
density at outlet, not the value at inlet which is that normally 
measured. Formulas for the estimation of density changes in 
steady flows of liquids do not, however, appear to be given 
in the standard works of reference, either on thermodynamics 
oron pumps. Such formulas are derived in this note, and are 
used to determine the density and temperature changes occur- 
ring in typical pump and supply line systems. 


Basic Equations 


In any adiabatic steady flow process, from the first law of 
thermodynamics 


+ E, + 42/29 + W = {1} 


where 
p = absolute static pressure 
\’ = specific volume, the inverse of the density 
v = velocity 7 
E£ = internal energy 
W = external work done on unit mass of the fluid in the 
process 


(Suffixes 1 and 2 refer to conditions of the fluid at inlet and 
outlet respectively.) 

The velocity terms may be converted to dynamic pressure 
at the appropriate density, and the equation for the total 
pressure P = p + v?/2Vg substituted into Equation [1]. 
This gives the equation for the change in internal energy as 


AE BE, -E, 

= W+P,V; — PV: 

= W — APY; — P,AV 


[2] 


Two basic processes may be considered: (a) steady flow 
pumps, in which W is large and known; (b) lines, injectors, 
etc., in which W = 0. 

In both processes the information from test results, as to 
work done, pressure differences, initial specific volume, etce., 
is sufficient to enable the change in internal energy to be 
calculated to a reasonable accuracy from Equation [2]. 

Since only two properties, such as internal energy and pres- 
sure, are needed to define the state of any stable liquid, the 
change of density in any process can be determined uniquely 
from the changes in internal energy and pressure. That is 


AV = ar + (or) at [3] 
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where the subscript denotes the property that is kept constant 
during partial differentiation. 

We therefore need to evaluate these partial derivatives in 
terms of the normally measured physical properties of the 
fluid. 

Evaluation of Partial Derivatives, (OV/OP), and (OV/OE)> 


In any small change in which heat 6Q is added to a sub- 
stance 


OE = 6Q — P-6V 


Thus at constant volume, from the definition of the specific 
reat 


\ar/y \aP/y \or/y 


Similarly, at constant pressure 
o¢ 
o7 P oT P oT P 
E V 
oT) p 


The physical characteristics of a liquid that are normally 
measured are the coefficient of volumetric expansion, a = 
(OV /OT)>/V, and the coefficient of isothermal compressibility 
B, = —(0V/dP),/V. 

Remembering that for any three variables XY, Y, Z 


(32), (sx), (SF), 


Equation [4] reduces to the form 
oT 


~°-(35),/ 


Similarly, from Equation [5] 


OV / p al al 
since C, > a PV for all the liquids with which we are con- 
cerned. 

Although C, is the physical quantity that is normally 
measured, the corresponding value of C, may be obtained 
from the fundamental relation a 


2 (x) 
oT os oP T Br 


Substituting from Equation [6] 


(x (2 Bl, /C, 
V\oP/, V a/av 
[8] 


Bz ( 4 , 
- — = — 
Cp Br Br Cp 


Finally, from Equations [7 and 8], the proportional changes 
in density (p) and specific volume in any process are given by 


p | V 


2VT E 


ll 


[5] 


e 

If we assume that the changes in the physical properties a, 
Cp, V, T and Br are negligible for a comparatively large 
change of E and P, Equation [9] becomes 


A AV oi 
ap (-6, ) + 
p J Cp 


Cy 


[10] 


a 
—1 
= 
te 
= [6] 
a 
r 
r 
| 
| 
O45 


temperature 7', K 
density 1/V, Ib/ft* 
coefficient of thermal expansion a, 1/ deg C 
specific heat C,, ft Ibf/lb deg C 
isothermal compressibility coefficient 87, 1/psi 
density change in 70 per cent efficient pump, per 
cent per 100 psi rise 
temperature change in 70 per cent efficient pump, 
deg C per 100 psi rise 
density change in lines, per cent per 100 psi fall 
temperature change in lines, deg C per 100 psi fall 


Liquid 


Oxy 


90 

71.4 

570 

13 xX 1076 


+0.07 0.285 
—0.036 —0.22 
+0. 165 +0.22 


+0.013 


+0.30 
—0.070 
+0.33 


‘Phe 


Density Change in a Steady Flow Pump 


In the case of steady, adiabatic flow through a pump, 
Equation [2] becomes 


AE = W — V-AP [11] 


since the term P; AV is V2AP. 

(It may be noted in passing that the form of this equation 
shows that pump characteristics should be based on the 
specific volume, or density, at outlet, not the normally meas- 
ured inlet value.) 

The work done JW’; in an ideal pump is given by the equation 
W, = V.AP, and the pump efficiency 7 is defined as the ratio 
W,/W. Thus aly 


W; 
AK = — —V.AP = V.AP —— [12] 
n n 


Substituting into Equation [10] to obtain the proportional 
change in density through the pump 


Cp 


a:V-AP 1 — 7 


” [13] 


- ( "+ ar) 6r| 
Cp ” 
Temperature Change in a Steady Flow Pump _ | ' 


The density change may also be expressed in terms of the 
temperature change AT since 


whence AV/V = a-AT — BrAP. Substituting from 
Equation [13] 


a-AT = AP E (= + ar) ~ Br |+ B,AP 


whence 


=- aT [14a] 
ti 
‘This equation could be rewritten in terms of AE from Equa- 


tion [12] as 
AT = (AE + aVT-AP)/C, [14b] 

The differences between Equations [13, 14a and 14b] and 
those applicable to changes in gases arise because the internal 
energy of a liquid is not a unique function of the temperature, 
as is assumed for a thermodynamic gas; and also because in 
the gas equations the term aPV equals R, the gas constant, 
which cannot be neglected in comparison with C,. 
Density Change in Feed Pipe Lines 

As the liquid delivered by the pump flows down the feed 
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lines, the pressure falls, and a further gain in internal energy 
oecurs. No external work is done or heat added, and since 
the proportional pressure change AP/P is much greater than 
the volume change AV/V, Equation [2] reduces to 


AE = —VAP 
(Note that the pressure rise in the pipes AP is negative.) 
Substituting into Equation [10] 


Ap AV ( aeVT aVAP 
— =— = AP(- -- 
Cp ) Cy 


= —AP at) + Br| 


Pp 


[15] 


Similarly, the temperature change in this process is again ob- 
tained from the relation 
AV 
aAT = V + B,AP = 
In a complete cycle in which the pump draws liquid from a 
tank and returns it to the same tank through an orifice, the 
overall temperature rise is obtained from Equations [14a 


VAP VAP W 
+ (1 —aT)=-— 
) Cy Cyn C, 


as one would expect. 
The corresponding density change is, as one would also 


expect, from Equations [13 and 15] a 


J 


AP [a —aT)+ 
Cp 


= 


Practical Applications 


Values for the physical constants of liquid oxygen and RP-1 
fuel were given in a paper by Chao.? Inserting these in the 
above equations, we get the values shown in Table 1 on the 
basis of a pump efficiency of 70 per cent. 

It will be seen that the density and temperature changes, 
although small, are not insignificant. Except in the case of 
water, for which a = 0, the density change in the pumps is 
much less than would have occurred on isothermal compres- 
sion and that in the lines is correspondingly greater. 


2 Chao, G. T. Y., ‘Isothermal Compressibility of Liquid Oxy- 
gen and RP-1,”’ ARS Journat, vol. 29, no. 3, March 1959, pp. 
199-203. 
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HEAT TRANSFER ... studies at JPL with a camera using a Kerr cell 
shutter taking photos at 20,000 frames per second were the first high- 
speed, high-resolution motion pictures successfully recording the action 
of nucleate boiling. 


L/QUID PROPELLANT SYSTEMS...were pioneered at JPL. Development 
work began in 1943 and led to the first practical rocket power-plant in the 
United States in which spontaneous ignition took place upon mixing of 
the oxidizer and fuel, 


MATERIALS RESEARCH AND TESTING .../is one of many supporting 
research programs under way at the Laboratory and are considered a 
“‘must'’, in providing needed data for engineers concerned with the 


SOLID PROPELLANT SYSTEMS. ..received momentous impetus in 1947 design and development of propulsion systems. 


with the successful flight of the Thunderbird, a test rocket. This JPL 
pioneering achievement demonstrated a new technique which has since 
revolutionized the field of solid propellant rockets. 


DEVELOPMENT .... of efficient rocket power plants involves large scale TESTING ... of rocket engines resulted in the establishment of a center 
testing and the application of knowledge from many scientific and engi- for recording rocket engine measurements when in 1948 the Lab estab- 
neering fields—thermodynamics, combustion, heat transfer, fluid me- lished the first system serving five engine test cells. This has now expanded 
chanics, and metallurgy. to a complex multi-channel system. 
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New 


Reaction propulsion device (2,875,577). 
: W. W. Odenkirchen, San Gabriel, Calif. 


} 2,675, S17 


casing of malleable sheet 
material and containing a_ propelling 
charge. One end is closed and the other 
is flattened to define an elongated outlet. 
An ignition member is located in the 
intermediate portion of the outlet. 

Device for controlling the flow direction 
from a reaction jet (2,875,578). M. 
Kadosch, F. G. Paris and R. H. Marchal, 
Paris, France, assignors to SNECMA Co. 


Tubular 


70 PRESIORE SOURCE 
yo 


2.675.876 | 
j 


Two pairs of opposite slot-like passages 
at 90 deg with each other extending along 
the periphery of a circular outlet. Each 
passage is in the form of an are. Oppo- 
site passages are selectively connected 
to sources of higher or lower pressures. 
Flameholding (2,875,580). T. 
K. Moy, and F. W. Henning, Merriam, 
Kans., assignor to W estinghouse E lectric 
Corp. 

First and second members pivotally 
connected at their upstream edges, and 
mounted for rotation in a plane trans- 
verse to the direction of flow for a high 
velocity combustible mixture. The flame- 
holder is disposed in the flow path for 
creating a degree of stagnancy sufficient 
to effect flame retention. 

Rocket motor (2,876,620). C. Kk. Wein- 
land and A. C. Ellings, Inyokern, Calif. 

Solid internal-burning propellant grain 
within a thin-walled motor tube. A 
restraining member, having vibration 
damping properties, at the forward end 
of the grain resiliently positions the grain 
and prevents leakage of gases between 
the outer surface of the grain and the 
motor tube. 


Eprror’s Nore: Patents listed above 
were selected from the Official Gazette of 
the U.S. Patent Office. Printed copies of 
patents may be obtained from the Com- 
oe of Patents, Washington 25, 

D. C., at a cost of 25 cents each: design 
patents, 10 cents. 
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Variable area jet engine inlet duct 
(2,876,621). R. C. Bogert and G. A. 
Olson, Anaheim, Calif., assignors to North 
American Aviation, Inc. 

Half-conical body rotatably mounted 
to form a semi-annular orifice between a 
lip and body. Rotation of the body out 
of the air flow path increases the inlet area. 
Static pressure measuring device 
(2,876,640). L. R. Beach Jr., K. Radey, 
D. J. Gildea and W. E. Strohmeyer, 
Fullerton, Calif., assignors to North 
American Aviation, Inc. 

Body with two flat outer surfaces and 
three internal chambers. One chamber 
has apertures leading to one outer surface 
and another chamber has apertures 
leading to the other outer surface. A 
pressure line from the third chamber leads 
to an indicator. 

Plumbing mirror (2,876,673). R. A. 
Hamilton, Glendale, Calif., assignor to 
Lockheed Aircraft Corp. 

Pool of light reflecting liquid at the 
bottom of socket. transparent 
cover is spaced between the pool and upper 
end of the socket. A layer of transparent 
liquid occupies the space between the 
cover and pool. The reflecting pool and 
liquid layers have different indexes of 
refraction. 

Telemetering transmitter for a projectile 
(2,877,452). A.V. Astin, Bethesda, Md., 
assignor to the U.S. Navy. 

Signaling device within a missile for 
transmitting to a remote location signals 
controlled by the functioning of missile 
components during flight. 
Airborne missile (2,876,677). 
and K. F. Finlay, Los Angeles, 
assignors to Northrop Aircraft, Inc. 


J. R. Clark 
Calif. 


\ 


2.876.677 


Body with a pair of hinged wings that 
pivot to and from a horizontal position, 
supported by an aircraft when the mem- 
bers are vertical. A time delay system 
changes the wings to horizontal after 
release. 

Circular wing aircraft (2,876,964). H. 
F. Streib, Chula Vista, Calif. 

Central opening in a circular wing with 
the airfoil section symmetrical about a 
vertical axis. An impeller in the opening 
draws air inward radially across the wing, 
and downward through the opening. 
Circular wing aircraft with ducted power- 
plant (2,876,965). H. F. Streib, Chula 
Vista, Calif. 

« Impeller in a central opening within a 
circular wing, mounted for movement 
through an angle of 90 deg. In vertical 
position, forward thrust is supplied; in 
horizontal position, vertical thrust is 
supplied. A shroud ring duct tilts with 
the impeller. 

Porous area-section flap (2,876,966). 

W. L. Cook, Mountain View, Calif., as- 
signor to the U.S. Navy. 


George F. McLaughlin, Contributor 


Means for applying suction to a porous 
area on the upper camber of an airfoil 
flap. Air is drawn through the porous 
area into the interior of the flap. 
Radiation tracker for aiming at center of 
multiple targets (2,877,354). A. F. Fair- 
banks and C. M. W olfe, Glendora, Calif., 
—" to North American Aviatic n, 

ne 

Amplification of the signal output from 
a scanning detector sensitive to radiation. 
Detector position compensates the sign:ils 
according to displacement of the detector 
from a predetermined reference position. 


Rocket gun (2,876,678). G. A. Lyon, 
Detroit, Miche 


2,676,676 


Missile launcher with a fixed barrel and 
a tubular cartridge support registrable 
with the breech of the barrel. A thrust 
block affords a thrust backing and center- 
ing means for the firing aspect of the 
launcher. 

Apparatus for measuring velocity and 
acceleration (2,877,415). S. D. Rollé, 
Philadelphia, Pa. 

System for determining acceleration of 
a mass ejected from an aircraft. Discrete 
voltage pulses are applied to a magnetic 
recording head at a constant frequency, 
and spacing of the pulses is recorded as a 
function of acceleration. 

Pyrotechnic circuit tester (2,877,419). 
W. J. East and H. M. Landers, Washing- 
ton, D. C. 

Device for recording the presence of 
electrical energy. A _ probe receiving 
energy greater than a_ predetermined 
magnitude ignites and changes color of 
pyrotechnic cartridges stored in the 
hand grip of the tester. 

Design for a space missile (183,835). Art 
Zola Fox, Verona, N. J., assignor to 
Thiokol Chemical Corp. 
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outstanding attribute of storable missiles fueled with — 


 DIMAZINE 
unsym-Dimethyihyarazine, UDMH 


THE STORABLE FUEL 


DIMAZINE helps to expand the orbit of stor- | DIMAZINE is stable and non-corrosive dur- 
age and launching sites for both tactical and _ing storage . . . is not shock sensitive . . . has 
mobile strategic weapons because it is easy to high thermal stability, low freezing point, 
ship and handle safely. minimum susceptibility to contamination and 


Missiles powered by DIMAZINE may be high compatibility with most metals and ap- 
transported fully fueled with only the normal _ propriate sealing materials. 
safety precautions applicable to flammable, 


moderately toxic industrial chemicals. In large 
missiles, transportation and erection problems able, high performance fuel that gives fast 


can be greatly simplified by shipping the mis- hypergolic starts, smooth, stable combustion 
sile “empty” and fueling at the launching site —_and easily-controlled shutdowns. 


for either instant use or ready storage. Both 
DIMAZINE and its storable oxidizers are DIMAZINE is amply available. We will be 
readily transported to the site in standard _ pleased to furnish trustworthy data on its sup- 


tank cars or trucks. ply logistics, properties and handling. 


Additionally, DIMAZINE, is a highly reli- 


Putting ldeas to Work 


FOOD MACHINERY AND CHEMICAL CORPORATION 
Westvaco Chior-Alkali Division 


General Sales Offices: 


dint otis) 161 E. 42nd STREET. NEW YORK 17 
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Book Reviews 


Dynamical Analogies, by Harry F. Olson, 
Van Nostrand Co., Princeton, 2nd 
edition, 1958, xi + 278 pp. $6.75. 

Reviewed by James J. KavzLaRIcH 
Worcester Polytechnic Institute 


“Dynamical Analogies’’ is an exhaustive 
study of three classical analogies— 
acoustical, electrical and mechanical vi- 
brating systems. The book was first 
published in 1943; the second edition, 
published in 1958, has been brought up-to- 
date by the additions of a chapter on 
magnetic analogies and two chapters on 
analogies in the fields of noise, distortion 
and feedback. 

The author is principally concerned with 
carrying the reader through a comparative 
study of an acoustical, electrical and two 
mechanical (rectilinial and rotational) 
systems. The text is subdivided into 
chapters representing the most elementary 
study of these four systems and progres- 
sively develops more complex aspects of 
the four systems, chapter by chapter. 
It assumes the reader’s familiarity with the 
elements of alternating current circuit 
theory, physics and ordinary differential 
equations. 

Of the three systems on which this book 
concentrates, the acoustical may be the 
most unfamiliar to the engineer. Formal 
training in acoustics at the undergraduate, 
and even the graduate, level is not very 
common. For the most part, it is likely 
that acoustics problems must be under- 
taken by engineers whose training has been 
in mechanical or electrical engineering. 
The importance of this book lies in its 
treatment of the acoustical vibrating 
system. 

The author is especially well qualified 
to discuss acoustical systems having been 
associated with the Radio Corporation of 
America, in charge of acoustical research, 
since 1934. His contributions to this field 
are extensive. 


Advanced Mechanics of Fluids, by D. W. 
Appel, P. G. Hubbard, L. Landweber, 
E. M. Laursen, J. S. McNown, H. 
Rouse, T. T. Siao, A. Toch and C. Ss. 
‘ih, edited by Hunter Rouse, John 
Wiley and Sons, Inc., New York, 1959, 
444 + xiv pp. $9.75. 
Reviewed by Victor L. STREETER 
University of Michigan 


“Advanced Mechanics of Fluids’’ is 
written as a text for the first or second year 
graduate student who is interested in in- 
compressible flow. He should have had 
at’least one course in fluid mechanics, and 
he should be well grounded in mathematics, 
The text is well-organized to cover a 
wide field of interest and has integrated 
recent research findings into a modern, 
unified treatment. 

The material is introduced by a chapter 
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Ali Bulent Cambel, Northwestern University, Associate Editor 


on research in fluid motion which dis- 
cusses analytical vs. experimental methods 
and develops dimensional considerations. 
This is followed by chapters on funda- 
mental concepts and equations, irrota- 
tional flow, conformal mapping, laminar 
motion, turbulence, boundary layers and 
free-turbulence shear flow. 

In a text of this length with its broad 
coverage, certain treatments must neces- 
sarily be abbreviated and incomplete, 
and in a few places the mathematical 
steps have been omitted. In spite of 
the many revisions of manuscript before 
publication, there is still a lack of uni- 
formity of viewpoint in the various 
chapters. 

The material has been carefully selected 
and illustrated with worked out examples. 
Two hundred and seventy-six problems 
are given, with answers for most of them. 
The text, with its broad coverage, fills 
an important place in fluid mechanics 
literature and should be very useful to 
the advanced student. 


Dynamics of Flight, by Bernard Etkin, 
John Wiley and Sons, Inc., New York, 
1959, 519 pp. $15. 

Reviewed by A. WILEY SHERWOOD 
University of Maryland 


This book shows the advantages of 
careful organization and thoughtful pres- 
entation of an involved subject. An 
unusual degree of clarity and readability 
are achieved by consistently proceeding 
from the general to the particular and 
from the simple to the more complicated. 
The author almost invariably charts his 
course verbally before making the mathe- 
matical excursion. 

The first third of the book covers 
static and dynamic stability in the classi- 
cal manner of Bryan. The use of the 


LaPlace transform and the transfer- 


function concept are then introduced 


nd used intermittently throughout the 


remainder of the book. A chapter de- 
voted to these techniques should prove 
valuable as a refresher to those who have 
had some manipulative experience with 
LaPlace transforms. Similarly, a chap- 
ter is included on machine computation, 
which covers analogue computers in some 
detail. A chapter is devoted to servo- 
mechanism theory and a portion of the 
appendix is devoted to vector analysis. 
This background material should all 
be valuable to supplement prior ex- 
perience. 

The treatment is/from the standpoint 
of the aeronautical engineer; i.e., mathe- 
matical analysis is used to achieve prac- 
tical results/rather than as an end in it- 
self. There are extensive, up-to-date 
bibliographical references and a few 
numerical examples to guide the reader. 
This book is highly recommended for 


i 


ment. 


graduate students and practicing engineers 


interested in flight stability. 4 | 


Heat and Mass Transfer, by E. R. G. 
Eckert and Robert M. Drake Jr., 
McGraw-Hill, New York, 1959, 530 -+ 
xiv pp. $12.50. 

Reviewed by R. L. Youne 
University of Tennessee 


“Heat and Mass Transfer’ is the second 
edition of the widely-used ‘Introduction 
to the Transfer of Heat and Mass” 
(1950). However, such extensive changes 
and additions have been made that the 
reviewer fully subscribes to Dr. Eckert’s 
statement in the preface, “Those familiar 
with the first edition will see that the 
major part of it was rewritten completely 
and that this edition can be considered 
in this respect practically a new book.” 
Additions have been made to the extent 
that the new edition contains approxi- 
mately twice as many pages as the old. 
Compared to the old “Heat and Mass 
Transfer’ represents a very successful 
attempt to present a more complete 
survey of the many techniques for solving 
problems of engineering interest while 
also providing a much improved book for 
classroom use. 

The addition of problems at the end of 
each chapter will be welcomed by those 
who plan to use the book as a text. As 
Dr. Eckert warns the reader in the 
preface, many of the problems are quite 
thought-provoking; thus the teacher 
would be well-advised to inspect them 
carefully before blithely assigning all 
of the problems at the end of chapter so 
and so for next week’s homework. As 
in the case of the first edition, the example 
problems are pertinent and well-chosen. 

The material on conduction was con- 
tributed by Prof. Drake and has been 
much expanded. In less than 100 pages, 
the reader is whisked from the simplest 
steady-state conduction situation to the 
timely problem of heat conduction with 
moving boundaries. In addition to the 
mathematical analysis of conduction prob- 
lems, information is given on the flux 
plot, the relaxation method and the 
Schmidt graphical method. 

The major portion of the book is devoted 
to heat transfer by convection. In 
approximately 200 pages, topics considered 
include: Flow along surfaces and in 
channels, forced convection in laminar. 
turbulent and separated flows, special 
heat transfer processes (including high 
velocity, rarified gas and liquid metal 
heat transfer and transpiration and film 
cooling), free convection and condensa- 
tion and evaporation. Throughout, basic 
principles are emphasized, and sufficient 
experimental results are presented to show 
comparisons between theory and experi- 
Contrary to the practice in the 
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immediate Openings in Advanced Areas for Engineers 
and Scientists at all Levels of Experience: 
ELECTRONICS: Inertial Guidance & Navigation e Digital Computer 
Development Systems Engineering Information Theory 
Telemetry-SSB Technique « Doppler Radar Countermeasures 
ee Radome & Antenna Design « Microwave Circuitry & Components 
¢ Receiver & Transmitter Design ¢ Airborne Navigational ae 
Systems ¢ Jamming & Anti-Jamming ¢ Miniaturization- 
Transistorization Ranging Systems Propagation Studies* 


Ground Support Equipment « Infrared & Ultra-Violet Techniques 


During the early years of this century the airplane was only THERMO, AERODYNAMICS: Theoretical Gasdynamics ¢ Hyper- 
the dream of a few dedicated men. Yet in the short span of Velocity Studies « Astronautics Precision Trajectories « Air Load 
5 decades this dream has evolved into such advanced aircraft and Aeroelasticity « Airplane/Missile Performance « Stability and 
as Republic’s F-105 — the free world’s most powerful Controls « Flutter & Vibration « Vehicle Dynamics and 


te ight in the } System Designs High Altitude Atmosphere Physics Re-entry 
— which is capable of Sight in the Mach 2 regime Heat Transfer Hydromagnetics « Ground Support Equipment 


The same holds true for missiles and space vehicles. Thirty brief PLASMA PROPULSION: Plasma Physics + Gaseous Electronics « 
age ae rom existed ae a few nee? Today at Hypersonics and Shock Phenomena « Hydromagnetics « Physical 

— act sa imaginations of many men are working to create Chemistry * Combustion and Detonation « Instrumentation « 
the vehicles that will allow man to explore the last frontier — High Power Pulse Electronics 
space. Included in this far-ranging research and development , : 
effort are plasma propulsion systems, electronic and hydraulic 
A eapons Effects « Radiation Environment in Space* Nuclear 
ens ages that will operate efficiently in a war Seas, Power & Propulsion Applications « Nuclear Radiation Laboratories 
and the calculation of super-accurate space flight trajectories. 
Working across the total technology of flight, Republic engineers Send resume in confidence to: 
and scientists see their ideas become realities because the novel, Mr. George R. Hickman 


the unique and the revolutionary in technical thinking are ; , 
Engineering E 10J 
appreciated and encouraged by management. New gineering Employment Manager, Dept 


investigations and new contracts mean you can put STEP 


your ideas in motion at Republic Aviation. FARMINGDALE, Lonc IsLAND, NEw YORK 
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A DIVISION 


OF VITRO 


CORPORATION 


OF AMERICA 


PREC 


ELECTRONICS SINCE 


NEMS -CLARKE 
PM-406 


PREAMPLIFIER- 
MULTICOUPLER 


Tue PM-406 is a combination instrument com- 
prised of a preamplifier and multicoupler located 
on a single chassis for rack mounting. This unit 
permits the coupling of eight receivers to a 
single antenna and is designed for use where 
short runs of cable from the antenna ore em- 
ployed. These units are of the “bath tub" chassis 


design having a 3-inch recessed front panel! to 


permit patching cables to hang without protrud- 


jing. Modifications of this design are available 


on special order. 


SPECIFICATIONS 


215-260 megacycles 

Within 3db 

22db (approximate) 

. Designed to operate in a 50-ohm system 


less than 4.5db 
37db minimum 
x 7” x 164%” 


Noise Figure 
Isolation between outputs 


919 JESUP-BLAIR DRIVE 


SILVER SPRING, MARYLAND 
MUNICIPAL AIRPORT 
MARTINSBURG, WEST VIRGINIA 
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first edition of basing virtually all con- 
vection heat transfer analysis on the von 
K4rman integral momentum and energy 
concept, the approach in the new edition 
is to start with general forms of the 
continuity, momentum and energy equa- 
tions and reduce them to expressions 
suitable for specific convection situations. 
The valuable integral concept is used in 
the new edition, but the present approach 
puts the entire discussion on a much 
sounder foundation. 

The material on thermal radiation has 
been extended to include topics on ra- 
diation pyrometry, radiation error in 
temperature measurements, the network 
method and solar radiation. Two chap- 
ters on the transfer of mass are included. 
Here again, the treatment is more ex- 
tensive than that in the first edition. In 
the final chapter titled Heat-Exchanger 
Calculations, both conventional and stor- 
age-type exchangers are discussed in some 
detail. Prof. Drake has again contributed 
an appendix of property values of com- 
pleteness sufficient for student use. 

In summary, this work is enthusiasti- 
cally recommended as a textbook and as a 
reference work for those engaged in the 
field of heat and mass transfer. 


Guided Missile Engineering, edited by 
Allen E. Puckett and Simon Ramo, 
McGraw-Hill, New York, 1959, viii 
497 pp. $10. 

Reviewed by JoHN GusTAVSON 
Grand Central Rocket Company 


This book represents a very good 
survey of the basic sciences which we 
usually include in missile system consider- 
ations. The two editors have understood 
the art of laying out the work, thereby 
providing continuity in the reading. 
Subsequently, they have managed to 
acquire good contributors for the 17 
chapters of the book. 

After a general discussion of the guided 
missile system, each author goes into 
detail outlining the fundamentals of 
aerodynamics, propulsion, navigation, in- 
strumentation and other related fields. 
The individual chapters serve as in- 
troductions to the particular fields, and 
since each contributor writes about 
his specialty only, the completed book 
offers a very solid introduction, indeed. 

The book is well-equipped with ref- 
erences, mostly to the classical works. 
The only serious gap is the exclusion of 
chapters on structures and materials for 
missiles, and electronic components. 
These subjects are of importance to a 
large part of our missile industry and 
have already become specialized fields. 

In general, the book is. well-illustrated 
and includes many excellent schematics 
and graphs. Quoting the editors, I be- 
lieve that the surveys of the various fields 
will assist each engineer in appreciating 
the nature of the problems encountered by 
his collaborators, and provide him with a 
better basis for the profitable discussion of 
mutual or interacting probiems. This 
book is certain to serve as a textbook 
for senior level courses in guided missile 
engineering. 
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October 1958, when the Thor-Able lunar probe soared 
79,000 miles, was a time of quiet pride for Clay Boyce. 
Design engineer Boyce was responsible for successfully 
predicting the in-flight performance of the Aerojet second 
stage of the Able vehicle. 

Clay Boyce has gone on to become an Aerojet Systems 
Division group leader, in charge of design and installation 
for the next generation of Able upper-stage vehicles for 


AZUSA AND NEAR SACRAMENTO, CALIFORNIA * A SUBSIDIARY OF THE GENERAL TIRE & RUBBER COMPANY 
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scientific and military applications. You’ll agree, a mighty 
important assignment for a BSME still in his twenties. 
Clay Boyce, with Aerojet since 1955, exemplifies the 
possibilities that exist at Aerojet for professionally gifted 
younger men to perform tasks of engrossing interest. 
We'd be delighted to hear from you. Write: Director of 
Scientific and Engineering Personnel, Box 296D, Azusa, 
California, or Box 1947D, Sacramento, California. 
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ENGINEERS 
and 
SCIENTISTS 


Research Opportunities 


Aeronutronic, a new divi- 
sion of Ford Motor Com- 
pany, has immediate need 
for qualified people to 
staff senior positions at 
its new $22 million Re- 
search Center in Ne rt 
Beach, Southern “Eali- 
fornia. 


The Space Technolo; 
Operation offers the 
highly desirable combina- 
tion of new facilities and 
advanced equipment, lo- 
cated in California’s finest 
environment for living 
and raising a family. In- 
vestigate these excep- 
tionally rewarding posi- 
tions now: 


VEHICLE TECHNOLOGY 
Aerodynamic design 
and testing 
Rocket engine 
development 
Rocket nozzle and re- 
entry materials 
High temperature 
chemical kinetics 
Combustion and 
detonation theory 
Combustion 
thermodynamics 
High temperature 
structural plastics 
& ceramics 
Advanced structures 
Rocket vehicle 
systems 


MISSILE DEFENSE 

Supersonic 
aerodynamics 

Aerothermodynamics 

High temperature heat 
transfer 

Space physics 

Re-entry programs 


ASTRO SCIENCES 


Space electronics 

Guidance & control 

Communications 

Instrumentation 

Experimental physics 

Plasma and magneto- 
hydrodynamics 
studies 


Qualified applicants are invited 
to send resumes and inquiries to 
Mr. R. W. Speich, Aeronutronic, 
Dept. 19, Box 451, 
Beach, California. 


AERONUTRONIC 


a Division of 


Newport 


FORD MOTOR COMPANY 
Newport Beach 
Santa Ana Maywood, California \ 


Technical Literature Digest—— 
M. H. Smith, Associate Editor, and M. H. Fisher, Contributor 
The James Forrestal Research Center, Princeton University 


Note on the Scavenge of 6-inch-diam- 
eter Ram-Jet Exhaust in Mach 3.1, 
1- by 1-foot Supersonic Tunnel, by 
Warren C. Burgess Jr. and L. Eugene 
Baughman, NACA RM E54K30, Feb. 
1955, 15 pp. (Declassified by authority 
of NASA Pub. Ann. no. 3, 2/4/59, 
p. 23.) 

Preflight Tests and Flight Performance 
of a 6.5-inch-diameter Ram-jet Engine, 
by Arthur H. Hinners Jr. and Douglas H. 
Foland, NACA RM L53H28, Nov. 1953, 
55 pp. (Declassified by authority of 
NASA Pub. Ann. no. 3, 2/4/59, p. 27.) 


Summary of Free-flight Performance of 
a Series of Ram-jet Engines at Mach 
Numbers from 0.80 to 2.20, by Warren 
J. North, NACA RM E53K17, Feb. 
1954, 37 pp. (Declassified by authority 
of NASA Pub. Ann. no. 3, 2/4/59, p. 21.) 

Effect of Angle of Attack and Exit 
Nozzle Design on the Performance of a 
16-inch Ram Jet at Mach Numbers from 
1.5 to 2.0, by Eugene Perchonok, Fred 
Wilcox and Donald Pennington, NACA 
RM E51G26, Oct. 1951, i (De- 
classified by authority of NASA Pub. 
Ann. no. 3, 2/4/59, p. 19.) 

Investigation at Zero Angle of Attack 
of a 16-inch Ram Jet Engine in 8 by 6 
Foot Supersonic Wind Tunnel, by 
Nussdorfer and E. Perchonok, NACA 
Res. Mem. E50L04, March 1951, 33 pp. 
(Declassified by authority of NASA 
Pub. Ann. no. 3, 2/4/59, p. 19.) 

Free Jet Investigation of a 16-inch Ram 
et at Mach Numbers of 1.35, 1.50and 1.73, 
y Fred Wilcox, Sol Baker and Eugene 

Perchonok, NACA Res. Mem. E50G19, 
Oct. 1950, 59 pp. (Declassified by authority 
of NASA Pub. Ann. no. 3, 2/4/59, p. 19.) 

Application of Gas Analysis to Com- 
bustor Research, by R. R. Hibbard and 
Albert Evans, NASA Memo 1-26-59E, 
Feb. 1959, 22 pp. 

Prophecy for Propulsion—One Man’s 
Guess, by T. F. Dixon, Aero-Space 
Engng., vol. 18, Feb. 1959, pp. 35-39. 

Problems of Research in Annular 
Combustion Chambers, by E. Caffier, 
Gesellschaft fir Luft- 
fahrit, Jahrbuch, 1957, Braunschweig, 
Friedr. Vieweg & Sohn, 1958, 512 pp., 
pp. 410-416. (In German.) 

Propellant Vaporization as a Criterion 
for Rocket-engine Design; Experimen- 
tal Performance, Vaporization, and Heat- 
transfer Rates with Various Propellant 
Combinations, by Bruce J. Clark, Martin 
Hersch and Richard J. Priem, NASA 
Mem. 12-29-85E, Jan. 1959, 31 pp. 

Optimum Burning Time of High Alti- 
tude Rockets, by L. E. Giogieri, Aero- 
tecnica, vol. 38, Aug. 1958, pp. 195-209. 
(In Italian.) 


Eprror’s Nore: Contributions from Pro- 
fessors E. R. G. Eckert, J. P. Hartnett, T. 
F. Irvine Jr. and P. J. Schneider of the 
Heat Transfer Laboratory, University of 
Minnesota, are gratefully acknowledged. 


Aerodynamics of Jet 
Propelled Missiles 


Hypersonic Viscous Flow over Slender 
Cones, by L. Talbot, T. Koga and ?’. M. 
Sherman, Inst. Aeron. Sci., Rep. 591, 
Jan. 1959, 21 pp. 

An Experimental Investigation of the 
Applicability of the Hypersonic Similarity 
Law to Bodies of Revolution, by Stanford 
E. Neice and Thomas J. Wong, NACA 
RM A52K07, Jan. 1953, 32 PP. (De- 
classified by authority of NASA Pub. 
Ann. no. 3, 2/4/59, p. 18.) 

The Cone Sphere in Hypersonic He- 
lium above Mach Number Twenty, 
by Robert H. Johnson, Inst. 
Sci., Rep. 59-2, Jan. 1959, 9 pp. 

Supersonic Drag and Base Pressure 
of a 70° Cone Cylinder, by G. D. Kahl, 
Aberdeen Proving Ground, Ballistic Res. 
Labs., Mem. Rep. 1178, Nov. 1958, 42 pp. 

Determining Air Reactions on Moving 
Vehicles, Part II: Methods of Rocketry, 
by M. Z. Krzwoblocki, Wright Air Deo. 
Center, Tech. Rep. 56-51 (ASTIA AD 
161004), July 1958, 92 pp. 

Two Non-liner Problems in the Flight 
Dynamics of Modern Ballistic Missiles, 
by John D. Nicolaides, Inst. Aeron. 
Sci., Rep. 59-17, Jan. 1959, 44 pp. 

Boundary Layer Transition on Blunt 
Bodies with Highly Cooled Boundary 
Layers, by Kenneth F. Stetson, /nst. 
Aeron. Sci., Rep. 59-36, Jan. 1959, 27 pp. 

The Cone-sphere in Hypersonic He- 
lium above Mach 20, by Robert H. 
Johnson, Aero-Space Engng., vol. 18, 
Feb. 1959, pp. 30-34, 53. 

Performance Measurements from a 
Rocket-powered Exploratory Research 
Missile Flown to a Mach Number of 10.4, 
by Robert O. Pilaud, NACA RM L54L29a, 
March 1955, 12 pp. (Declassified by 
authority of NASA Pub. Ann. no. 3, 
2/4/59, p. 28.) 

Flow in the Shock Layer Region of 
Supersonic and Hypersonic Blunt Bodies: 
(1) A Survey, by C. A. Andrade, Chrysler 
Corp., Missile Div., Tech. Rep. AME- 
RIJ, June 1958, 59 pp. 

A Note on the Interpretation of Base 
Pressure Measurements in Supersonic 
Flows, by R. C. Hastings, Gt. Brit., 
Aeron. Res. Council, Current 
409 (formerly ARC Tech. Rep. ; 
Roy. Aircr. Estab., TN Aero 2572), June 
1958, 13 pp. 

Some Variational Solutions to Rocket 
Trajectories over a Spherical Earth, 
by Angelo Miele and James O. Cappellari 
Jr., Purdue Univ., School of Aeron. 
Engng., Rep. A-48-9 (AFOSR-TN-58-1012; 
ASTIA AD 206155), Nov. 1958, 29 pp. 

Aerodynamic Studies in the Aero- 
ballistics Range, by G. V. Bull, Wissen- 
schaftlichen Gesellschaft fur Luftfahrt, 
Jahrbuch, 1957, Braunschweig, Friedr. 
Vieweg & Sohn, 1958, pp. 247-260. 

Development of Aeroballistic Ranges, 
by R. E. Kutterer, Wissenschaftlichen 
Gesellschaft fiir Luftfahrt, Jahrbuch, 1957, 
Braunschweig, Friedr. Vieweg & Sohn, 
1958, pp. 261-278. (In German.) 
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Bs, MS, PhD 


at 2000° F to 5000° F in Nuclear Flight Program 


Ceramic problems facing engineers and 
scientists at General Electric’s Aircraft 
Nuclear Propulsion Department equal — 
and transcend —the materials require- 
ments created by missile re-entry 
conditions. 
Consider this contrast: a nose cone 
_ plunging to earth at 15,000 mph reaches 
extreme temperatures for less than one 
minute... but ceramic components envis- 
_aged for nuclear flight will encounter 
_ operating cycles in the 2000°F to 5000°F 
range for periods of hours or days. 
Current opportunities for ceramists in- 
clude both applied research and materials 
development: 


® LEAD TECHNICAL DIRECTION of major mate- 
rials development activity utilizing ceramics. PhD. 


@ SUPERVISE APPLIED CERAMIC RESEARCH 
_ including fabrication processes, testing, develop- 
~ ment of novel ceramic materials. PhD. 


= CONDUCT RESEARCH in ceramic technology in 
_ particular fields of fabrication, chemical reactivity, 
high temperature testing, development of high 


strength ceramic materials. MS, BS. 
CERAMISTS who value the opportunity to do original work 


TAKE with a company that fosters free inquiry and initiative, are 
'..-including financial estimates... for major ce- invited to inquire about positions now open in the areas to the 


PROVIDE entire Applied Materials Research op- wes Me. PW. Chas 
eration with scientific consulting advice in the 
fields of ceramics and physical chemistry related 


_ramic materials development activity. PhD, MS, BS. left. Please include salary requirements with resume. 


AIRCRAFT NUCLEAR PROPULSION DEPARTMENT : 


plus graduate study essential.) 
P.O. Box 132 Cincinnati 15, Ohio 
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by NEWBROOK 


Application of Power Spectral Methods 
in Airplane and Missile Design, by Ken- 
neth R. Thorson and Quentin R. hne, 
a Aeron. Sci., Rep. 59-42, Jan. 1959, 
4 PP. 


MOTOR 


Solid and Liquid 


PLENUM CHAMBERS 


Heat Fluid 
Flow 


Propellants 


JATO CASES 
NOZZLES 


BLAST TUBES 
FUEL INJECTORS 


The newest addition to the Quality Control facilities at Newbrook 


is the Hydrostatic Test Cell illustrated below. All controls are on 


the outside. A T.V. Camera inside the cell enables the engineers 


to watch the test on a T.V. screen. This is only one of many 


projects at this modern plant manned and equipped to produce 


the finest in missile components. 


Finest Welding Facilities 
Certified Welders 


Below: Machining Motor Cases 


Temperature Measurements of Shock 
Waves by the Spectrum-line Reversal 
Method, by J. G. Clouston, A. G. Gaye 


don, F. R. S., and I. I. Glass, Proc. /’oy. 
Soc. London, vol. A248, Dec. 9, 1458, 
pp. 429-441. 


Magnetohydrodynamics Opens up New 
Electronic Vistas, by Irsin Stamler, 
vol. 31, Jan. 1959, pp. 28- 

Conversion of Decibels to 
Albert T. Isaacs, Space/Aeron., vol. 
Jan. 1959, p. 73. 


An Experimental Investigation of the 
Stability of the Hypersonic Laminar 
Boundary Layer, by Anthony Demetria:les 
Calif. Inst. Tech., Guggenheim Aeron. 
Lab., Hupersonic Res. Proj., Mem. 43, 


May 1958, 85 pp. 

Transport and Thermodynamic Proper- 
ties in a Hypersonic Laminar Bound 
Layer, by Sinclaire M. Scala and Charles 
W. Baulknight, General Electric Co., 
Missiles and Ordn. Sustems Dept., Docu- 
ment 588d232 (Aerophys. Res. Mem. 10), 
April. 1958, 49 pp. 

Interaction of Slow Electrons with 
Atomic Oxygen and Atomic Nitrogen, by 
Milton M. Klein and Keith A. Brueckner, 
General Electric Co., Missile and Ordn. 
Systems Dept., Document 58Sd222 (Aero- 
phys. Res. Mem. 11), April 1958, 17 pp. 

figs. 

Cooling of JP-4 Fuels in Electrically 
Heated Metal Tubes, by Arthur L. Smith, 
William P. Cook and Vincent F. Havin, 
NACA Res. Mem. E56H21, Nov. 1956, 
17 pp. (Declassified from Confidential 
by authority of NASA PA 1, p. 16, 
11/14/58.) 

The Heat Transfer and Flow Friction 
Characteristics of Three Double and 
Triple Sandwich High Performance Heat 
Transfer Surfaces, by W. H. Kays, 
Stanford Univ., Dept. Mech. Engnq., 
Tech. Rep. 37, Oct. 1958, 42 pp. 

The Thickness of Cylindrical Shocks 
and the PLK Method, by H. C. Levey, 
Australia. Aeron. Res. Labs., Rep. A 1-8, 
Jan. 1958, 22 pp., 5 figs. 

Some Methods of Evaluating Imperfect 
Gas Effects in Aerodynamic Problems, by 
G. A. Bird, Gt. Brit., Aeron. Res. Council. 
Curr. Paper 397. (Formerly ARC Tech. 
Rep. 19547; Roy. Aircr. Estab., Tech. 
Note Aero 2488), 18 pp., 13 figs. 

The Exact Flow Behind a Yawed Conical 
Shock, by G. Radhakrishnan, Cranfield. 
Coll. Aeron., Rep. 116, 22 pp., 15 figs. 

On the Mixing of Two Parallel Streams, 
by Lu. Ting, Brooklyn Polytech. Inst., 
Dept. Aeron. Engng. & Appld. Mech., 
Rep. 441 (AFOSR TN 58-628; ASTIA 
AD 162158), July 1958, 23 pp., 3 figs. 

Structure of Shock Fronts in Various 
Gases, by W. H. Anderson and D. F. 
Hornig, Brown Univ., Metcalf Res. Lal, 
Tech. Rep. 8, May 1958, 12 pp. 

Shock Formation from Strong Com- 
pression Waves, by Reynold A. Shunk, 
Lehigh Univ., Inst. Res., Tech. Rep. 11, 


July 1958, 60 pp. 
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WYMAN-GORDON !S FORGING RE-ENTRY SHIELDS 


e Copper in production 
e Beryllium in limited production 
e Reinforced plastics in development 


FORGING 
STEEL TITANIUM BERYLLIUM 
AND OTHER UNCOMMON MATERIALS 


WORCESTER, MASSACHUSETTS 


et 
A: 
HARVEY, ILLINOIS GRAFTON. MASSACHUSETTS FRANKLIN PARK, ILLINOIS 
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BOEING 
SCIENTIFIC RESEARCH 
LABORATORIES 
announces 
supervisory and senior 
staff positions in its 
GAS DYNAMICS 
LABORATORY 


Qualifying scientists are provided 


with exceptional opportunities to 


work on problems of their choice 
and to publish their findings. 
Positions are available now for 


physicists, physical chemists and 


aerodynamicists to conduct inves- 


tigations in the general areas of 


combustion, detonation, shock 


dynamics, high temperature 


gaseous reactions and magneto- 


hydrodynamics. 


For further information, write: 


Dr. Y. A. Yoler 
Head, Gas Dynamics 
Laboratory 
Boeing Airplane Company 
P.O. Box 3822 - JTC 
Seattle 24, Washington 


EMR 


.. environment for 
dynamic career growth 


Measured Heat Flow to the Wall of a 
Shock Tube, by W. Paul Thompson, 
Lehigh Univ., Inst. Res. Tech. Rep. 10, 
July 1958, 36 pp. 

A Study of the Processes Induced by the 
Reflection of Strong Shocks in Nitrogen, 
by Robert G. Jahn and Fred A. Grosse, 
Lehigh Univ., Inst. Res., Tech. Rep. 13, 
Sept. 1958, 38 pp. 

Measurement of Wall Heat Transfer 
and of Transition to Turbulence during 
Hot Gas and Rarefaction Flows in a 
Shock Tube, by Albert J. Chabai, Lehigh 
Univ., Inst. Res., Tech. Rep. 12, Sept. 1958, 
86 pp. 

ypersonic Heat Transfer to Surfaces 
Having Finite Catalytic Efficiency, by 
Sinclaire M. Scala, Gen. Electric Co., 
Missile & Ordn. Systems Dept., Document 
57SD646 (Aerophys. Res. Mem. 4), 
July 1957, 24 pp., 12 fig. 

Dynamic Stabilization of Hydrodynamic 
and Hydromagnetic Instabilities, by How- 
ard D. Greyber, Gen. Electric Co., Missile 
& Ordnance Systems Dept., Document 
R58SD208 (Aerophys. Res. Mem. 7), 
Jan. 1958, 23 pp. 

A Shock Tube Investigation of Heat 
Transfer in the Wake of a Hemisphere 
Cylinder, with Application to Hypersonic 
Flight, by William E. Powers, Kenneth F. 
Stetson and Mac C. Adams, AVCO 
Mfg. Corp., AVCO Res. Lab., Res. Rep. 
30, Aug. 1958, 39 pp. 

On the Instability Thecry of the Soft 
or Melted Surface of an Ablating Body 
When Entering the Atmosphere, by 
Saul Feldman, AVCO Mfg. Corp., AVCO 
Res. Lab., Res. Rep. 34, Aug. 1958, 38 


Electric and Gas Heaters for Re-entry 
Simulation and Space Propulsion, by 
Thomas R. Bragan, AVCO Mig. Corp., 
AVCO Res. Lab., Res. Rep. 35, Sept. 
1958, 29 pp. 

Radiative Transport within an Ablating 
Body, by Leo P. Kadanoff, AVCO Mfg. 
Corp., AVCO Res. Lab., Res. Rep. 37, 
Oct. 1958, 30 pp. 

On Flow of Electrically Conducting 
Fluids over a Flat Plate in the Presence of 
a Transverse Magnetic Field, by Vernon 
J. Rossow, NACA Rep. 1385, 1958, 20 
pp. (Supersedes 7'N 3971.) 


Combustion, Fuels and 
Propellants 


Thermal Stresses in Long Cylindrical 
Case Bonded Propellant Grains, by Eric 
KE. Unger and Bernard W. Shaffer, NY 
Univ., Coll. Engng., Res. Div., Sept. 1958, 
30 pp. 

Experimental Evaluation of Boron- 
hydrocarbon Slurry in a 16-inch Ram-jet 
Combustor, by William B. Kerslake, 
E. E. Dangle and A. J. Cervenka, NACA 
Res. Mem. E55C07, June 1955, 37 pp. 
(Declassified from Confidential by author- 
ity of NASA PA 1, p. 15, 11/14/58.) 

Review of the Toxicological Properties 
of Pentaborane, Diborane, Decaborane, 
and Boric Acid, by Joseph M. Lamberti, 
gen Res. Mem. E56H13a, Dec. 1956, 

(Declassified by authority of 
A PA, 1, p. 16, 11/14/58.) 

pew ae of NACA Research on Igni- 
tion Lag of Self-igniting Fuel—Nitric 
Acid Propellants, by Gerald Morrell, 
NACA Res. Mem. E57G19, Oct. 1957, 
48 pp. (Declassified from Confidential 
by authority of NASA PA 1, p. 16, 


Infrared Method for Determining 2,4- 
and 2,6-dinitrotoluene; 2,4 ,6-trinitrotol- 
uene; and o- and nitrotoluene in 
Admixture, by Walter E. Fredericks and 
Frank Pristera, Picatinny Arsenal, Felt- 
man Res. and Engng. Labs., Tech. Rep. 
2546, Oct. 1958, 9 pp. 

Study of Supersonic Flows with Chem- 
ical Reactions: Measurement of Rate Con- 
stants of the Reaction N. + N.04 s N, + 
2NO. at Low Reactant Concentration, 
by Peter P. Wegener, Jack E. Marte 
and Carl Thiel, Calif. Inst. Tech., Jet 
Prop. Lab., Progr. ep. 20-349, ‘J ine 
1958, 26 pp. 


Materials of Construction 


On the Linear Relation between the 
Softening Temperature and the Melting 
Point of Ceramics, by George A. Hoff - 
a. William J. Knapp, Rand Corp., 

2263 (ASTIA AD 306554), Oct. 
12 pp. 

Materials for Space Flight, by George 
A. Hoffman, Rand Corp., P-1420, July 
1958, 26 pp. 

Roll-spot Welding for Ballistic Mis- 
siles, by James K. Dawson, Welding J/., 
vol. 38, May 1959, pp. 475-484. 

Stability of Ceramics in Hydrogen 
Between 4000° and 4500° F, by Charles 
E. May, Donald Koneval and George C. 
Fryburg, NASA Memo 3-5-59E, March 
1959, 13 pp. 

Reinforced Plastics in High Tempera- 
ture Applications, by Edwin H. Jaffe, 
Industrial Labs., vol. 10, March 1959, 
pp. 67-68, 70-72, 74. 

Aeronautical and Space Applications 
for Ceramics, by James D. Welterlen, 
Aero/Space Engng., vol. 18, April 1959, 
pp. 72-75. 

Molybdenum Coating Stands up at 
2000 Deg F, by Irwin Stambler, Space/- 
Aeron., vol. 31, April 1959, pp. 97, 98, 100. 


Instrumentation, Data 
Recording, Telemetering 


Diffusion Effects in the Transpiration 
Method of Vapor Pressure Measurement, 
by Ulrich Merten, J. Phys. Chem., vol. 
63, March 1959, pp. 443-445. 

Interior Ballistic Telemetering, by 

. M. Kendrick and L. A. Peters, 
Aberdeen Proving Ground, Ball. Res. 
Labs., Mem. Rep. 1183, Dec. 1958, 29 pp. 

Transistors Telemeter Trans- 
mitter, by Enemark, Electronics, 
vol. 32, Mareh 13, 1959, pp. 136-137. 

New Power Sources for Space-age 
Electronics, by D. Linden and F. 
Daniel, Electronics, vol. 32, March 20, 
1959, pp. 43-47. 


Guidance Systems and 
Components 


Investigation of the Possibility of 
Simplifying Missile Guidance Systems 
by the Use of Free-floating Flaps and 
Spring-mounted Control Surfaces, by 

atsumi Kikido, Paul H. Hayashi and 
Henry C. Lessing, NACA RM A55L09, 
May 1956, 72 pp. (Declassified from 
Confidential by authority of NASA 
Pub. Ann. 4, p. 17, 3/16/59.) 

Application of Nonlinear Aerodynamic 
Stability Characteristics to Simplify Mis- 
sile Control and Guidance Systems, 
by Howard J. - Curfman Jr., NACA RM 
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C8, 
Du Pont Lino-Writ 4 is the fastest, toughest, whitest photorecording 
ge paper you can use. It accommodates writing speeds as high as 5000 
: cps at maximum amplitude; has exceptional exposing and processing 
’ latitude; its ultra-thin, all-rag stock has unusual wet or dry strength. 
Our other outstanding papers include: 
Lino-Writ 1. Orthochromatic paper designed for test requirements 
at lower writing speeds. 
Lino-Writ 2. Moderately high-speed orthochromatic paper for use in 
of intermediate speed recording requirements. 
i Lino-Writ 3. For high-frequency oscillographic traces without loss 
y of detail. 
‘ These three papers are available in two types: Type B, standard- 
7 Ths advertisement was prepared exclusively by Phototypography 


' 


SrPTEMBER 1959 


Better Things for Better Living... through Ch 


4 
weight stock when greater opacity is re- 
quired and Type W, a thin 100% all-rag stock. 


Lino-Writ Rapid Processing Chemicals 
Kit is Du Pont’s high-speed, low-odor for- 
mulation designed for rapid stabilization 
processing of all Du Pont photorecording 
papers. 

For more information on our complete 
photorecording line, write: E. I. du Pont 
de Nemours & Co. (Inc.), Photo Products 
Department, 2430-A Nemours, Wilming- 
ton 98, Delaware. In Canada: Du Pont of 
Canada Limited, Toronto. 
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orbiting 
around 


paddles 


capturing 


power 


, The scientific data that will some day enable 
7 us to probe successfully to the very fringes of 
4% the universe is being recorded and transmitted 
at this moment by the space laboratory 
Explorer VI, a satellite now in orbit around 
the earth @ This project, carried out by Space 
Technology Laboratories for the National 
Aeronautics and Space Administration under 
the direction of the Air Force Ballistic Missile 
Division, will advance man’s knowledge of: 
The earth and the solar system...The magnetic 
field strengths in space ...The cosmic ray 
intensities away from earth...and, 
The micrometeorite density encountered in | 
inter-planetary travel @ Explorer VI is the" 
- most sensitive and unique achievement ever 
launched into space. The 29” payload, 
STL designed and instrumented by STL in 
cooperation with the universities, will remain 
“vocal” for its anticipated one year life. 


7 Space Technology 


How ? Because Explorer VI’s 132 pounds of 
electronic components are powered by storage 
batteries kept charged by the impingement 
of solar radiation on 8,000 cells in the four 
sails or paddles equivalent to 12.2 square feet 
in area @ Many more of the scientific and 
technological miracles of Explorer VI will be 
reported to the world as it continues its epic 
flight. The STL technical staff brings to this 
space research the same talents which have 
provided systems engineering and over-all 
direction since 1954 to the Air Force Missile 
Programs including Atlas, Thor, Titan, 
Minuteman, and the Pioneer I space probe 


Important staff positions in connection with 
these activities are now available for 
scientists and engineers with outstanding 
capabilities in propulsion, electronics, 
thermodynamics, aerodynamics, structures, 
astrophysics, computer technology, and 
other related fields and disciplines. 


Inquiries 

are 
Laboratories, Inc. 
P.O. Box 95004 

Los Angeles 45, California 
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L54L06a, Feb. 1955, 23 pp. (Declas- 
sified from Confidential by authority of 
NASA Pub. Ann. 4, p. 25, 3/16/59.) 

A Discussion of Space Vehicle Guidance 
Problems, by F. T. Smith, Rand Corp., 
P-1568, Dec. 1958, 34 pp. 

Doppler Radar Navigation, by F. B. 
Berger, Electronics, vol. 32, May 8, 1959, 

62. 


Orientation and Control, by T. B. 
Garber, Rand Corp., P-1430, Feb. 1958, 
22 pp. 

Improvement of Radar Tracking, by 
A. I. Hastings and J. E. Meade, Naval 
Res. Lab., Rep. R-3424, Feb. 1949, 7 

. (Declassified from Confidential by 
authority of NRL ltr. Ser. 11214, dtd. 
Nov. 18, 1958, signed H. W. Ottenstroff.) 

The Synthesis of Optimum Homing 
Missile Guidance Systems with Statistical 
Inputs, by Elwood C. Stewart and Gerald 
L. Smith, NASA Memo 2-13-59A, 
April 1959, 57 pp. 

The Effects of Target and Missile 
Characteristics on Theoretical Minimum 
Miss Distance for a Beam-rider Guidance 
System, by Frank Druding and Togo 
Nishiura, NASA Memo 2-12-59A, March 
1959, 32 pp. 


Servomechanisms and 
Controls 


Spectral Methods for the Study of 
Servo Circuits Defined by Linear Differen- 
tial Equations with Stochastic Coeffi- 
cients, by Mauritz Sundstrom, Sweden, 
Roy. Inst. Tech., Dept. Aeron. Sci., TN 
48, 1958, 11 pp. 

Magnetic Amplifiers for Servo Sys- 
tems, by S. Davis, Electronics, vol. 32, 
March 13, 1959, pp. 134-135. 

Ground Simulator Studies of a Non- 
linear Linkage in a Power Control Sys- 
tem, by Arthur Assadourian, NASA 
Memo 2-15-59L, April 1959, 18 pp. 

A Class of Optimum Control Systems, 
by C. W. Merriman III, J. Franklin Inst., 
vol. 267, April 1959, pp. 267-281. 


Flight Vehicle Design and 
Testing 


Modern Launch Facility, by Warren 
F. Opitz, ARS Preprint 745-59, March 
1959, 12 pp. 

Real-time Analysis—New Approach in 
Flight Testing, by Guenther Hintze, 
ARS Preprint 764-59, March 1959, 15 pp. 

Project Space Track at Air Force 
Missile Test Center, by Walter H. Mann- 
ing Jr., ARS Preprint 766-59, March 
1959, 7 pp. 

Missile Countdowns, by John S. Harri- 
son, ARS Preprint 753-59, March 1959, 
11 pp. 

Thor Countdown Concepts, by T. J- 
Gordon, ARS Preprint 754-59, March 
1959, 12 pp. 

Automatic Countdown Checkout for 
Defensive Missile Systems, by Herbert 
W. Bethel, ARS Preprint 755-59, March 
1959, 3 pp. 

Missile Handling Equipment for Flight 
Test Programs, by C. M. Wendertz, 
ARS Preprint 759-59, March 1959, 8 pp. 

Philosophy of Prelaunch Preparation 
for Research and Development Missiles, 
by Robert E. Moser, ARS Preprint 760- 
59, March 1959, 7 pp. 

Propulsion Ground Equipment, by 
Andrew Pickett, ARS Preprint 746-59, 
March 1959, 5 pp. 
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Use of Tolerance Limits in Missile 
Evaluation, by Albert Madansky, Rand 
Corp., P-1423, July 1958, 13 pp. 

Designing a Jet Deflector for the Hawk 
Launcher, by J. A. Quinville, Space/ 
Aeron., vol. 31, April 1959, pp. 65, 67, 71. 

New Light Construction in the Guided 
Oerlikon Rocket, by F. B. Stencel, 
Fusees et Recherche Aeron., vol. 3, no. 3, 
1958, pp. 143-151. (In French.) 

A Simulation Study of a Wingless 
Missile, by Henry C. Lessing and David 

. Reese Jr., A RM _ A55L06, 
Feb. 1956, 18 pp. (Declassified from 
Confidential by authority of NASA 
Pub. Ann. 4, p. 17, 3/16/59.) 

Missile Components Tested at New 
Gaseous Oxygen Laboratory, by A. L. 
Spangenberg, Industrial Labs., vol. 10, 
March 1959, pp. 100-104. 

Topics in Dynamic Programming for 
Rockets, by Angelo Miele and James O. 
Cappellari Jr., Zeitschrift fiir Flugwissen- 
schaften, vol. 7, Jan. 1959, pp. 14-21. 
(In English.) 


Space Flight 


Man Machine Integration in Space 
Vehicles, by George W. Hoover, Jnst. 
Aeron. Sci., Rep. 59-24, Jan. 1959, 14 pp. 

The Problem of Escape from Satellite 
Vehicles, by C. V. Carter and W. W. 
Huff Jr., Inst. Aeron. Sci., Rep. 59-41, 
Jan. 1959, 19 pp. 

Systems Using Solar Energy for Aux- 
iliary Space Vehicle Power, by Albert 
E. von Doenhoff and James M. Hallissy 
Jr., Inst. Aeron. Sci., Rep. 59-40, Jan. 
1959, 20 pp. 

Studies of Pilot Control During Launch- 
ing and Reentry of Space Vehicles, Utiliz- 
ing the Human Centrifuge, by C. H. 
Woodling and C. C. Clark, Inst. Aeron. 
Sci., Rep. 59-39, Jan. 1959, 16 pp. 

Wave Propagation through Ionized 
Gas in Space Communications, by F. J. 
Tischer, Inst. Aeron. Sci., Rep., 59-34, 
Jan. 1959, 13 pp. 


Aerophysics, Astrophysics 


Radar Observations of the Second 
Russian Earth Satellite (Sputnik II 
19578), by J. G. Davies, J. V. Evans, 
8. Evans, J. 8S. Greenhow, J. E. Hall, E. 
L. Neufeld and J. H. Thomson, Proc. Roy. 
Soc., vol. A 250, March 24, 1959, pp. 367- 
376. 

Meteorites and the Origin of Life, 
by M. H. Briggs, Spaceflight, vol. 2, 
April 1959, pp. 39-43. 

The Extent of the Lunar Atmosphere, 
by Geoffrey Turner, Spaceflight, vol. 2, 
April 1959, pp. 57-58. 

Minor Planets, by G. E. B. Stephenson, 
Spaceflight, vol. 2, April 1959, pp. 59-60. 

The Rotational Period and Tilt of the 
Planet Venus, by F. C. Wykes, Space- 
flight, vol. 2, April 1959, pp. 63-64. 


Atomic Energy 


Problems in Servicing Nuclear Rock- 
ets, by Sidney G. Rumbold, ARS Pre- 
print 758-59, March 1959, 9 pp. 

On an Atomic Airplane, by V. Artamkin, 
Rand Corp., T-106, Oct. 14, 1958, 7 pp. 

Survey of Reacting Mixtures Employ- 
ing and U?3 for Fuel and H.O, 
D.O, C, Be, and BeO for Moderator, 
by G. Safonov, Rand Corp., R-259, Jan. 
1954, 13 pp. 


Expanding the Frontiers 
of Space Technology in 


FLIGHT TEST 


W@ Flight testing is conducted by 
Lockheed Missiles and Space 
Division in a unique manner. All 


. components and sub-systems of 


a new project are initially tested 
on known-performance, produc- 
tion missiles. Thus, when the 
final system is ready for first 
flight, its individual components 
already possess flight-tested reli- 
ability. 

This new concept of flight test- 
ing is of major significance and 
has enabled Lockheed to pro- 
duce extremely complex missile 
systems in record time and at 
greatly reduced expense. 

Actual testing is conducted at 
Cape Canaveral, Florida; Van- 
denberg AFB, California, and 
Alamogordo, New Mexico. 
Underwater launch tests for the 
Navy POLARIS FBM —includ- 
ing studies of cavitation, wave 
simulation and skip motion—are 
carried on at the Sunnyvale facil- 
ity and at the Navy test base 
on San Clemente Island. In 
addition, structural and certain 
restricted flight tests are per- 
formed at Hunter's Point Naval 


Shipyard, San Francisco. 


ENGINEERS AND 
SCIENTISTS 


Lockheed Missiles and Space 
Division programs reach far into 
the future and deal with un- 
known and challenging environ- 
ments. If you are experienced in 
one or more of the above areas, 
or have background in related 
work, we invite you to share in 
the future of a company that has 
an outstanding record of achieve- 
ment and make an-important 
individual contribution to your 
nation’s progress in space tech- 
nology. Write::Research and 
Development staff, Dept..I-25, 
962 W. El Camino Real, Sunny- 
vale, California. U.S. citizenship 
required. 


Lochheed 


MISSILES AND SPACE 
DIVISION 


Systems Manager for the 
Navy POLARIS FBM; 
DISCOVERER, SENTRY 
and MIDAS; Army KINGFISHER; 
Air Force Q-5 and X-7 


SUNNYVALE, PALO ALTO, VAN NUYS, 
SANTA CRUZ, SANTA MARIA, CALIFORNIA 
CAPE CANAVERAL, FLORIDA 
ALAMOGORDO, NEW MEXICO © HAWAII 
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ADVANCED 
PRELIMINARY 


DESIGN 
ENGINEERS 


The creation of the new Preliminary 
Design Department of our Solid Rocket 
Plant has led to several unusual and chal- 
lenging positions f for or experienced engineers 


and scientists. 


We are expanding: 


into new fields of rocket propulsion and 
space technology, far beyond “‘classical” 
solid propellant rocketry. For this we are an 
internationally acknowledged leader. 
We are looking for; 
mature, experienced and highly versatile 
engineers, preferably over 30 years’ old. 
These men will be mechanical, aeronautical, 
or missile engineers. They may also be 
engineering-minded physicists or chemists. 
An MS. degree is essential, a Ph.D. wel- 
come. A high degree of proficiency in such 
fields as thermodynamics, aerodynamics, heat 
transfer, stresses, and physics is required, as 
well as a practical understanding of manu- 
facturing fundamentals. Each of the men we 
are seeking must be capable of integrating 
the essentials of a new design, based on 
the more quantitative work of our several 
analytical groups. 


Eight positions as “Technical Specialist’ 
in our Preliminary Design Department are 
available. Each successful candidate will 
work on this team as an equal among equals. 
These top engineers will be given enough 
responsibility and freedom to work as they 
see fit. Excellent salaries are offered, com- 
mensurate with demonstrated ability and 
experience. The Sacramento, California, 
area offers a favorable, healthy climate and 
living conditions which ere among the 
finest in the country, one and a half hours 
driving time from either the Sierra Nevada 
Mountains or San Francisco. 


Please send your resume and direct any 
detailed questions to: 


Mr. Emil L. Eckstein 
Head, Department for Preliminary 
Design 


Through: 

Mr. E. P. James, Supervisor 
Technical and Scientific Placement 
AEROJET-GENERAL CORPORATION 


Box 1947E 


Sacramento, California | 


AEROJET-GENERAL 
CORPORATION 


a subsidiary of 
the General Tire & Rubber Company 


Azusa & near Sacramento, California 


Fundamental Study of Atomic Bat- 
tery, by Chiyoe Yamanaka, Hiroshi 
Wada and Yutaka Yamamura, Osaka 
Univ., Tech. Reps., Vol. 8, nos. 309-336, 
Oct. 1958, pp. 233-241. 

Effects of Chemical Binding on the 
Neutron Cross Section of Hydrogen, 
by W. L. Whiteemore and A. W. Mc- 
Reynolds, Phys. Rev., vol. 113, Feb. 1, 
1959, pp. 806-808. 

Nuclear Rocket Effort Acceleration 
Urged, Aviation Week, vol. 70, March 
23, 1959, pp. 28-29. 

A 20,000-Kilowatt Nuclear Turboelec- 
tric Power Supply for Manned Space 
Vehicles, by Robert E. English, Henry O. 
Slone, Daniel T. Bernatowicz, Elmer H. 
Davison and Seymour Lieblein, NASA 
Memo 2-20-59E, March 1959, 65 pp. 

Some Design Parameters of a Simpli- 
fied Ion Rocket, by H. Preston-Thomas, 
J. Brit. Interplan. Soc., vol. 16, Nov.-Dec. 
1958, pp. 575-583. 

Solid Nuclear Materials as Constit- 
uents of Propulsive Systems, by Jean 
A. Ternisien, Fusées et Recherche Aéron., 
vol. 3, no. 3, 1958, pp. 127-142. 

An Investigation of Barium as an 
Ionic Fuel, by Robert M. Norem, [AS 
Student Branch Paper Competition, 1968, 
First Award Papers, Inst. Aeron. Sci., 
Minta Martin Aeron. Student Fund, 1958, 
127 pp., pp. 48-53. 

First Airborne Atomic Engines Will 
Be Nuclear-Chemical Hybrids, by G. N. 
Nesterenko, A. I. Sobolev and Yu. N. 
Sushkov, Space/Aeron., vol. 31, April 
1959, pp. 47-49. 
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Russian Technical | 


Articles* 


On the Thickness of the Turbulent 
Mixing Zone at the Boundary Layer of 
Two Gas Streams Having Different 
Speeds and Densities, by O. V. Iakovlev- 
skii, Akademiia Nauk SSSR Izvestiia, 
Otdelenie Tekhnicheskikh Nauk, no. 10, 
Oct. 1958, pp. 153-155. (In Russian.) 

Investigation of the Motion of a Vis- 
cous Fluid with Heat Transfer in the 
Domain of Greater Prandtl Numbers, 
by S. I. Kosterin and A. D. Magomelov, 
Akademiia Nauk SSSR Izvestiia, Otde- 
lenie Tekhnicheskikh Nauk, no. 10, Oct. 
1958, pp. 155-157. (In Russian.) 

Thermal Waves in a Layer and a Plate 
Resting on It, by E. I. Chernigovskaia, 
Akademiia Nauk SSSR Izvestiia, Ot- 
delenie Tekhnicheskikh Nauk, no. 10, 
Oct. 1958, pp. 91-93. (In Russian.) 

Convective Heat Transfer in the Non- 
Stationary Regime, by E. A. Sidorov, 
Akademiia Nauk SSSR Izvestiia, 
delenie Tekhnicheskikh Nauk, no. 10, 
Oct. 1958, pp. 116-117. (In Russian.) 

Gas Flow in a Cylindrical Tube with 
Friction and Heat Transfer, by I. I. 
Mezhirov, Akademiia Nauk SSSR Izvestita, 
Otdelenie Tekhnicheskikh Nauk, no. 10, 
Oct. 1958, pp. 118-120. (In Russian.) 

Swirling Flow in Compressible Fluid 
in Circular Pipes, by M. A. Gol’dshtik, 


*The Editors of Technical Literature 
Digest are making a systematic search for 
pertinent Russian articles, both in the orig- 
inal Russian and in translated form. For 
the balance of 1959, these will appear in a 
separate section with this heading, mainly 
to draw the readers’ attention to the — 
listing. After that, they will ap ite 
by item under the proper subject bas. 


ACE TECHNO 


THIS IS THE SYMBOL 


of a Wiley program which deals 
in detail with the many problems 
of outer space 


THIS BOOK IS THE KEYSTONE 


SPACE 
TECHNOLOGY 


Edited by HOWARD S. SEIFERT, 
Space Technology Laboratories, 
Inc., with 38 contributors 


For the first time in one collection, 
all phases of space technology receive 
serious analytical attention. Ma- 
terial is grouped into five major areas: 
1) ballistics and flight dynamics; 2) 
3) communications and 
guidance; 4) man in space; 5) sci- 
entific uses of space. The 500 ref- 
erences form a valuable guide to im- 
portant parts of a vast literature 
1959. 1156 pages. $22.50 
SEND NOW FOR AN ON-APPROVAL COPY. 


JOHN WILEY & SONS, Inc. 


440-FOURTH AVE., NEW YORK 16, N. Y. 


Akademii Nauk SSSR, Izvestiia Otdelenie 
Tekhnicheskikh Nauk, no. 12, Dec. 1958, 
pp. 24-31. (In Russian.) 

Kinetic Theory of Magnetohydrodynam- 
ics, by K. N. Stepanov, Soviet Physics- 
JETP, vol. 34 (7), no. 5, Nov. 1958, 
pp. 892-899. 

On the Theory of the Stabili - 
Liquid Jets in an Electric Field, by 
Glonti, Soviet Physics-JETP, vol. s(n), 
no. 5, Nov. 1958, pp. 917-923. 

On Certain Unsteady Motions of a 
Compressible Fluid, by L. I. Sedov, 
Rand Corp. Translation T-57 (from 
Prikladnaya Matematika i Mekhanika, 
vol. 9, no. 4, 1945, pp. 293-311), Feb. 
1956, 36 pp. 

Calculation of Axisymmetrical Jets 
at Minimum Gravity, by L. E. Sternin, 
Akademii Nauk SSSR, Izvestiia Otdelenie 
Tekhnicheskikh Nauk, Mashinostroenie 1 
Mekhanik, no. 1, 1959, pp. 41-45. (In 
Russian.) 

The Differential Form of the Kinetic 
Equation of a Plasma for the Case of 
Coulomb Collisions, by B. A. Trubnikov, 
Soviet Physics-JETP, vol. 34(7) no. 5, 
pp. 926-931. 

Propagation of Detonation Waves in 
the Presence of a Magnetic Field, by 
E. Larish and I. Shektman, Soviet Physics- 
JETP, vol. 8(35), no. 1, Jan. 1959, pp. 
139-144. 

Ionization of Molecular Hydrogen by 
Ht, H+, and H*+; Ions, by V. V. Afrosi- 
mov, R. F. Il’in and N. V. Fedorenko, 
Soviet Physics-JETP, vol. 34(7), no. 6, 
Dec. 1958, pp. 968-972. 

Method of Perturbations for Problems 
of Strong Detonation, by E. I. Andriankin, 
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Akademiia Nauk SSSR, Izvestiia, Ot- 
delenie Tekhnicheskikh Nauk, no. 12, 
Dec. 1958, pp. 5-14. 

Development of the Process of Com- 
bustion in Propellant Grains in Non- 
stationary Hydrodynamic and Thermal 
Conditions, by V. A. Popov, Akademiia 
Nauk SSSR, Izvestiia Otdelenie Tekhnich- 
eskikh Nauk, no. 12, Dec. 1958, pp. 90- 
95. (In Russian.) 

Regular Precession in a Gyroscope of 
Variable Mass, by V. S. Novoselov, 
Akademii Nauk SSSR, Izvestiia Otdelenie 
Tekhnicheskikh Nauk, no. 11, 1958, pp. 
98-99. (In Russian.) 

Pulse Multifrequency Telemetering De- 
vice, by F. A. Katkov, Abtomatika 7 
Telemekhanika, vol. 20, Jan. 1959, pp. 
54-61. (In Russian.) 

On the Application of Audible Signal 
to Telemetering Apparatus, 
A. Tsykhan, Akademii Nauk SSS 
Izve Otdelenie Tekhnicheskikh Nauk, 
no. 9, 1958, pp. 124-126. (In Russian.) 

Efficiency of Using a Frequency Band 
in Telemetering, by R. R. Vasiliev, Ab- 
tomatika t Telemekhanika, vol. 19, no. 11, 
1958, pp. 1066-1069. (In Russian.) 

Approximate Investigation of Transient 
Process in Second Order Non-linear 
Systems, by E. I. Khlypago. Akademii 
Nauk SSSR, Izvestiia Otdelenie Tekhnich- 
eskikh Nauk, no. 10, 1958, pp. 1-11. 
(In Russian.) 

Choice of Power Unit of Optimal Auto- 
matic Control Systems, by L. N. Fitsner, 
Abtomatika i Telemekhanika, vol. 19, 
Dec. 1958, pp. 1107-1117. (In Russian.) 

Determination of Optimal Pulse Trans- 
fer Function of a Servo-system for a 
Certain Class of Disturbances, by P. S. 
Matseev, Abtomatika i Telemakhanika, 
vol. 20, Jan. 1959, pp. 3-15. (In Rus- 
sian.) 

On Certain Oscillations in Pulse Sys- 
tems, by Da-Chuan Shao, Abiomatika 7 
Telemekhanika, vol. 20, Jan. 1959, pp. 
85-89. (In Russian.) 


Approximate Determination of Self- 
oscillations in a Synchro-motor Control 
System, by D. P. Petelin. Abtomatika ¢ 
Telemekhanika, vol. 20, Jan. 1959, pp. 
16-22. (In Russian.) 


Self-oscillations in a Control Single- 
loop System Having Two Symmetric 
Relays, by Syui-jan Tu and Lui-vy Tei, 
Abtomatika 1 Telemekhanika, vol. 20, 
Jan. 1959, pp. 90-94. (In Russian.) 


Optical Tracking of Soviet Artificial 
Satellites, by E. Z. Gindin et al., Russian- 
English Technical Translation Center, 
Hollywood 28, Calif., Translation no. 
Al, 1959, 70 pp. $2.65. 

USSR Biomedical Investigations in 
Outer Space, by A. M. Calkin et al, 
Russian-English Technical Translation 
Center, Hollywood 28, Calif., Translation 
no. A8, 1959, 80 pp. $2.90. 

From Sputniks to Interplanetary Flights, 
by A. Sternfeld, Russian-English Tech- 
nical Translation Center, Hollywood 28, 
te Translation no. B, 1959, 176 pp. 
5.50. 


Flight of Rockets into Interplanetary 
Space, by N. A. Merkulow, Russian- 
English Technical Translation Center, 
Hollywood 28, Calif., Translation no. C, 
1959, 125 pp. $4.85 


Results of Investigations of the Upper 
Ionosphere Based on Radio Signal Ob- 
servations from the First USSR Earth 
Satellites, by A. L. Alpert et al., Russian- 
Iinglish Technical Translation Center, 
Hollywood 28, Calif., Translation ‘no. A2, 
1959, 90 pp. ‘$2.80. 
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...and TOTAL systems 
in which the big bird and support 
equipment may rank only 
as a component. 


This difference between systems can make 
v a big difference in your career 


IF YOU ARE QUALIFIED and interested in contributing to programs of ; 
“total” scope, it will be of value to you to investigate current oppor- 
tunities with General Electric’s DEFENSE SYSTEMS DEPT., whose 
work lies primarily in providing total solutions to large scale defense 
problems of the next 5, 10 and 20 years. 


The work here lies almost entirely in the areas of systems engi- 
neering and systems management. 


Inquire about these positions: 


Guidance Equation Engineers Systems Test Evaluation Engineers 

Systems Logistics Engineers Engineering Psychologists 

Electronic Systems Radar Equipment Engineers 
Management Engineers Weapons Analysis Engineers 

Operations Analysis Engineers Weapons Systems Integration 

Systems Program Engineers Engineers 

Data Processing Engineers Engineering Writers 


Forward your confidential resume at an early 

Whereas the growth potential is evident — both for DSD and the ; 

yy)y)))) engineers who join us — the positions we fill during these early — 
months will carry significant “ground-floor” benefits. 


ALN Write fully to Mr. E. A. Smith, Room 9-M : 


DEFENSE SYSTEMS DEPARTMENT 


A Department of the Defense Electronics Division 


GENERAL ELECTRIC. 


300 South Geddes Street, Syracuse, N. Y. 
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SOUTHWEST 


ROD END 
TYPES 


PLAIN TYPES 


PATENTED U. S.A. 
World Rights Reserved 


CHARACTERISTICS 
ANALYSIS RECOMMENDED USE 


: For types operating under high 
Stainless Steel Ball and Race pera ant (800-1200 degrees F.). 
Chrome Alloy Steel Ball For types operating under yt radial 
and Race ultimate loads (3000-893,000 Ibs.). 


Bronze Race and Chrome For types operating under normal loads 
Steel Ball with friction requi 


Thousands in use. Backed by years of service life. Wide variety of 
Plain Types in bore sizes 3/16” to 6” Dia. Rod end types in similar 
size range with externally or internally threaded shanks. Our Engi- 
neers welcome an opportunity of studying individual requirements 
and prescribing a type or types which will serve under your demand- 
ing conditions. Southwest can design special types to fit individual 
specifications. As a result of thorough study of different operating 
conditions, various steel alloys have been used to meet specific 
needs. Write for Engineering Manual No. 551.Address Dept. ARS-59 


SOUTHWEST PRODUCTS CO. 


1705 SO. MOUNTAIN AVE., MONROVIA, CALIFORNIA 


RESEARCH ENGINEERS 


Heat and Mass Transfer 


Outstanding opportunities exist in our Heat 
and Mass Transfer Section for men with ad- 
vanced degrees and experience in this field. 
Duties will involve working on research pro- 
grams in Aerodynamic Heating, Thermal 
Radiation, Heat Transfer in Reactor Design, 
and Solar Energy. 


These positions are non-routine and require 
men with initiative and creative ability. Ex- 
cellent employee benefits and tuition free 
graduate study. Please send resume to: 


E. P. Bloch 


ARMOUR RESEARCH FOUNDATION 


10 West 35th Street 
Chicago 16, Illinois 
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